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OPTICAL  SIGNAL  PROCESSING 


1.0  INTRODUCTION 

As  the  bandwidth  of  signals  increase  and  as  the  electromagnetic 
environment  becomes  increasingly  dense,  processing  operations  such  as 
convolution,  spectrum  analysis,  correlation,  ambiguity  function  generation 
and  filtering  become  computationally  intensive  operations.  Optical  systems 
have  the  capability  to  perform  a  large  number  of  complex  multiplications  and 
additions  per  unit  time.  As  the  bandwidth  increases,  the  number  of 
computations  increase  as  the  square  of  the  bandwidth  because  the  degrees  of 
freedom  of  the  signal  is  linearly  proportional  to  bandwidth  and  the 
allowable  computational  time  interval  is  a  inversely  proportional  to 
bandwidth.  Optical  processing  provides  high-speed,  parallel  computations  so 
that  digital  post-processing  techniques  can  be  used  for  lower-speed,  serial 
computation. 

The  research  effort  described  in  this  report  has  resulted  in  several 
innovative  optical  processing  techniques  for  improved  performance;  it  covers 
the  three-year  period  from  1  October  1980  to  30  September  1983.  The  major 
accomplishments  can  be  divided  into  three  areas:  (1)  interferometric 
spectrum  analysers,  (2)  acoustic  spreading  in  both  single  and  multichannel 
Bragg  cells,  and  (3)  adaptive  optical  processing  based  on  transversal 
filtering  with  feedback.  In  the  following  paragraphs,  we  summarize  the  key 
results  in  each  area;  further  details  can  be  found  in  the  referenced  journal 
articles  that  have  been  published. 


2.0  INTERFEROMETRIC  SPECTRUM  ANALYZERS 


A  conventional  Bragg  cell  power  spectrum  analyzer  that  measures  the 
instantaneous  frequency  content  of  a  wideband  signal  generally  has  a  dynamic 
range  of  25-30  dB.  We  developed  an  interferometric  spectrum  analyzer 
technique  that  uses  a  unique  reference  wavefront  which  serves  as  a 
distributed  local  oscillator.  As  a  result,  the  photodetectors  measure  the 
instantaneous  amplitude  of  the  frequency  content  rather  than  the 
instantaneous  power.  The  dynamic  range  is  thereby  doubled  in  dB  so  that, 
for  an  equivalent  set  of  system  parameters,  the  dynamic  range  is  50-60  dB. 

In  addition  to  providing  a  significant  increase  in  the  dynamic  range, 
the  distributed  local  oscillator  provides  a  fixed  offset  temporal  frequency 
at  each  photodetector  position.  The  post-detection  circuitry  for  each 
photodetector  element  is  therefore  identical;  furthermore,  both  the 
amplitude  and  phase  of  the  instantaneous  spectrum  can  be  measured  if 
desired.  Because  the  scattered  light  is  not  frequency  shifted,  the  system 
is  much  less  sensitive  to  scattered  light.  Another  advantage  of  this 
technique  is  that  short  pulses,  such  as  those  from  radar  systems,  can  be 
detected  even  though  their  duration  is  less  than  the  Bragg  cell  transit 
time. 


Further  details  are  provided  in  Reference  1  in  which  the  results  of  the 
analytical  effort  was  published.  Since  that  time,  Harris  has  completed  a 
feasibility  study  and  prototype  development  effort  for  AFAL.  The 
experimental  results  are  in  excellent  agreement  with  the  theory. 


3.0  ACOUSTIC  SPREADING  IN  BRAGG  CELLS 


An  acousto-optic  device  operated  at  high  RP  frequencies  is  often  called 
a  Bragg  cell.  A  piezoelectric  transducer  is  bonded  to  a  suitable 
interaction  medium  to  convert  an  electrical  into  a  traveling  acoustic 
pressure  wave.  This  pressure  vave,  in  turn,  causes  an  index  of  refraction 
change  that  modulates  the  light  in  space  and  time.  The  acoustic  wave,  as  it 
propagates  away  from  the  transducer,  spreads  in  much  the  same  way  as  does 
light  from  a  small  source.  We  analyzed  the  effects  that  this  acoustic 
spreading  has  on  some  optical  processing  operations. 

For  a  single  transducer  Bragg  cell,  the  principal  effort  of  acoustic 
spreading  is  a  curvature  of  the  Fourier  spectral  components.  The  locus  of 
this  curve  is  a  parabola  whose  vertex  is  at  the  position  corresponding  to 
the  RF  frequency  and  whose  focus  is  at  the  optical  axis.  The  degree  of 
curvature  is  dependent  on  the  anisotropic  nature  of  the  interaction  medium; 
a  special  case  is  an  isotropic  medium  for  which  the  parabola  degenerates 
into  a  circular  function. 

The  effects  of  this  curvature  on  spectrum  analysis  is  a  function  of  the 
size  of  the  photodetector  in  a  direction  orthogonal  to  that  of  the  acoustic 
propagation.  For  a  photodetector  haling  a  very  large  size  (one  which  would 
be  acceptable  if  the  curvature  were  cot  present)  the  effect  is  to  cause  the 
peak  value  of  the  diffraction  pattern  to  move  toward  the  optical  axis  by  4Z 
of  its  main  lobe  width,  to  increase  the  main  lobe  half  power  points  by  3Z, 


and  to  increase  the  side lobe  levels  by  5  dB  on  the  side  of  the  main  lobe 
toward  the  optical  axis.  Analysis  shows  that  these  effects  can  be 
controlled  by  reducing  the  sire  of  the  photodetector  with  a  penalty  in  the 
light  collecting  power  of  about  10Z.  Further  details  of  this  analysis  and 
experimental  results  are  given  in  a  paper  which  was  published  in  Applied 
Optics  (Reference  2). 

The  acoustic  spreading  manifests  another  problem  in  Bragg  cells  having 
multiple  transducers.  The  phenomena  is  that  the  spreading  acoustic  wave  from 
a  given  transducer  overlap  with  waves  from  adjacent  transducers.  An 
analysis  of  this  phenomena  showed  that  the  optical  effects  of  acoustic 
spreading  can  be  compensated  by  a  holographic  optical  element  placed  in  a 
combination  image/Fourier  plane.  At  an  image  plane  of  the  Bragg  cell,  we 
find  that  the  corrected  light  distribution  is  constrained  to  lie  within 
channels  whose  height  is  nearly  equal  to  that  of  the  transducer.  He 
desionatrated  that,  for  a  12-channel  Bragg  cell  operating  over  an  80  MHz 
bandwidth,  the  usable  aperture  (defined  as  that  region  for  which  the  channel 
information  does  not  overlap)  could  be  increased  from  10Z  to  100Z. 
Furthermore,  the  channel  packing  density  can  be  increased  by  a  factor  of 
four.  Details  of  this  analysis  and  the  supporting  experiments  will  be 
published  in  the  December  1983  issue  of  Applied  Optics  (Reference  3). 

4.0  ADAPTIVE  OPTICAL  PROCESSING 

The  transversal  filter  is  widely  used  in  digital  data  processing.  An 
even  wider  range  of  applications  are  possible  when  feedback  loops  are  used; 
we  have  developed  a  method  for  implementing  such  operations  optically.  The 
basic  concept  is  to  consider  a  Bragg  cell  as  being  equivalent  to  a  delay 
line  which  can  be  tapped  optically  instead  of  electrically.  In  this 
fashion,  very  wide  bandwidth  signals  can  be  processed. 


Fundamental  to  the  notion  of  transversal  filtering  is  the  need  to 
produce  the  appropriate  tap  weights.  An  innovative  way  to  produce  them 
optically  is  to  use  two  orthogonally  oriented  Bragg  cells  arranged  so  that 
they  are  mutually  imaged  at  a  given  plane.  Ve  show  that  if  the  light  in 
this  plane  is  integrated  along  one  of  the  diagonals,  the  tap  weights  are 
produced.  A  convenient  way  to  perform  the  integration  is  to  create  the  two- 
dimensional  Fourier  transform  of  the  tap  weight  image  plane  and  to  evaluate 
the  transform  along  a  line  normal  to  the  diagonal  over  which  the  integration 
is  to  be  made. 

The  processing  is  completed  by  interferometrically  adding  the  Fourier 
transform  of  the  signal  to  be  processed  to  the  Fourier  transform  of  the  tap 
weights.  This  total  light  distribution  is  then  square-law  detected  by  a 
single  element  photodetector  to  produce  an  estimate  of  the  received  signal. 
Depending  on  the  application,  the  estimate  is  further  processed 
electronically  and  the  difference  between  it  and  the  received  signal 
provides  the  feedback  signal.  Thus ,  the  system  adapts  to  a  changing  signal 
environment . 

The  chief  advantage  of  this  optical  processing  technique  is  that 
transversal  filtering  concepts  can  now  be  applied  to  signals  having 
bandwidth8  in  the  SO  MHz  to  500  MHz  range.  The  system  has  an  excellent 
dynamic  range  because  it  is  interferometric  in  nature  and  is  linear  in  light 
amplitudes.  The  tap  weights  contain  both  amplitude  and  phase  information 
which  provides  for  increased  computational  accuracy.  Scattered  light  from 
the  undiffracted  beam  doe#  not  affect  the  performance  of  the  system.  Details 
of  this  concept  are  contained  in  a  paper  published  in  Applied  Optics 
(Reference  4). 

The  finite  length  of  the  Bragg  cells  impose  a  constraint  on  the  number 
of  signal  components  that  are  integrated  to  produce  the  tap  weights.  We 
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performed  several  computer  simulations  of  the  adaptive  filtering  system  to 
characterize  and  determine  the  performance  due  to  this  limitation  (digital 
systems  can  integrate  signal  contributions  for  an  indefinite  period).  The 
application  ve  chose  to  illustrate  these  effects  is  that  of  notch  filtering 
or  frequency  excision.  Ve  found  that  ve  can  improve  the  performance  of  the 
system  significantly  if  the  tap  veight  plane  is  tapered  so  that  the 
readaptation  phenomena  is  minimized  as  signal  components  leave  the  tap 
veight  plane  (the  accumulator).  In  principle,  the  performance  could  be 
improved  if  the  taper  vere  extended  to  also  veight  the  signal  components  as 
they  enter  the  accumulator;  this  vas  not  tested  in  the  simulations.  The 
taper  is  equivalent  to  using  a  "leaky"  integrator  in  conventional  analog 
systems. 

The  simulations  shov  that  the  performance  improves  as  the  number  of  tap 
veights  increase;  this  argues  for  using  Bragg  cells  having  a  large  time- 
bandvidth  product.  We  also  found  that  the  system  adapts  rapidly  so  that 
notches  can  be  quickly  formed  to  track  agile  jammers.  Notch  depths  of  the 
order  of  30-40  dB  vere  obtained,  depending  on  the  t ime-bandvid th  product  and 
the  number  of  taps  used.  Since  the  t ime-bandvid th  product  of  Bragg  cells 
are  exceptionally  high,  a  much  larger  number  of  tap  veights  are  available 
as  compared  to  digital  processing  systems.  The  results  shov  that  the  finite 
length  of  the  accumulators  is  not  a  serious  dravback  for  many  applications, 
particularly  those  in  vhich  the  system  must  operate  in  a  rapidly  changing 
signal  environment.  The  performance  of  the  system  then  approaches  that  of 
one  having  an  infinite  accumulator  and  the  feedback  gain  value  needed  to 
give  equivalent  tracking  performance.  Details  of  these  simulations  are 
given  in  a  paper  that  has  been  accepted  for  publication  by  Optical 
Engineering  (Reference  5).  A  paper,  summarizing  these  results  vas  given  at 
the  SPIE  Conference  on  Advances  in  Optical  Information  Processing  in  January 
1983  at  Los  Angeles.  A  short  paper  also  appears  in  the  conference 
proceedings  (Reference  6). 
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Interferometric  spectrum  analyzer 

A.  Vander  Lugt 


Dynamic  range  is  a  key  performance  parameter  for  spectrum  analyzers.  The  dynamic  range  of  a  Bragg  cell 
power  spectrum  analyzer  is  generally  limited  by  the  dynamic  range  of  self-scanned  photodetector  arrays. 
Interferometric  techniques  can  be  used  to  increase  the  dynamic  range:  but  it  is  at  the  expense  of  increasing 
the  number  of  photodetectors  required,  when  the  interference  is  introduced  in  the  spatial  domain,  or  a  Large 
photodetector  bandwidth,  when  the  interference  is  introduced  in  the  temporal  domain  In  this  paper  we 
describe  an  interferometric  approach  wherein  a  second  Bragg  cell  generates  a  spatially  modulated  reference 
waveform  to  produce  an  interference  term  that  has  a  constant  temporal  frequency  for  all  spatial  frequencies. 
The  advantages  of  this  approach  are  lower  photodetector  bandwidth,  improved  dynamic  range,  improved 
cross  talk  suppression,  more  efficient  use  of  the  Bragg  cell  time-bandwidth  product,  immunity  to  scattered 
noise,  and  improved  short  pulse  detectability.  The  chief  disadvantage  is  the  need  for  a  discrete  element 
photodetector  array*,  when  such  arrays  become  available  in  hybrid  or  integrated  packages,  an  additional  ad¬ 
vantage  will  be  that  of  parallel  postdetection  processing. 


I.  Introduction 

Optical  processing  techniques  can  be  used  to  produce 
the  instantaneous  spectrum  of  wide  bandwidth  signals. 
Lambert1  described  how  acoustooptic  Bragg  cells  can 
be  used  to  convert  an  electrical  waveform  /(£  i  to  an 
optical  waveform  that  is  a  function  of  both  space  and 
time.  If  the  cell  is  coherently  illuminated,  the  optical 
system  displays  the  Fourier  transform  of  that  segment 
of  the  waveform  present  within  the  cell.  The  complex 
valued  Fourier  transform  can  be  described  as 

Xi-r 

f(u)  expl-ywul  du.  (1) 

where  T  is  the  processing  time  of  the  cell  and  is  a 
temporal  radian  frequency.  Generally,  Fr(«,£ )  is  called 
the  instantaneous  spectrum  of  the  input  waveform  fit); 
we  note  that  it  is  a  function  of  the  present  time  £  as  well 
as  a  segment  of  its  past  history. 

A  photodetector  array  at  the  Fourier  plane  senses  the 
instantaneous  energy  spectrum  | Fji  w.£ )  | 2;  such  a  sys¬ 
tem  is  generally  called  a  power  spectrum  analyzer.  A 
key  performance  parameter  of  any  spectrum  analyzer 
is  the  achievable  dynamic  range  which  is  a  function  of 
the  laser  power,  the  dynamic  ranges  available  from  the 
Bragg  cell  and  the  photodetector  array,  the  time- 


The  author  is  with  Hams  Corporation.  Advanced  Technology 
Department.  Government  systems  Group.  P  0.  Box  17.  Melbourne. 
Florida  :29<H. 

Received  t*  rebruarv  1981. 

HJ03-89:' >  81  162770- IOSOO..W 0. 

•:  .981  Omicai  Society  oi  America. 


bandwidth  product,  and  the  overall  efficiency  of  the 
system.  Power  spectrum  analyzers  that  use  self- 
scanned  photodetector  arrays  are  generally  limited  to 
a  25-35-dB  dynamic  range  due  to  the  squaring  opera¬ 
tion  on  F 7*(u).£  )  and  the  inherent  dynamic  range  limi¬ 
tation  of  the  array.  As  the  photodetector  dynamic 
range  improves,  the  laser  power  is  the  next  important 
limitation  to  achieving  a  large  dynamic  range  since 
Bragg  cells  have  been  developed  having  very  large  dy¬ 
namic  ranges. 

One  method  for  increasing  the  dynamic  range  is  to  use 
an  interferometric  spectrum  analyzer  whose  output  is 
proportional  to  the  instantaneous  magnitude  spectrum 
|Friw.£)|.  King  et  ai-  describe  heterodyning  tech¬ 
niques  for  recovering  both  the  amplitude  and  phase 
information  in  a  light  distribution.  In  their  system  the 
interference  of  an  unmodulated  reference  beam  with 
Fr((i>,£)  produces  a  temporal  frequency  proportional  to 
the  input  signal  frequency  u>.  Since  the  fractional 
bandwidth  of  the  input  signal  is  usually  ~50%.  centered 
on  a  frequency  of  several  hundred  megahertz,  the  in¬ 
terference  term  occurs  at  a  rather  high  frequency  that 
varies  as  a  function  of  the  spatial  frequency. 

In  this  paper.  I  describe  an  interferometric  spectrum 
analyzer  in  which  a  spatially  modulated  reference  beam 
is  used  to  reduce  the  temporal  interference  frequency 
to  a  small  and  fixed  value  over  the  entire  spectrum. 
Discrete  element  photodetectors  having  a  small  band¬ 
width  and  a  low  noise  equivalent  power  can  then  be 
used:  an  additional  benefit  of  discrete  detectors  is  that 
the  postprocessing  operations  are  more  flexible  and  can 
be  performed  in  parallel  to  reduce  the  output  data 
rate. 
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Fig.  1.  Interferometric  spectrum  analyzer. 

Self-scanned  arrays  cannot  be  used  with  this  tech¬ 
nique  because  they  integrate  over  a  time  period  Large 
compared  with  2 ir/w.  and  the  information  contained  in 
the  interferometric  output  is  lost.  Turpin3  and  Bader4 
desc'ibe  interferometric  spectrum  analyzers  which 
produce  a  spatial  fringe  structure  instead  of  a  temporal 
fringe  structure;  self-scanned  arrays  can  be  used  in  such 
systems,  but  the  number  of  elements  must  be  increased 
to  resolve  the  spatial  fringes.  Although  discrete  pho¬ 
todetectors  are  not  the  most  elegant  for  use  in  systems 
having  a  large  time-bandwidth  product,  fairly  large 
arrays  have  been  implemented,5  and  advanced  photo¬ 
detector  fabrication  techniques8  may  produce  inte¬ 
grated  devices  having  attractive  operational  features. 

Some  additional  advantages  of  interferometric 
spectrum  analysis  are  improved  cross  talk  rejection, 
immunity  to  scatter  noise,  short  pulse  detectability,  and 
uninterrupted  evaluation  of  the  spectrum.  In  Sec.  II 
we  describe  the  basic  theory  of  the  interferometric 
spectrum  analyzer  and  establish  the  required  charac¬ 
teristics  of  the  reference  beam.  In  Sec.  QI  we  determine 
the  photodetector  geometry  and  postdetection  band¬ 
width  required  to  achieve  a  given  frequency  resolution. 
In  Sec.  IV  we  analyze  and  compare  the  performance  of 
candidate  reference- beam  waveforms.  In  Sec.  V  we 
compare  the  laser  power  required  and  the  dynamic 
range  obtained  by  this  interferometric  method  with 
those  of  a  power  spectrum  analyzer. 

II.  Basic  Theory 

Consider  the  interferometric  system  shown  in  Fig.  I; 
this  system  does  not  suggest  how  a  practical  spectrum 
analyzer  would  be  confijcured  but.  is  used  to  explain  the 
theory  of  operation.  The  two  Bragg  cells  are  illumi¬ 
nated  by  a  collimated  source  of  monochromatic  light  at 
the  Bragg  angle.  The  signal  waveform  fit)  is  applied 
to  the  tiansducer  of  the  Bragg  cell  located  in  the  lower 
leg  of  the  interferometer  at  plane  P\.  Lens  L-z  produces 
the  Fourier  transform  of  the  complex  light  amplitude 
a\ix,t )  leaving  the  cell:  for  a  low  modulation  index  we 
have  a  series  expansion  for  the  signal  given  by 

atix.:i  «  ,AIU  +  -  x/v)  cos(ai„it  -  x/o )|  +  (2) 

where  m  ( is  the  modulation  index,  u  is  the  velocity  of  the 
acoustic  wave  within  the  cell,  and  is  the  center  fre¬ 
quency  of  the  applied  rf  signal.  The  higher  order  terms 
can  generally  be  neglected  if  the  fractional  bandwidth 
is  <50°r.  Further,  we  are  usually  interested  in  only  the 
positive  diffracted  order  so  that  the  constant  as  well  as 


the  negative  diffracted  order,  which  is  highly  suppressed 
by  the  Bragg  mode  of  operation,  can  be  ignored.  A 
fraction  a  of  the  laser  power  P  is  directed  into  the  signal 
beam  so  that  the  amplitude  factor  Ai  is  equal  to  ( aP / 
L ) 1/2,  where  L  is  the  length  of  the  cell.  The  effects  of 
weighted  illumination  and  optical  losses  will  be  con¬ 
sidered  later. 

The  amplitude  light  distribution  for  the  First  dif¬ 
fracted  order,  if  we  ignore  a  time  delay  T /2,  is 

Xui 

fit  -  x/v)  expl/WcU  -  x/v>]  etp< -jpx )dz. 
■Ln 

(3) 

where  p  is  a  radian  spatial  frequency  variable  related 
to  a  physical  distance  £  in  plane  P>  by  p  *  2 v^/XF,  \  is 
the  wavelength  of  the  source,  and  F  is  the  focal  length 
of  lens  L-z .  Through  a  change  of  variables,  we  can  re¬ 
write  Eq.  (3)  as 

f(u)  exp \Juv(p  -  pc)\du. 

t*rn 

(41 

The  integral  portion  of  Eq.  (4)  is  similar  to  Eq.  ( 1)  ex¬ 
cept  that  it  is  centered  at  a  position  corresponding  to  pc 
in  plane  P>,  and  the  limits  of  integration  are  slightly 
different  because  we  ignored  a  time  delay  equal  to  772 
in  the  representation  given  by  Eq.  ( 2).  The  exponential 
phase  factor  shows  that  the  light  frequency  is  shifted 
linearly  as  a  function  of  p.  Since  pv  ■  ui,  the  light  dif¬ 
fracted  by  any  continuous  waveform  has  the  same 
temporal  frequency  as  the  signal  component. 

Suppose  that  the  input  signal  waveform  has  a  cw 
component  of  frequency  w*  and  amplitude  C*.  The 
light  distribution  in  plane  P»  can  then  be  calculated 
from  either  Eq.  (3)  or  (4c 

.4i(p.i)  a;'mtA|LC*  exp< -)<■>*< ) sinc[(p  —  p»)L/2ir],  (5) 

which  reveals  that  the  spectrum  of  a  sinusoidal  input 
is  a  (sinxx)/Tx  function,  centered  at  p*.  whose  ampli¬ 
tude  is  proportional  to  CV  The  entire  function  (in¬ 
cluding  the  sidelobes)  is  multiplied  by  a  phasor  of  fre¬ 
quency  w*.  Equation  (5)  further  shows  that  the  fre¬ 
quency  resolution  is  2vv/L  if  we  use  the  Rayleigh  cri¬ 
terion  for  resolution.  We  will  use  the  function  Aiip.t) 
as  given  by  Eq.  (5)  extensively  in  subsequent  analyses; 
first,  however,  we  consider  the  spectrum  of  short  pulses 
whose  duration  may  be  less  than  or  equal  to  the  pro¬ 
cessing  time  of  the  cell. 

Consider  the  case  of  a  pulse  having  a  duration  T0  less 
than  or  equal  to  T  and  having  a  carrier  frequency  ai*. 
Let  the  time  of  arrival  be  at  t  *  0:  at  some  later  time  t . 
the  leading  edge  of  the  pulse  will  have  moved  to  a  po¬ 
sition  — VoL  +  vt.  Consider  the  output  for  the  time  in¬ 
terval  0  <  t  <  To  so  that  the  trailing  edge  of  the  pulse 
has  not  yet  entered  the  cell:  we  have  that 

■4iip.fl  »  jm\A\  J  exp(-ru)*ir  -  xii-  -  T  21] 

X  exp<—  ;px  utx.  ' 61 

where  we  have  now  included  the  delay  factor  772  in  the 
integral.  After  a  change  of  variables,  we  find  that 
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.Adp.t)  ■  exp ijpL/2)  exp(-/ip  +  p*  >t'f ,'2| 

X  sinc[(p  -  p*k’£/2x|.  iTl 

This  result  shows  that  the  spectrum  of  the  pulse  is 
centered  at  p*  and  that  as  t  increases,  the  amplitude  of 
the  sine  function  increases  while  its  width  decreases. 
We  are  particularly  interested  in  the  temporal  fre¬ 
quency  behavior  of  this  function.  Note  that  at  p* ,  the 
centroid  of  the  sine  function,  the  temporal  frequency 
is  a;*  which  is  the  frequency  of  the  carrier.  In  contrast 
to  the  case  of  a  cw  signal,  however,  the  sidelobes  of  the 
sine  function  have  a  frequency  (p  +  p*)i7 2ir  which 
varies  continuously  as  p  varies.  In  particular,  note  that 
at  some  spatial  frequency  p*  +  Ap,  the  temporal  fre¬ 
quency  is  hi,  +  Aoj/2. 

The  last  situation  we  consider  is  that  of  a  pulse  whose 
leading  and  trailing  edges  are  both  within  the  aperture 
of  the  Bragg  cell.  Let  £  ■  0  be  the  time  at  which  the 
trailing  edge  of  the  pulse  has  just  entered  the  cell;  then, 
for  0  <  t  <  ( T  -  To),  we  have 

X—  L,  2+ec-*u7*o 
■L/2*ut 

X  exp{-;oi*(t  -  x/v  -  772)! 
x  exp<  —jpx )dx 

m  jm\AivToexp[jp\L  -  t’T0)/2| 
x  expljpki'To/2)  exp(  -jpvt  I 
X  jinc|(p  -  pk)vTJ2r\.  18) 

We  see  that  the  sine  function  has  an  amplitude  and  a 
width  that  are  determined  by  vT^  both  are  independent 
of  time.  As  before,  the  temporal  frequency  at  p*  is  w*. 
The  frequency  of  the  sidelobes  at  some  spatial  fre¬ 
quency  p  =  pfc  +  Ap  is  equal  to  w*  +  Aw. 

We  can  now  summarize  the  three  cases  discussed  so 
far.  In  each  case,  the  temporal  frequency  at  the  cen¬ 
troid  of  the  sine  function  is  w*.  At  some  incremental 
spatial  frequency  Ap  away  from  the  centroid,  we  find 
that  the  temporal  frequency  of  the  sidelobes  differs 
from  that  at  the  centroid  by  zero  for  the  case  of  a  cw 
signal  ior  a  very  long  pulse),  by  Aw/2  for  a  pulse  with 
only  its  leading  edge  in  the  aperture,  or  by  Aw  for  a  pulse 
with  both  leading  and  trailing  edges  within  the  aperture. 
The  explanation  for  this  phenomenon  is  that  in  the  first 
case  neither  edge  defining  the  pulse  moves  in  time  since 
it  is  the  aperture  of  the  Bragg  cell  that  determines  the 
signal  length.  In  the  second  case  one  edge  moves,  but 
the  other  is  stationary,  in  the  third  case,  both  edges 
move  at  the  same  velocity.  The  temporal  frequencies 
in  the  sidelobe  structure,  then,  contain  the  information 
about  where  the  centroid  of  the  pulse  is  located  within 
the  Bragg  cell.  For  example,  if  we  collect  the  terms  in 
p  from  Eq.  (3),  we  have  a  term 

»  exp{— jpi—L.2  ft  +  '-7V2t|.  (9) 

which  shows  that  the  centroid  of  the  pulse  is  at  -Li 2  + 
vt  +  c T.,/2:  the  centroid  moves  at  the  velocity  of  the 
acoustic  wave.  The  residual  phase  term  in  Eq.  1 8>  is 
then  expt;p,fT.)/2)  which  is  a  constant;  the  argument 
PbC  r.„  i  is  equal  to  one  ha if  the  number  of  carrier  cycles 
contained  in  the  pulse.  A  similar  calculation  for  the 
case  of  a  pulse  with  only  the  leading  edge  in  the  aperture 

<*»**-•/* 
tm!  /  d. 


is  obtained  from  Eq.  (7)  which  shows  that  the  centroid 
of  the  pulse,  located  at  —Li 2  +  cr/2.  is  moving  at  one 
half  the  acoustic  velocity.  The  number  of  cycles  in  the 
pulse  is  pkt't/2  which  increases  as  the  effective  pulse 
width  increases.  For  a  very  long  pulse  the  centroid  is 
fixed  at  the  optical  axis  and  the  number  of  cycles  in  the 
pulse  is  fixed  at  p*t)T/2x  -  pkL/2tt. 

We  now  turn  our  attention  to  the  question  of  gener¬ 
ating  a  reference  beam  that  will  produce  an  interfero¬ 
metric  fringe  structure  having  a  fixed  temporal  variation 
at  all  spatial  frequencies.  A  spatially  modulated  ref¬ 
erence  beam  will  produce  a  fixed  offset  frequency  be¬ 
cause  the  light  distribution  at  the  Fourier  plane  will 
have  a  temporal  and  spatial  frequency  relationship 
similar  to  that  of  the  signal  spectrum.  The  reference 
distribution  must,  however,  be  displaced  in  plane  P± 
relative  to  the  signal  spectrum  so  that  the  spatial  dis¬ 
placement  does  not  affect  the  temporal  frequency  dis¬ 
tribution.  The  spatial  displacement  can  be  achieved 
by  rotating  the  combining  beam  splitter  through  a  small 
angle  or  by  placing  a  prism  immediately  after  the  ref¬ 
erence-beam  Bragg  cell.  In  either  case,  the  reference- 
beam  waveform  at  plane  P>  will  be  displaced  to  produce 
a  fixed  frequency  offset  over  the  entire  spectrum  as  we 
shall  now  show. 

We  denote  the  bandwidth  of  the  Bragg  cell  by  W  so 
that  the  time-bandwidth  product  is  TW.  We  begin  by 
choosing  a  reference- beam  waveform  that  is  a  pulse  of 
unit  amplitude  and  length  L/TW.  Such  a  signal  is 
equivalent  to  one  resolution  element  at  the  plane  of  the 
Bragg  cell.  The  Fourier  transform  of  this  narrow  pulse 
can  be  obtained  directly  from  Eq.  (8)  if  we  replace  dT0 
by  L/TW  and  replace  p  by  p  +  pd,  where  pd  represents 
the  small  displacement  between  the  reference  and  signal 
beams  at  plane  P>.  The  reference- beam  distribution 
in  plane  P  <  then  becomes 


A->\p.t ) 
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expO®>  exp(-;<p  +  pa)cf  I 


x  sinc((p  +  pa  -  pc  >Li2xTW\,  1 10) 


where  <t>  is  a  phase  term  that  is  not  a  function  of  time  and 
pc  is  the  same  center  frequency  as  that  used  in  the  signal 
beam  Bragg  cell.  The  envelope  of  the  reference  beam 
is  not.  of  course,  constant  in  p  as  we  desire,  but  is  a  sine 
function  that  is  a  factor  of  TW  broader  than  the  signal 
distribution.  For  the  moment,  we  shall  assume  that  the 
reference  beam  is  constant  over  the  frequency  band. 

The  intensity  at  plane  Pi  is  given  by  the  squared 
magnitude  of  the  sum  of  the  signal  and  reference-beam 
distributions: 


ftp.:)  m  | .A itp.r )  +  Ai'p.t  i|- 


■  11) 


■  I  a  ii  p.t  i|J  +  |a^p.£i|J  +  2  Re(.4|ip..'  i.-t»ip.n|.  1 12) 

where  *  denotes  complex  conjugate.  Each  photode¬ 
tector  will  integrate  light  over  an  interval  p*  -  Ad  to 
pi,  Ap.  where  pi,  denotes  the  center  position  and  2Ap 
the  width  of  the  photodetector.  Let  1  be  one  half  the 
ratio  of  the  width  of  the  photodetector  to  the  spacing 
between  the  photodetectors  so  'hat  Ap  =  2~  i  L.  The 
integrated  intensitv  for  the  fir«t  term  of  Eq.  T2) 
calculated  by  using  Eq.  *  5 1 : 
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Q\m—  |.4iip.M|*dp  ■  2Jmf.4fl.C2,  tl3l 

-T  %/p*--r“ 

which  is  valid  for  3  <  0.4.  The  integrated  intensity  for 
the  second  term  of  Eq.  ( 12)  is  obtained  from  Eq.  ( 10): 


1 

Qi  *  —  |  ,rJ  |.4;lp.i)|2<ip  *  2dmu4d./(7'W>:!  (141 

The  contribution  to  the  output  from  the  third  term  of 
Eq.  (12)  is  obtained  by  combining  Eqs.  (5)  and  ( 10): 


1  ym  » 

QsUi  »  —  2  Rel<4llp.fia:(p,f  )|<ip 

2r  Jpk--r~ 
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We  see  that  the  third  term  is  proportional  to  the 
magnitude  of  the  signal  multiplied  by  a  constant  tem¬ 
poral  offset  frequency  Pd<-’  *  u.'a  which  can  be  set  at  a 
convenient  value.  This  term  can,  therefore,  be  sepa¬ 
rated  from  the  signal  and  reference-beam  bias  terms, 
as  given  by  Eqs.  ( 13)  and  ( 14),  by  postdetection  filtering. 
Since  Q:iU>  is  not  a  function  of  p*,  the  same  temporal 
frequency  offset  is  provided  by  the  reference  beam  for 
any  input  signal.  The  reference  beam,  in  this  sense,  is 
a  distributed  local  oscillator  whose  temporal  frequency 
can  be  changed  by  a  simple  geometric  adjustment. 

The  total  optical  power  collected  by  a  photodetector 
is  the  sum  of  Eqs.  ( 13),  (14),  and  (15).  When  the  optical 
power  is  multiplied  by  the  photodetector  sensitivity  S. 
expressed  in  A/W,  and  we  account  for  the  optical  ef¬ 
ficiencies  of  the  signal  and  reference  beams  (<,  and  er), 
we  find  that  the  photodetector  output  current  is 


ittl  -  ti  +  it  +  i3<t> 

”  IdmiAU.LSCi  +  2d  miA'UrLS/(TW )2 

4dmlm2AiA-!  \  (,t.LS  „ 
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Each  photodetector  in  the  Fourier  plane  produces  a 
signal  similar  to  i  ( t )  except  that  the  value  of  C*  will  vary 
depending  on  the  strength  of  that  frequency  component 
in  the  input  signal.  Because  the  output  current  is 
proportional  to  C*  instead  of  |  C*  J 2.  as  it  would  be  for 
a  power  spectrum  analyzer,  the  improvement  in  dy¬ 
namic  range  is  significant.  We  now  consider  the  factors 
that  determine  the  value  of  the  parameter  3  associated 
with  tie  photodetector  geometry  and  postdetection 
bandwidth. 


ill.  Photodetector  Geometry  and  Bandwidth 

In  Fig.  2  we  show  the  light  distribution  in  the  Fourier 
plane  caused  by  two  cw  signals  that  are  just  resolved  by 
the  Rayleigh  criterion.  For  sake  of  clarity,  only  the 
central  lobe  of  the  sine  functions  as  given  by  Eq.  (5)  are 
shown:  we  can  also  represent  the  signals  by  sinc(/T) 
which  show  that  the  signals  are  separated  by  I  T  Hz. 
The  reference  beam  is  shown  as  a  uniform  amplitude 
tight  distribution  whose  temporal  frequency  is  offset  by 
a  fixed  amount  from  that  of  the  central  frequency  of  the 
sinci  f'T)  functions.  That  is.  if  /  is  the  temporal  fre¬ 
quency  of  a  cw  -signal.  the  reference-beam  frequency  at 
the  centroid  is  ;r  *  f  +  where  fa  *  p,3c/‘2~. 


If  we  have  R  photodetectors  per  resolvable  frequency, 
the  center  spacing  is  l/RT.  and,  if  the  duty  cycle  is  0  < 
d  <  1.  the  width  of  each  photodetector  is  d/RT.  Sup¬ 
pose  that  the  postdetection  bandwidth  is  a  rectangular 
function  of  width  b/T  centered  at  fa-  The  effect  of  the 
postdetection  bandwidth  is  to  accept  light  if  the  dif¬ 
ference  frequency  between  the  cw  signal  and  the  refer¬ 
ence  beam  falls  within  the  bandpass  and  to  reject  it 
otherwise.  As  an  example,  suppose  that  the  cw  signal 
has  a  frequency  of  300  MHz  and  that  the  reference- 
beam  frequency  is  adjusted  geometrically  so  that  it  is 
310  MHz.  The  interference  beat  frequency  will  then 
be  fd  m  10  MHz,  and  the  narrowband  postdetection 
filter  will  accept  all  the  light  near  the  centroid  of  the  sine 
function  associated  with  the  signal.  As  noted  in  Sec. 
II.  the  entire  sine  function  oscillates  at  a  temporal  fre¬ 
quency  of  300  MHz  so  that,  as  we  move  the  photode¬ 
tector  away  from  the  300- MHz  position,  we  continue  to 
collect  light  until  the  difference  between  the  reference 
and  signal  frequency  exceeds  U  ±  b/2T.  For  example, 
if  the  frequency  resolution  of  the  system  is  1/T  *  3  MHz 
and  if  b  is  chosen  to  be  1,  we  find  that  when  the  photo¬ 
detector  is  placed  at  a  position  corresponding  to  301.5 
MHz  (or  298.5  MHz),  the  output  of  the  narrowband 
filter  rapidly  falls  to  a  low  value. 

The  narrowband  filter  is  therefore  highly  effective  in 
reducing  cross  talk  in  the  spectrum  analyzer.  The 
sidelobe  levels  decay  rather  slowly  so  that,  in  the  ab¬ 
sence  of  filtering,  they  contribute  to  an  erroneous 
measurement  of  the  true  signal  levels  at  nearby 
frequencies.  Since  the  filter  can  suppress  out-of-band 
signals  by  30-40  dB  in  addition  to  the  inherent  roll-off 
of  the  sidelobes,  it  is  not  necessary  to  control  the  side- 
lobe  levels  by  input  aperture  weighting. 

The  postdetection  filtering  does  not,  however,  pre¬ 
vent  the  detection  of  narrow  pulses.  In  Sec.  II  we 
showed  that  the  spectrum  of  a  pulse  whose  duration  is 
«T  is  a  sine  function  whose  sidelobes  oscillate  at  a 
frequency  proportional  to  their  distance  from  the  cen¬ 
troid.  The  interference  frequency  is  therefore  fixed 
over  a  large  spatial  range  so  that  any  photodetector 
within  the  sine  function  envelope  will  respond  with  the 
proper  output.  This  result  can  be  seen  by  noting  that 
Eq.  (8)  contains  a  term  exp i-jpvt)  that,  when  multi¬ 
plied  by  the  term  exp[+;(p  +  Pd)ut]  from  the  reference 
beam  as  given  by  Eq.  (10),  produces  a  fixed  interference 
frequency. 

A  second  benefit  of  narrowband  filtering  is  in  re¬ 
solving  two  signals  closely  spaced  in  frequency.  As 
shown  in  Fig.  2,  two  received  signals  spaced  by  1/T  Hz 

—  T  *H*f*Ct  t*Ml 


Fig.  2.  Relerence  and  sienal  beam  nistribution  in  the  Fourier  plane 
tor  two  cw  signals. 
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may  occur  at  any  position  relative  to  a  set  of  photode- 
tector  elements.  We  must  ensure  that  the  response 
from  a  photodetector  is  sufficiently  high  when  it  is 
sampling  the  light  at  the  centroid  of  any  sine  function 
to  provide  accuracy  and  sufficiently  low  at  the  midpoint 
between  sine  functions  to  provide  resolution.  Since  the 
received  signals  may  occur  at  any  position  relative  to  the 
photodetectors,  we  must  find  a  way  to  relate  the  number 
of  photodetectors  per  frequency,  the  spatial  duty  cycle, 
and  the  poetdetection  bandwidth  to  the  required  dip 
between  resolvable  frequencies. 

The  first  step  is  to  determine  the  response  of  a  fixed 
photodetector  toacw  signal  as  a  function  of  its  fre¬ 
quency.  The  output  of  the  photodetector  will  be  the 
convolution  of  a  sinc(/T)  function  with  the  rectangular 
function  of  width  d/RT  (representing  the  photodetec¬ 
tor),  followed  by  a  multiplication  with  the  bandpass 
fn  "er  function.  The  postdetection  filter  function  is 
apprv  vimated  by  a  uniform  response  for  |/j  <  b/2T  and 
a  lineariy  decreasing  response  from  |/f  ”  b/2T  to  j/j  ■ 
q/2T  as  shown  in  Fig.  3.  This  function  is,  aside  from 
pass  band  and  stop-bond  ripple,  similar  to  that  produced 
by  an  elliptic  filter.  Since  the  convolution  is  nearly 
constant  over  the  defined  range  of  frequencies,  the  net 
result  can  be  represented  by 

Ci/)  «  d/RT;  0  S|/|  S  b/2T, 

mi?r~k<~~b]i,~bnT)''  bnTs  1/1  iq/2T-cl  *b- 

-0 ;  /  £  <j/2T.  (1?) 

The  smallest  dip  between  two  frequencies  occurs  when 
the  midpoint  between  the  frequencies  falls  between  two 
photodetectors.  This  condition  is  shown  in  Fig.  4  where 
we  show  C(/)  centered  at  the  positions  corresponding 
to  two  resolvable  frequencies.  Since  we  have  already 
accounted  for  the  convolutional  effects  introduced  by 
photodetectors  having  a  finite  size,  we  can  now  deter¬ 
mine  the  response  from  a  given  photodetector  by  using 
the  value  of  C(/)  at  a  point  corresponding  to  the  center 
of  the  photodetector.  The  output  of  the  photodetector 
nearest  the  center  of  C (/)  is  equal  to  C (/o).  The  largest 
possible  value  of  |/o|  is  1/2RT,  and  we  require  that  the 
output  be  constant  for  any  frequency  less  than  |/o| .  By 
referring  to  Fig.  3,  we  see  that  this  condition  is  satisfied 
if 

Rsl/b.  (18) 

The  outputs  from  the  two  photodetectors  symmet¬ 
rically  positioned  about  the  midpoint  are  equal  and  have 
values  given  by  C (f\).  The  second  relationship  can  then 
be  obtained  by  requiring  that 

Cifo  *  hClfa).  (19) 

where  0  <  h  <  1  is  the  desired  dip  between  the  two 
frequencies.  Since  f\  »  1/2T  -  1/2 RT.  we  use  Eqs.  ( 17) 
and  <  19)  to  obtain  the  relationship  that 

R  > - - - •  >20> 

1  -  </i  1  -  h  i  -  hb 

The  inequalities  as  given  by  Eqs. '  IS)  and  *  20)  must 
be  satisfied  to  achieve  a  given  dip  between  frequencies 


Fig.  3.  Convolution  of  a  sinc(f  7b  function  and  a  photodetcctor  of 
width  d/RT. 


Fig.  4.  Worst -ca»*  signal  position  relative  to  the  photodetectors  for 
resolving  two  frequencies. 


as  well  as  to  ensure  that  the  magnitudes  of  the 
frequencies  are  accurately  measured.  In  the  limit  as  the 
slope  of  the  filter  becomes  very  high,  we  have  that  q  * 
b  and  Eq.  (20)  reduces  toR  >  1/(1  -  6).  If  Eq.  (18)  is 
also  satisfied,  we  find  that  R  *  2  is  the  minimum  num¬ 
ber  of  photodetectors  per  resolvable  frequency. 

For  h  »  0.5.  which  corresponds  to  a  3-dB  dip  between 
frequencies.  Eq.  (20)  becomes 


Again,  if  we  use  the  equality  from  Eq.  (18),  and  let  q  * 
2b.  we  have  that  b  »  0.4  and  R  »  2.5. 

One  of  the  key  advantages  of  using  narrowband 
postdetection  filtering,  then,  is  to  provide  the  necessary 
dip  between  cw  signals  with  a  number  of  photodetectors 
per  frequency  that  is  close  to  the  theoretical  limit  The 
natural  dip  between  the  two  Rayleigh  resolved 
frequencies  will  vary  from  zero  to  1.3  depending  on  the 
relative  phases  of  the  two  spectral  components. 
Equation  (20)  shows  that  a  fixed  dip  of  value  h  can  be 
achieved  independently  of  the  phase  relationship 
without  increasing  the  length  of  the  Bragg  cell  to  sepa¬ 
rate  physically  the  spectral  components. 

We  now  calculate  the  relationship  between  R  and  the 
Bragg  cell  length  for  a  power  spectrum  analyzer. 
Gaussian  illumination  is  generally  used  to  control  the 
sidelobe  levels  to  meet  the  cross  talk  requirements. 
This,  in  turn,  causes  a  loss  in  resolution  which  must  be 
restored  by  increasing  the  Bragg  cell  length.  We  denote 
the  increased  cell  length  by  L  i  30  that  the  input  aperture 
illumination  function,  in  terms  of  amplitude,  is  given 
by 
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am 


-  .4,exp<-4rf  xVL{).  (22) 

where  T\  is  a  truncation  ratio  whose  value  is  chosen  to 
control  the  sidelobe  level.7  The  intensity  response  of 
the  system  to  a  cw  input  frequency  is  then  obtained  by 
squaring  the  Fourier  transform  of  a{x ): 

/(pi  a  exp( ~p2L'l/9T\).  (23) 

By  following  the  same  procedure  as  before,  we  find  that 
the  convolution  of  I(p)  with  a  small  photodetector  of 
width  d/RT  yields  a  normalized  result  that  is  very  close 
to  being  Gaussian: 


where  we  have  replaced  pv  with  2  iff.  The  worst-case 
condition  for  resolving  two  frequencies  occurs  when  the 
midpoint  between  those  frequencies  lies  at  the  midpoint 
between  two  photodetectors.  Since  the  intensities  add, 
the  relationship  that  must  be  satisfied  is 

c<fi)  +  c(/i+^;)  <«?(/<>).  125) 

where  /i  *  (R  - 1  )/2 RT  and  the  maximum  value  of  |/o| 
is  1/2RT  as  before.  For  a  dip  of  h  =  0.5  between 
frequencies,  we  find  that  Eq.  (25)  is  satisfied  when 

Li  -O-STiLR/iR  -  1).  (26) 

Suppose  we  set  R  =  2.5  for  comparison  purposes  and  set 
T i  «  1.7  which  ensures  that  the  maximum  sidelobe  is 
down  by  at  least  40  dB.7  We  then  find  that  L  i  *  2.15L 
which  shows  that  the  length  of  a  Bragg  cell  used  in  a 
power  spectrum  analyzer  must  be  more  than  twice  as 
long  as  that  required  for  an  interferometric  spectrum 
analyzer.  The  impact  of  the  larger  aperture  required 
in  a  power  spectrum  analyzer  is  significant  since  the 
increased  length  translates  directly  into  the  need  for  an 
increased  time-bandwidth  product  to  achieve  equiva¬ 
lent  performance. 

Another  advantage  of  the  interferometric  approach 
is  that  any  scattered  noise  caused  by  stationary  ele¬ 
ments  in  the  system  is  rejected  by  the  narrowband  filter 
because  it  does  not  interfere  to  produce  in-band  com¬ 
ponents.  This  feature  also  causes  the  light  distribution 
at  the  optical  axis  to  be  rejected. 

We  conclude  this  section  by  calculating  the  value  of 
the  parameter  3  which  is  equal  to  one  hall''  of  the  ratio 
of  the  width  of  the  photodetector  d/RT  to  the  spacing 
1/T  between  resolvable  frequencies.  We  choose  d  =  0.7, 
which,  for  R  *  2.5.  gives  a  value  for  3  equal  to  0. 14;  this 
value  justifies  several  assumptions  made  in  Sec.  II  in 
connection  with  deriving  Eq.  ( 16). 

IV.  Reference  Waveform* 

In  Sec.  II  we  used  a  narrow  pulse  as  the  reference- 
beam  waveform  to  develop  the  basic  theory  of  a  fixed 
offset  frequency  interferometric  spectrum  analyzer. 
The  narrow  pulse  approximates  a  delta  function  that 
provides  most  of  the  desired  properties  of  the  reference 
waveform  in  the  Fourier  plane  which  are  that  1 1  >  the 
amplitude  should  be  uniform  in  spatial  trequencv;  t'2) 
the  spatiai  and  temporal  frequencies  should  be  coupied 


so  that,  with  a  geometric  displacement  of  the  reference 
waveform,  a  fixed  offset  frequency  is  produced  at  each 
photodetector  location;  (3)  the  amplitude  should  not  be 
a  function  of  time;  (4)  the  duty  cycle  should  be  high  so 
that  short  duration  signals  are  not  missed;  and  (5)  the 
light  is  efficiently  used.  The  major  drawback  to  the  use 
of  a  narrow  pulse  as  the  reference  signal  is  that  a  fraction 
of  only  l/TW  of  the  light  at  the  Bragg  cell  is  intercepted 
by  the  pulse. 

An  equivalent  way  to  generate  a  traveling  pulse  is  to 
use  a  chirp  waveform  to  activate  a  Bragg  cell  situated 
so  that  the  focused  beam  scans  across  plane  P3.  If  we 
project  the  rays  produced  by  a  traveling  impulse  from 
plane  P3  toward  the  source,  we  find  that  the  Bragg  cell 
must  have  a  length  equal  to  2 L,  given  that  the  chirp  rate 
is  equal  to  the  bandwidth  of  the  signal  divided  by  the 
processing  time  T.  The  focal  length  of  the  chirp  is  then 
equal  to  L  -/\TW,  which  is  of  the  order  of  several  meters 
for  typical  values  of  these  parameters.  The  advantages 
of  this  approach  are  that  sharply  formed  impiilses  are 
produced  at  plane  P. j  and  that  the  resultant  reference 
beam  in  the  Fourier  plane  is  a  plane  wave.  The  disad¬ 
vantages  are  that  a  Bragg  cell  of  twice  the  time-band¬ 
width  product  is  required  and  the  length  of  the  optical 
system  is  excessive. 

Another  way  to  generate  the  reference  beam  is  to 
drive  the  Bragg  cell  in  plane  P.i  with  a  chirp  directly.8 
Lens  L-i  then  focuses  the  chirp  at  plane  P4  so  that,  at  the 
Fourier  plane,  the  expanding  beam  covers  the  signal 
spectrum.  The  chirp  waveform  can  be  represented 
by 

I  rW  \ 
r(t)  *  cos|wf£  +  —  f2l 

for  a  chirp  of  increasing  frequency  on  a  carrier  ivc.  The 
Fourier  transform  of  the  chirp,  valid  over  the  spatial 
frequency  range  corresponding  to  the  bandwidth  of  the 
signal,  is 

I  [  T  W 

^  ^  expj-;  ujc  +  -^7-  ( C  -  T  -  r/v  I 

X  (t  -  P  -  x/ioj  exp(-/lp  +  pt)x\dx 

jmnA^L  .  . 

-  — ==-  expijtfl  exp  -;(p  +  pc  )vt | 

v  TAT 

(  uL 

x  exp  ->  —  (p  +■  Pd  -  Pc'2  ,  127) 

|  )»lv 

where  2 irW  is  the  bandwidth  of  the  rf  signal  and  the 
chirp  rate  is  2xH7T.  We  see  that  Eq.  (27)  is  similar  to 
Eq.  1 10)  but  with  two  notable  differences:  (1)  the  am¬ 
plitude  is  larger  by  v  TW.  which  increases  the  efficiency 
of  the  system:  and  (2)  the  sine  function  envelope  is  re¬ 
placed  bv  a  Fresnel  amplitude  diffraction  pattern. 

Figure  5(a)  shows  the  focused  spot  at  plane  P»  for  the 
instant  in  time  when  the  Bragg  cell  is  just  filled  with  one 
chirp  function.  At  plane  Pj,  the  spatiai  extent  of  the 
reference  beam  is  equal  to  that  of  the  signal  spectrum: 
the  amplitude  is  uniform  except  near  the  edges  of  the 
spectrum  where  the  ripple  effects  of  the  Fresnel  dif¬ 
fraction  pattern  are  evident.  At  a  time  T  2  sec  later, 
the  reference  waveform  consists  of  parts  of  two  chirp 
functions  that  produce  two  focused  spots  in  plane  P4  as 


15  August  1981  V01.  20.  No  '6  AP°LlED  OPTICS 


27-5 


Fig.  5.  Chirp  reference  beam:  fa)  Bragg  ceil  with  one  chirp  segment, 
and  ib>  Bragg  cell  with  one  half  of  two  chirp  segments. 


shown  in  Fig.  5(b).  Each  spot  has  increased  in  width 
by  a  factor  of  2  but  each  expands  to  cover  only  one  half 
of  the  spectral  range.  Thus,  as  the  chirp  waveforms 
move  through  the  Bragg  cell,  the  reference  waveform 
in  the  Fourier  plane  will  generally  consist  of  two  seg¬ 
ments  of  Fresnel  patterns  which  move  at  the  same  ve¬ 
locity.  If  the  chirp  duty  cycle  is  not  100%,  there  will  be 
a  small  gap  between  segments  which  should  not  ad¬ 
versely  affect  the  performance  of  the  system. 

When  the  chirp  waveform  is  used,  the  integrated  in¬ 
tensity  of  the  desired  output  term  is  similar  to  Eq.  (15), 
except  that  TW  is  replaced  by  v  TW,  and  the  chirp 
produces  a  residual  spatial  fluctuation  that  is  a  potential 
source  of  time-amplitude  modulation  on  the  output 
current  of  the  photodetector.  Since  the  Fresnel  pattern 
travels  at  a  high  velocity,  the  modulation  is  well  outside 
the  postdetection  bandwidth  so  that  there  is  no  tem¬ 
poral  modulation  on  the  output  due  to  spatial  fluctua¬ 
tions. 

A  second  type  of  time-amplitude  modulation  of  the 
reference  function  at  the  Fourier  plane  arises  from  the 
movement  of  the  input  chirp  waveform  under  the 
Gaussian  input  illumination.  A  frequency  wii.hin 
the  bandwidth  of  the  chirp  will  appear  at  a  fixed  posi¬ 
tion  pj  in  the  Fourier  plane;  its  amplitude  as  a  function 
of  time,  however,  is  determined  by  its  position  in  the 
input  plane.  Since  a  severe  truncation  of  the  Gaussian 
illumination  is  not  needed  to  control  the  sidelobe  levels, 
we  can  reduce  the  amplitude  variation  at  the  expense 
of  some  loss  of  light.  A  reasonable  compromise  is  to 
truncate  the  Gaussian  illumination  at  the  exp<— ¥2) 
points  at  the  edges  of  the  aperture  which  leads  to  a 
±25%  variation  in  the  output  current  over  the  time  in¬ 
terval  T.  This  time-amplitude  modulation  is  not  a 
concern  because  it  can  be  eliminated  by  the  bandpass 
filter  if  fd  »  1  IT.  The  chirp  reference  waveform 
therefore  satisfies  most  of  the  requirements  set  forth  at 
the  beginning  of  this  section. 

An  alternative  reference-beam  waveform  that  is 
continuous  in  time  is  a  pseudonoise  sequence  iPNS) 
whose  Fourier  transform  also  satisfies  many  of  the 
stated  requirements.  These  sequences  and  similar 


codes  have  been  studied  for  applications  such  as  gen¬ 
erating  random  phase  masks  for  holography,  recovering 
imagery  from  coded  aperture  systems,  and  spectral 
shaping.9-16  Generally,  only  the  spatial  properties  of 
these  codes  have  been  considered;  we  are  interested  in 
both  the  spatial  and  temporal  properties  of  their  Fourier 
transforms  to  determine  whether  they  are  better  than 
those  of  a  chirp  waveform. 

A  shift  register  sequence  is  one  whose  period,  or  frame 
length,  is  N  <  2r  -  1,  where  .V  is  the  number  of  ele¬ 
ments  in  the  frame  and  r  is  an  integer.  If  the  period  is 
equal  to  2r~l,  the  sequence  is  said  to  have  maximal 
length  and  has  the  properties  (1)  that  in  every  period 
the  difference  between  the  number  of  +  l’s  and  —  l’s  is 
eque t  >  one;  (2)  that  in  every  period  the  longest  run  of 
+l’s  and  -l’s  is  equal  to  r,  and  for  each  run  of  length  m 
+  1  there  are  two  runs  of  length  m;  and  (31  that  the  au¬ 
tocorrelation  function  has  only  two  values.  Such  se¬ 
quences  can  be  produced  by  an  r -stage  linear  shift  reg¬ 
ister  having  the  appropriate  feedback.17 

Figure  6  shows  the  magnitude  of  the  Fourier  trans¬ 
form  of  an  .V  *  26  —  1  *  31  element  sequence  at  time  t 
*  0  and  t  *  T/31;  that  is,  in  a  time  interval  corre¬ 
sponding  to  one  shift  position  of  the  PNS  in  the  input 
plane,  the  magnitude  changes  by  the  amount  shown. 
At  the  integer  frequencies  the  normalized  magnitude 
is  exactly  equal  to  one  for  any  shift  position.17  The 
Fourier  transform  is  symmetric  about  the  spatial  fre¬ 
quency  iV/2;  in  an  optical  system  the  sequence  modu¬ 
lates  a  carrier  so  that  the  spectrum  is  centered  at  pr.  In 
practice  the  input  waveform  will  consist  of  a  sequence 
of  narrow  pulses  so  that  the  spectrum  must  be  shifted 
and  multiplied  by  a  sine  function  as  given  by  Eq. 
GO). 

As  shown  in  Sec.  III.  we  need  at  least  2.5  photode¬ 
tectors  per  resolvable  frequency  which  means  that  some 
photodetectors  will  be  located  where  the  spectrum  is  not 
uniform  in  either  space  or  time.  F urther.  since  photo- 
detectors  of  finite  size  must  be  used  to  collect  optical 


Fig;,  o.  Magnitude  of  the  Fourier  ’ranstorm  if  i  pseudonoise  se¬ 
quence  for  two  adjacent  shift  positions. 
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Fig.  7.  Magnitude  as  a  function  of  shift  position  for  four 
frequencies. 


power  efficiently,  we  investigated  the  spatial  and  tem¬ 
poral  variations  in  amplitude  for  finite  photodetectors. 
Figure  7  shows  the  magnitude  variation  as  a  function 
of  shift  position  (equivalent  to  time)  for  a  photodetector 
whose  width  is  d/RT  »  0.28,  centered  at  four  different 
spatial  frequencies.  The  variation  in  magnitude  as  a 
function  of  time  is  of  the  order  of  ±15%,  which  is  caused 
by  using  a  Gaussian  input  illumination  beam  truncated 
at  the  expi-1^)  points  in  amplitude  at  the  edges  of  the 
sequence.  We  note  that  the  functions  are  the  sums  of 
sinusoidal  functions  whose  frequencies  are  1  IT  and 
multiples  thereof;  these  components  can  also  be  re¬ 
moved  by  the  narrowband  filter.  As  noted  before,  the 
instantaneous  output  of  a  spectrum  analyzer  is  generally 
integrated  for  at  least  T  sec.  Figure  8  shows  the  average 
value  of  the  amplitude  as  a  function  of  spatial  frequency 
for  a  full  frame  of  the  PNS.  Except  at  zero  spatial 
frequency,  where  the  magnitude  is  low  because  the  av¬ 
erage  value  of  the  sequence  is  low,  the  average  magni¬ 
tude  varies  by  ±5%.  These  variations,  along  with  others 
such  as  Bragg  cell  frequency  roll-off  and  the  sine  func¬ 
tion  weighting,  can  be  compensated  when  the  system 
is  calibrated. 

The  question  of  whether  the  temporal  frequencies  of 
these  candidate  waveforms  are  continuous  functions  of 
the  spatial  frequency  variable  depends  on  the  waveform 
and  the  postdetection  process.  Equations  (8)  and  (27) 
suggest  that  the  spectrum  of  a  moving  impulse  or  a  chirp 
waveform  is  indeed  continuous  so  that  the  theory  is 
valid  if  the  output  is  sampled,  after  narrowband  filtering 
and  envelope  detection,  once  for  each  time  period  T. 
The  results  are  less  clear  for  the  PNS  because  the 
Fourier  transform  cannot  be  calculated  in  the  same  way. 
Some  insights  ran  be  gained  for  the  case  of  longer  ob¬ 
servation  times  such  that  the  reference  waveform  is 
repetitious  with  period  T. 

Let  fix)  represent  the  amplitude  of  the  reference 
beam  over  the  time  interval  T  and  let  Ftp >  be  its  Fourier 
transform.  Suppose  that  this  distribution  is  repeated 


at  intervals  of  L  and  that  the  entire  distribution  moves 
with  velocity  u.  The  Fourier  transform  of  the  extended 
signal  g(x)  is  then 

Gip)  ■  u  exp i—jpvt)  F(pv)  — ■  £  5(pi>/2x  -  n/T).  (28) 

T 

We  now  account  for  the  finite  aperture  in  the  input 
plane  by  convolving  Gip)  with  the  aperture  function 
A(p)  ■  L  sinc(pl/2ir)  to  get 

Gip)  *  £  F  |~^”|  ttpi-jlxnt/T)  *inc(pii772»  —  a >.  (29) 

Equation  (29)  shows  that  the  Fourier  transform  consists 
of  a  set  of  sine  functions  spaced  at  intervals  of  2 ir/L, 
which  is  consistent  with  our  notion  of  the  spatial  reso¬ 
lution  of  the  optical  system.  These  sine  functions  take 
on  the  sample  values  F(2tnx/T)  and  oscillate  at  a  tem¬ 
poral  frequency  of  2 irn/T.  The  final  step  is  to  multiply 
G(p)  by  sinc(pL/2irTW)  to  account  for  the  finite  du¬ 
ration  of  each  pulse  in  the  PNS  or  by  a  rectangular 
function  for  the  case  of  a  chirp. 

If  we  observe  the  output  of  a  given  photodetector  for 
a  long  time  period  T'  »  T,  we  will  see  the  discrete  na¬ 
ture  of  the  temporal  spectrum  of  the  reference  beam. 
Furthermore,  if  we  increase  the  frame  length  of  the  PNS 
so  that  it  is  greater  than  T  but  less  than  T',  we  will  note 
that  the  discrete  temporal  frequency  resolution  will 
increase  whereas  the  spatial  frequency  resolution  is 
fixed.  The  exact  nature  of  the  output  for  the  candidate 
waveforms  and  an  observation  time  of  the  order  of  T 
will  be  deferred  to  a  later  paper  where  we  will  discuss 
the  experimental  results. 

V.  Laser  Power  and  Dynamic  Range 

The  current  produced  by  a  photodetector  is  given  by 
Eq.  (16);  the  values  of  Ai  and  A2  are  given  by 


•At  - 

CIO) 

w 

e  ) 

L  1 

where  P  is  the  power  of  the  laser  and  a  is  the  fraction 
of  the  power  directed  into  the  signal  beam.  The  output 
current  can  then  be  arranged  in  the  form 


a  ?  a  s  ;a 


SP»TI*l  raEOUENC'' 

Fig.  ■).  Magnitude  averaged  over  a  frame  period  as  a  function  if 
spatial  frequency. 
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Fif.  9.  Dynamic  range  for  interferometric  and  power  spectrum 
analyzers. 


iit)  ■  i\  +  ij  +  13(0 


*  aP (1  —  a)  +  bPa  +  cP(a(l  —  a)]l/2  coeluat  +  <#),  (31) 

where  the  coefficients  a,  6,  and  c  contain  the  values  of 
parameters  associated  with  various  system  components. 
Equation  (31)  represents  the  general  form  of  the  output 
of  any  interferometric  system,  and  the  SNR  is  given  by 
the  ratio  of  the  output  signal  power  to  the  sum  of  the 
shot  noise  and  thermal  noise  powers: 


SNR  * 


_ <ii2<t))R _ 

2 cfl(ia  +  ii  +  +  4 kTB 


(32) 


where  the  brackets  denote  time  average,  R  is  the  load 
resistance  of  the  photodetector,  e  is  the  charge  of  an 
electron,  B  is  the  bandwidth,  ig  is  the  photodetector 
dark  current,  k  is  Boltzmann’s  constant,  and  T  is  the 
equivalent  system  temperature  in  degrees  kelvin.  By 
using  Eq.  (31)  in  (32),  we  obtain 


SNR  -  - 


0.oc2Pla(  1  —  a) 


(33) 


D[aP(l  —  <*)  +  6Pa|  +  F 
where  D  *  2 eB  and  F  m  (2 eBi*  +  4 kTB/R). 

We  want  to  find  the  minimum  laser  power  required 
to  achieve  a  given  dynamic  range.  The  procedure  is  to 
solve  Eq.  (33)  for  the  laser  power  required  to  produce 
a  detectable  signal  (a  SNR  of  1)  when  c  has  its  minimum 
value  as  determined  from  the  required  dynamic  range. 
After  some  algebraic  manipulations,  we  find  that  the 
minimum  laser  power  is  obtained  when 

g  4*  1 

a  ■  — -  (34) 

iff  +  1 

where  g  *  ( Fc  2/2D  -a  2) 1/2.  The  minimum  laser  power 
is  then 

„  20a  1 2g  1 ) 

P™  n-  - - -  '35) 


We  are  now  in  a  position  to  make  some  numerical 
calculations  using  Eq.  <  16)  as  our  basis  for  determining 
the  values  of  the  constants  in  Eq.  1 31).  We  will  assume 
that  the  efficiency  of  each  optical  path  is  1/2. 


that  the  sensitivity  of  the  photodetector  is  S  »  0.5  A/W, 
and  that  d  *  0.14.  The  modulation  indices  m  i  and  mo 
are  chosen  to  keep  the  intermodulation  products  at  an 
acceptable  level.  If  we  require  third-order  two-tone 
intermodulation  products  to  be  down  by  at  least  50  dB 
when  20%  of  the  possible  cw  signals  in  the  Bragg  cell  are 
saturated,  we  must  keep  the  diffraction  efficiency  per 
frequency  to  <0.01.  Since  the  modulation  index  is 
equal  to  the  square  root  of  the  diffraction  efficiency,  we 
have  that  m  1  ■  0.1.  The  diffraction  efficiency  of  the 
reference  beam  can  be  much  higher  because  inter¬ 
modulation  products  are  not  a  factor;  we  choose  m2  “ 
0.7.  We  will  also  use  a  time-bandwidth  product  of  TW 
=  100,  and,  if  we  use  the  chirp  reference-beam  wave¬ 
form,  we  replace  TW  by  \'TW  in  Eq.  (16)  so  that 

(It)  ■  «i  +  i*  +  13U)  *36) 

-  7(  lO~*)CiPa 

+  3.5(10“4)P(l  -  «>  +  9.8<10-4>C„Plcr(l  -  a)|1/2  •  (37) 

from  which  we  note  that  a  =  3.5(l0-4),<>  =*  7(10-4)C», 
and  c  *  9.8(10-4)C*. 

The  value  of  g  is  a  function  of  the  parameters  asso¬ 
ciated  with  the  photodetector  elements  and  the  required 
dynamic  range.  A  dynamic  range  of  60  dB  in  terms  of 
input  signal  power  yields  a  minimum  value  of  C*  =  10“3 
from  which  we  have  that  cmm  =  9.8(  10-T).  We  chose  a 
photodetector  having  a  dark  current  of  id  »  10~9  A,  a 
bandwidth  of  B  *  1  MHz,  and  a  load  resistance  value 
of  R  «  50  kO.  From  these  parameters  we  calculate  that 
g  *  7.2  and,  from  (35),  that  the  minimum  laser  power 
is  P  *  3.2  mW. 

We  now  calculate  the  dynamic  range  produced  by  a 
power  spectrum  analyzer  that  has  the  same  system 
parameters.  The  SNR  for  a  power  spectrum  analyzer 
can  be  given  in  a  form  similar  to  Eq.  (32): 

SNR  -  — - - ^ - -  .  ,38) 

2eBU*  4-  :,)R  +  4 kTB 

where  i  1  ■  7(  10-4)C*P;  we  have  set  a  *  1  because  all  the 
laser  power  can  be  directed  into  the  signal  beam.  The 
dynamic  range  can  be  found  by  setting  the  SNR  ■  1:  we 
then  have 

DRp  -  10  lo«(7(10-4)P/v ?!,  (39) 

which,  for  a  laser  power  of  3.2  mW,  gives  a  dynamic 
range  of  33  dB.  We  find  that  the  interferometric 
spectrum  analyzer  provides  27  dB  more  dynamic  range 
than  a  power  spectrum  analyzer  for  the  set  of  parame¬ 
ters  given. 

It  is  frequently  cited  in  the  literature  that  an  inter¬ 
ferometric  spectrum  analyzer  provides  twice  the  dy¬ 
namic  range  in  decibels  as  that  produced  by  a  power 
spectrum  analyzer.  We  now  examine  in  more  detail  the 
relationships  between  the  interferometric  and  power 
spectrum  analyzer  dynamic  ranges.  In  Fig.  9  we  graph 
the  output  signal  power,  as  given  by  the  numerators  of 
Eq9.  (32)  and  (38),  as  a  function  of  the  rf  input  signal 
power.  The  noise  floor  is  established  by  the  noise 
equivalent  power  of  the  photodetector,  except  at  high 
input  signal  levels  where  the  system  has  a  slightly  higher 
noise  floor  caused  by  signal  dependent  shot  noise. 
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As  expected,  the  output  vs  input  line  for  a  power 
spectrum  has  a  slope  of  two,  whereas  the  slope  is  one  for 
an  interferometric  spectrum  analyzer.  For  a  given  laser 
power,  the  vertical  axis  intercept  point  is  lower  for  the 
interferometric  spectrum  analyzer  because  <  1 )  half  of 
the  laser  power  is  lost  by  the  recombining  beam  splitter, 
and  (2)  the  reference-beam  power  is  spread  over  a  larger 
range  of  spatial  frequencies  so  that  the  reference-beam 
amplitude  at  any  particular  photodetector  is  low. 

The  dynamic  range  is  determined  by  the  points  at 
which  the  two  input/output  lines  intercept  the  noise 
floor.  A  relationship  between  the  interferometric  dy¬ 
namic  range  (DR/)  and  the  power  dynamic  range  (DR/»> 
can  be  obtained  by  substituting  the  laser  power,  as  given 
by  Eq.  ( 35 » .  required  to  provide  a  given  DR/  into  Eq. 
(39).  We  can  simplify  Eq.  (35)  for  the  case  of  g  very 
much  greater  than  one  to  the  form 

P»v8F7p  <40) 

Since  DR/  =>  20  log(C*)  and  c  ■  9.8(10~4)C*,  we  can 
rearrange  Eq.  (39)  in  the  form  of 

DR,  -  IDRp  -  10  logtTW/25).  (41) 

This  result  shows  that,  in  the  limit  of  small  time- 
bandwidth  products,  the  interferometric  system  does 
provide  twice  the  dynamic  range  in  decibels  as  expected. 
A  small  time-bandwidth  product  has  the  effect  of 
closing  the  gap  between  the  vertical  axis  intercept  points 
shown  in  Fig.  9  so  that  the  superior  performance  of  the 
interferometric  approach  is  evident  at  higher  input 
signal  levels.  We  also  note  from  Eq.  (41)  that  the  im¬ 
provement  factor  asymptotically  approaches  two  as  the 
absolute  performance  level  increases  for  a  fixed  value 
of  TW. 

VI.  Summary  and  Conclusions 

The  dynamic  range  of  a  power  spectrum  analyzer  is 
generally  limited  by  the  dynamic  range  available  from 
self-scanned  photodetector  arrays.  Furthermore,  the 
output  current  is  proportional  to  the  input  rf  power  so 
that  a  significant  amount  of  laser  power  is  needed  to 
achieve  a  high  dynamic  range.  Interferometric  spec¬ 
trum  analyzers  produce  an  output  that  is  proportional 
to  the  input  rf  magnitude  and  generally  provide  more 
dynamic  range  for  a  riven  laser  power. 

In  this  paper  a  technique  is  described  wherein  the 
reference  beam  is  spa  aally  modulated  and  geometrically 
shifted  so  that  the  interference  temporal  frequency  is 
constant  over  the  entire  spectrum.  This  fixed  fre¬ 
quency  offset  permits  a  narrowband  postdetection  filter 
to  separate  the  signal  term  from  the  bias  terms  at  the 
output.  Furthermore,  the  filter  rejects  sidelobe  con¬ 
tributions  from  nearby  cw  signals  so  that  reduced  cross 
talk  levels  can  be  achieved  without  the  need  to  weight 
the  input  illumination.  A  more  uniform  input  illumi¬ 
nation.  in  turn,  leads  to  better  short-pulse  detectability 
because  rhe  pulse  amplitude  is  more  constant  over  the 
processing  time  of  the  cell.  A  further  advantage  of  the 
interferometric  technique  is  that  scattered  light  does 
not  contribute  ro  the  output  because  it  fails  outside  the 
nandwidtn  of  rhe  filter. 


Two  types  of  reference- beam  modulation  have  been 
studied.  A  chirp  waveform,  having  the  same  bandwidth 
at  the  rf  signal,  provides  a  continuum  of  spatial  and 
temporal  frequencies  at  the  output  plane  and,  for  most 
applications,  has  adequate  amplitude  stability  as  a 
function  of  space  and  time.  We  also  investigated  the 
properties  of  pseudonoise  sequences;  the  same  post¬ 
detection  filtering  techniques  apply,  and,  for  very  short- 
pulse  or  cw  signals,  the  results  are  similar  to  those  ob¬ 
tained  from  the  chirp  waveform  except  that  the  sidelobe 
suppression  effect  may  not  be  as  strong  as  with  the  chirp 
waveform.  A  potential  advantage  of  the  pseudonoise 
sequence  waveform  is  that  it  can  be  recirculated  through 
the  Bragg  cell  and  provides  phase  continuity  at  the 
output  of  the  system.  The  choice  of  which  reference- 
beam  modulation  to  use  is  dependent  on  the  applica¬ 
tion. 

We  have  also  determined  the  minimum  laser  power 
required  to  achieve  a  given  dynamic  range  and  com¬ 
pared  the  performance  of  an  interferometric  and  a 
power  spectrum  analyzer  having  the  same  operating 
parameters.  The  dynamic  range  >f  the  interferometric 
system  is,  in  the  limit  of  small  time-bandwidth  prod¬ 
ucts,  a  factor  of  2  larger  in  decibels  than  that  of  a  con¬ 
ventional  system.  The  improvement  in  the  dynamic 
range  is  a  function  of  the  time-bandwidth  product  and 
the  absolute  dynamic  range  obtained. 
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Bragg  cell  diffraction  patterns 


A.  VanderLugt 


The  2-D  diffraction  patterns  produced  by  fully  illuminated  Bragg  cells  can  be  characterized  by  curved  singu¬ 
larity  functions.  The  degree  of  curvature  is  related  to  the  optical  wavelength,  the  acoustic  wavelength,  and 
the  degree  of  anisotropy  of  the  interaction  material.  Analytical  expressions  are  derived  for  the  singularity 
functions,  and  the  impact  of  the  nonideal  diffraction  pattern  in  certain  signal  processing  problems  is  calcu¬ 
lated.  The  most  notable  effect  is  an  increase  in  the  sidelobe  levels  and  a  slight  loss  in  resolution  when  Bragg 
cells  are  used  in  spectrum  analyzers.  Truncation  of  the  diffraction  pattern  or  the  use  of  the  line  illumina¬ 
tion,  when  possible,  reduces  the  degradation  of  the  diffraction  pattern. 


I.  Introduction 

Bragg  cells  are  widely  used  as  modulators,  beam  de¬ 
flectors,  and  analog  or  digital  delay  lines  for  signal 
processing.1  In  each  case  light  waves  interact  with 
sound  waves,  produced  by  an  electrical/acoustic 
transducer,  so  that  the  light  is  modulated  in  space  and 
time.  The  quality  of  the  optical  wave  front  is  particu¬ 
larly  important  in  applications  such  as  spectrum  anal¬ 
ysis,  wherein  the  resolution  of  the  system  is  significantly 
reduced  if  the  wave  front  is  distorted. 

One  method  for  testing  the  optical  quality  of  a  wave 
front  is  to  use  Fourier  analysis  techniques  which  give 
a  direct  measure  of  performance.  We  examined  the 
Fourier  transforms  of  various  Bragg  cells  operating  in 
the  longitudinal  wave  mode.  The  diffraction  pattern 
was  the  expected  2-D  sine2  function,  except  that  in  one 
dimension  the  function  was  curved  toward  the  optical 
axis. 

The  purpose  of  this  paper  is  to  analyze  this  phe¬ 
nomenon  for  both  isotropic  and  anisotropic  interaction 
materials,  to  illustrate  the  phenomenon  for  certain  types 
of  Bragg  cells,  and  to  calculate  its  effect  on  the  perfor¬ 
mance  of  some  optical  processing  systems. 

II.  Theory 

One  way  to  visualize  the  optical  wave  front  produced 
by  a  Bragg  cell  is  through  use  of  a  schlieren  imaging 
technique.  Figure  1  shows  a  schlieren  image  of  a  lead 
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molybdate  (PbMo04)  device  operated  in  the  longitu¬ 
dinal  mode  at  a  400  MHz  frequency.  In  the  near  field 
of  the  transducer,  the  acoustic  waves  interfere  to  pro¬ 
duce  an  intricate  pattern  that  has  been  studied  by  Cook 
et  al.2  A  null  is  formed  at  a  distance  ~H2/ 8 A  from  the 
transducer,  where  H  is  the  height  of  the  transducer  and 
A  is  the  acoustic  wavelength.  As  the  distance  from  the 
transducer  increases,  the  acoustic  waves  diverge  to 
produce  the  far-field  interference  pattern.  Several 
studies  have  been  made  in  which  the  sound  pressure 
fields  have  been  calculated  and  measured  by  optical 
probing  techniques.3-9  These  analyses  begin  with  the 
observation  that  the  transducer  is  equivalent  to  an  op¬ 
tical  source  so  that  the  sound  field  can  be  calculated  by 
using  the  Huygens-Fresnel  diffraction  formula.  The 
excitation  produced  by  a  finite  source  is  decomposed 
into  a  family  of  plane  waves,  and  the  far  field  pattern 
is  generally  expressed  in  terms  of  an  angular  spec¬ 
trum. 

We  take  a  related  approach  in  which  we  start  with  the 
sound  field  produced  by  a  transducer  of  infinitesimal 
height  and  account  for  its  finite  height  after  we  find  the 
Fourier  transform  of  the  coupled  optical  waveform. 
Our  model  for  the  acoustic  wave  front  is  described  with 
reference  to  Fig.  2.  The  Bragg  cell  is  located  in  plane 
Pi,  having  coordinates  x  and  y.  which  lies  in  the  front 
focal  plane  of  a  lens  having  focal  length  F.  The  Bragg 
cell  has  length  L  in  the  x  direction:  the  transducer  has 
a  height  H  in  the  y  direction  and  a  width  W  in  the  di¬ 
rection  of  the  optical  axis.  The  cell  is  illuminated  by 
a  uniform  collimated  beam  of  monochromatic  light. 
The  acoustic  wave  produced  by  a  wide  transducer 
having  a  very  small  height  can  be  represented  as  cylin¬ 
drical  sheets  that  diverge  in  the  y  direction  as  they 
propagate  in  the  x  direction.  We  shall  assume  that  the 
width  of  the  transducer  is  nearly  equal  to  the  width  of 
the  Bragg  cell  and  ignore  effects  due  to  reflections. 
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Fig.  1.  Schlieren  image  of  a  lead  molybdate  Bragg  cell  (400  MHz). 


Fig.  2.  Fourier  transform  setup  for  displaying  diffraction  patterns. 


The  acoustic  wave  front  propagation  within  the  cell 
is  similar  to  the  optical  wave  front  propagation  from  a 
line  source.  Given  an  ideal  interaction  between  light 
and  sound,  the  optical  wave  front  at  the  exit  side  of  the 
Bragg  cell  for  the  negative  first-order  diffracted  light 
is 


where  a  =  arctan(£/f’)  and  (5  =  arctan (tj/F);  the  vari¬ 
ables  £  and  rj  are  the  physical  coordinates  in  the  Fourier 
transform  plane.  We  now  change  to  polar  coordinates 
by  letting  x  =  r  cos 6  and  y  =  r  sin0  so  that  Eq.  (2)  be¬ 
comes 


J»  L  f  #72  \ 

(  .=  ezp(-;2rrr/A)  exp(-<ri 

0  J  —t/2  yj  j  Ar 

f  2irr 

X  exp  — —  (a  CO90  +  (3  sinl?)  rdrff,  (3) 

l  A 

where  the  region  of  integration  extends  over  the  semi¬ 
circle  in  the  right-hand  plane  of  the  Bragg  cell.  It  is  also 
convenient  to  convert  the  coordinates  in  the  Fourier 
plane  to  the  polar  form;  that  is,  let  a  -  a  cosy  and  0  = 
<r  siny.  Considering  the  integration  on  6  first,  we 
have 


a(xo')  “  i  :■■■■  exp(— y2xr/A)  exp(-«r),  (l) 

V]\r 

where  r2  *  x2  +  y2  and  t  is  a  parameter  that  accounts 
for  the  frequency  dependent  attenuation  of  the  acoustic 
wave  as  it  propagates.  We  have  suppressed  the  time 
dependence  of  this  propagating  wave  because  we  are 
interested  in  only  the  spatial  characteristics  of  the  re¬ 
sulting  diffraction  pattern.  The  optical  wave  front  as 
given  by  Eq.  (1)  does  not  account  for  the  Bragg  angle 
matching  conditions  which  change  the  amplitude  as  a 
function  of  the  acoustic  frequency.  These  effects  are 
not  of  fundamental  importance  here;  we  are  interested 
in  the  form  of  the  normalized  diffraction  pattern  of  a 
single  acoustic  frequency. 

We  first  consider  the  situation  for  a  Bragg  cell  made 
from  an  isotropic  medium  in  which  the  acoustic  veloc¬ 
ities  are  equal  in  all  directions.  The  Fourier  transform 
of  a(x,y)  is  given  by 


(i(r,a,y)  -  J  ^exp  J-;  —^—cos(7  -  0)|  d9,  4) 

which  is  independent  of  y  and  equal  to10 

G(r,a)  “  2*-«/o(2jrr<7/M,  (5) 

where  Jc(z)  is  the  zeroth-order  Bessel  function  of  the 
argument  z.  We  now  use  G(r,a)  in  Eq.  (3)  to  obtain 

A(<t)  »  Jo(2irr<r/X) 

Jo  \Jj\r 

X  exp(— ;2xr/A)  exp(— fr)dr,  (6) 

We  can  extend  the  upper  limit  to  infinity  if  we  multiply 
the  integrand  by  rect (r/L  — %),  and,  after  Eq.  (6)  has 
been  solved,  we  account  for  the  finite  length  of  the 
Bragg  cell  through  a  convolution  in  a  with  a  function 
W(<t),  which  is  the  Fourier  transform  of  the  product 
rect  (r/L  —  M2)  exp(— er).  The  solution  to  Eq.  (6)  »s  given 
by  the  relationship11 
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Fig.  3.  Phase  vector  K(4>)  and  Poynting  vector  Kty)  plotted  for  s 
=  -  ll2-  Both  vectors  follow  the  dashed  circle  for  isotropic  materials 
(s-0). 


where  F(-)  is  the  hypergeometric  function.  For  our 
application,  we  have  0=0  and  p  =  3/2  so  that 


.4(0)  ■ 


2*T(3/2) 


v7A^r(l)(a2  +  f>2)i'2 

where  a  =  j 2ir/A  and 


F  -1/4.3/4;!: 


if 

a2. 


(8) 


b  =  2ir<r/\. 


(9) 


Since  A(<r)  is  the  Fourier  transform  of  a  (nearly)  co- 
phasal  optical  wave  front  over  the  infinite  half-plane, 
we  expect  A  (<r)  to  be  a  singularity  function  that  behaves 
as  a  5  function.12  A  singularity  function  in  a  results  if 
b2  =  -a2,  provided  that  a  ^  0  and  that  F(- 1/4, 3/4;l;l) 
is  finite.  If  b2  =  —a2,  we  use  Eq.  (9)  to  obtain 


a2  =  a2  +  |32  =  (A/A)*.  UO) 

which  is  a  valid  solution  because  a  0,  and  we  have 
that11’ 


F(-l/4,3/4;l;l) 


0.54, 

r(5/4ira/4) 


so  that  the  hypergeometric  function  is  finite.  The 
equation  of  the  singularity  function  as  given  by  Eq.  (10) 
is  that  of  a  circle  having  a  radius  A/A  centered  on  the 
optical  axis.  As  the  acoustic  wavelength  changes  for 
different  rf  drive  frequencies,  the  singularity  function 
is  always  a  circle  centered  on  the  optical  axis,  but  its 
radius  varies  inversely  with  the  wavelength. 

To  complete  the  analysis  for  isotropic  crystals,  we 
now  consider  the  effect  of  the  Bragg  angle  matching 
conditions  which  have  not  been  explicitly  included  in 
the  analysis.  To  maximize  the  amplitude  of  the  dif¬ 
fracted  light,  the  input  illumination  must  be  at  the 
Bragg  angle  (Og  =  A/2A).  If  we  illuminate  the  cell  at 
the  Bragg  angle  with  respect  to  the  y-z  plane,  the  dif¬ 
fracted  light  will  be  maximized  in  the  direction  of  y  = 
;r,  and  the  singularity  occurs  at  a  =  —A/ A.  Since  the 
position  of  the  singularity  is  measured  with  respect  to 
the  point  at  which  the  undiffracted  light  is  focused 
(defined  here  as  the  optical  axis),  we  get  the  classical 
result  that  the  angle  between  the  diffracted  and  un- 
diffracted  light  is  twice  the  Bragg  angle. 

The  singularity  must  also  be  convolved  in  the  radial 
direction  with  W(<r),  which  is  the  Fourier  transform  of 
recti  r/L  —  Vo'  exp(—  er).  If  the  attenuation  is  small  over 


the  aperture  of  the  cell,  Wit r)  can  be  approximated  by 
sind^L/A).  Finally,  we  account  for  the  finite  height  of 
the  transducer  through  a  multiplication  bv  sinc(  iiH! 
A). 

In  an  anisotropic  crystal,  the  phase  velocity  of  the 
acoustic  wave  varies  as  a  function  of  the  direction  of 
propagation  and  the  orientation  of  the  crystal  relative 
to  the  transducer.  Waterman14  has  developed  rela¬ 
tionships  for  the  velocities  and  the  relative  changes  in 
velocity  for  waves  propagating  in  directions  close  to  the 
principal  axes  of  crystals.  He  shows,  for  example,  that 
velocity  u  in  the  [001]  propagation  direction  for  a  te¬ 
tragonal  crystal  is 

t,’  =  <C.J;l/p)1/2,  (111 

and  the  relative  change  in  velocity  in  a  direction  close 
to  the  direction  of  propagation  is 

At)  »-  (C|3  +  C.m)(2c44  +  C|;)  ~  C:l:l)  ( 

<•'  2ci:i  <c  a  -  c44) 

where  ctJ  are  the  elastic  constants,  p  is  the  density  of  the 
crystal,  and  <p  is  shown  in  Fig.  2.  For  small  values  of  <t> 
we  can  write  Eq.  (12)  as  Ac  =  —s<t>2v,  so  that  the  phase 
velocity  at  an  angle  <t>  is 


=  u(l  -  s«>2).  (13) 

We  use  a  negative  sign  in  the  relationship  for  A v  to  be 
consistent  with  Papadakis5  and  Cohen.8  Since  the 
acoustic  wavelength  is  A  =  2irv/d>,  where  u>  is  the 
acoustic  frequency,  the  wavelength  dependence  on  <t>  has 
the  same  form  as  Eq.  (13). 

We  are  now  in  a  position  to  see  how  the  solution  to  the 
isotropic  case  can  be  modified  to  account  for  the  an¬ 
isotropy  of  certain  materials.  The  coupling  of  the  op¬ 
tical  and  acoustic  wave  fronts  requires  that  we  first  Find 
the  acoustic  pressure  field  within  the  cell.  In  Fig.  3  we 
show  the  acoustic  K  vector  associated  with  the  phase 
velocity  of  a  propagating  wave.  The  magnitude  of  this 
vector  at  an  angle  <t>  is  given  by 


KW  = 


2t 


.4(1—  s<t>2) 


(14) 


The  acoustic  energy  in  an  anisotropic  material  propa¬ 
gates  at  an  angle  x-Q  +  'P-  Auld15  derives  a  relation¬ 
ship  between  the  magnitudes  of  the  phase  velocity  and 
the  energy  velocity  in  which  ce(x)  =  o(tf>)/ cosi/,  where 
\p  is  the  angle  between  the  phase  vector  and  the  energy 
(Poynting)  vector.  Furthermore,  Papadakis5  shows 
that  p  =  —2 s<j>  so  that 


which  is  valid  to  terms  of  the  order  of  x2-  The  magni¬ 
tude  of  the  acoustic  K  vector  associated  with  the  energy- 
flow  is  then 


a;  2jt 

K.ax)  = - =  — ; - - - ’• 

\)  .V(l  —  s\~/(  1  —  2s  1 1 

The  scalar  product  of  the  Poynting  vector  and  a  dis¬ 
tance  vector  defines  the  pressure  at  each  point  within 
the  cell.  The  scalar  product  is 
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K.T»Ylx+;y2/2*(1-2*l-  (17) 

where  we  have  used  the  relationship  that  x  -  arc- 
tan(y/x)  ~  y/x.  The  Fourier  transform  of  the  optical 
wave  front  in  the  anisotropic  case  is  then 

Ma&  *  X,  S  vfe eip  H  t  [i + y2/2ja  -  24,|I 

f  2  r 

X  exp  \—j  —  (ox  +  (3y)  exp (-tx)dxdy,  (18) 

where  we  have  replaced  r  by  x  in  the  denominator  and 
in  the  attenuation  factor  because  we  require  that  x  be 
small.  This  integral  is  considerably  more  difficult  to 
solve  than  Eq.  (2);  the  solution  is  given  in  the  Appendix 
where  we  find  that  the  singularity  function  has  the 
parabolic  form 


A(q  +  A/A) 
APfe  -  .») 


(19) 


Although  Eq.  (18)  is  the  more  complete  diffraction 
integral  because  it  accounts  for  the  pressure  field  via  the 
Poynting  vector,  the  same  singularity  function  can  be 
obtained  in  a  simpler  fashion  for  the  special  case  of  a 
continuous  wave  signal.  In  this  case  we  cannot  physi¬ 
cally  distinguish  the  direction  of  the  Poynting  vector 
since  we  assume  that  the  infinitesimal  transducer  ra¬ 
diates  energy  in  all  directions.  We  can  then  find  the 
singularity  function  more  directly  from  the  isotropic 
solution  by  using  the  scaling  property  of  the  Fourier 
transform.  If  a  2-D  aperture  function  is  extended  by 
a  scale  factor  in  a  particular  direction,  the  diffraction 
pattern  contracts  in  the  same  direction,  and  the  am¬ 
plitude  is  adjusted  by  the  same  scale  factor.16  From  Eq. 
(14)  it  ,s  easily  seen  that  the  acoustic  wavelength  is 
scaled  by  a  factor  (1  -  s<p2)  so  that  the  spatial  frequency 
<r  for  the  isotropic  case  is  scaled  by  (1  —  s02)-1  in  the 
direction  <b.  We  then  have  that  rr{<i>)  =  (,\/A)(l  —  s02)-1 
so  that 


a  ■ - —  ■  3  =* - —  ■  120) 

A(  1  —  s02)  All  -  s<t>2) 

The  amplitude  A(a)  as  given  by  Eq.  (8)  must  also  be 
scaled  by  (1  -  sdr)  to  obtain  the  result  for  the  aniso¬ 
tropic  case.  Since  the  physical  boundaries  of  the  ap¬ 
erture  do  not  change,  the  weighting  functions  are  the 
same  as  those  for  the  isotropic  case. 


The  parabolic  form  of  the  singularity  function  can 
also  be  obtained  from  Eqs.  (20)  if  we  note  that  3  *  \<p/A 
and  use  this  value  for  3  in  the  expression  for  a.  For  the 
negative  diffracted  order,  we  have 

a  *>  —  [1  —  (s  -  lh>02]  *  —  [1  -  (s  —  V2)(A/S/X)2), 

A  A 

from  which  we  get  the  result  that 

„  Ma  +  A/A) 

**1^7 

in  agreement  with  Eq.  (19). 

In  Fig.  4  we  show  the  singularity  function  of  aniso¬ 
tropic  materials  for  several  values  of  the  parameter  s. 
For  s  =  0,  the  singularity  function  is  a  sector  of  a  circle. 
For  s  <  0  the  singularity  function  curves  more  rapidly 
toward  the  optical  axis,  whereas  for  s  >  0,  the  singularity 
function  curves  less  rapidly  toward  the  source.  The 
horizontal  difference  between  the  spatial  frequency  at 
an  angle  0  and  that  at  0  =  0  is  obtained  from  Eq. 
(20): 

A <x(<t>)  *  ao  -  =*  U  -  V2)o2A/A.  (21) 

From  this  relationship  we  see  that  s  =  V2  results  in  a 
singularity  function  that  is  a  straight  vertical  line 
passing  through  a  =  -X/A.  Papadakis5  gives  some 
limits  on  the  value  of  s  in  terms  of  the  elastic  constants 
and  argues  that  s  must  <  '/2.  The  values  of  s  as  given  by 
Papadakis  and  by  Waterman  range  from  +0.38  for  KC1 
to  -5.23  for  Zn.  Cohen8  observed  that  a  crystal  for 
which  s  =  V2  is  one  in  which  the  acoustic  beam  diver¬ 
gence  is  zero.  Such  a  crystal  would  produce  a  self-col- 
limating  acoustic  beam,  and  the  diffrr  :tion  pattern 
would  be  identical  to  that  produced  by  a  2-D  rectan¬ 
gular  aperture  having  a  uniform  phase  weighting. 

Equation  (21)  also  shows  that  the  magnitude  of 
Aa(0)  can  be  large  if  s  is  highly  negative,  which  is  the 
case  for  TeO>  operated  in  the  slow-shear  mode.  Fur¬ 
thermore,  Aa(0)  gives  a  direct  measure  of  s  from  ex¬ 
perimentally  generated  diffraction  patterns  as  we  shall 
show  in  the  next  section. 


III.  Experimental  Results 

Figure  2  shows  the  experimental  setup  used  to  gen¬ 
erate  the  diffraction  patterns.  Light  from  a  He-Me 
laser  is  collimated  so  that  the  Bragg  cell  is  illuminated 
at  the  Bragg  angle.  Lens  Ln  produces  the  2-D  Fourier 
transform  of  the  fully  illuminated  Bragg  cell  at  plane 
P_>-  Figure  5  shows  the  resulting  diffraction  pattern  for 
a  PbMo04  cell  driven  by  a  single  frequency  at  400  MHz. 
The  intensity  I(a,3)  follows  the  singularity  function 
whose  coordinates  are  given  by  Eq.  (20).  The  intensity- 
weighting  in  the  3  direction  is  sin c2(tf  d/X>,  and  the  in¬ 
tensity  weighting  normal  to  the  singularity  function  is 
sinc-(L<r/X).  An  additional  weighting  in  the  horizontal 
direction  is  due  to  the  Bragg  angle  mismatch  conditions: 
it  is  not  obvious  unless  the  value  of  s  is  highly  negative. 
At  3  =  0  the  singularity  function  occurs  at  o  =  -  A 
relative  to  the  focal  point  of  the  undiffracted  light, 
where  A/A  is  twice  the  Bragg  angle. 
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Pig.  5.  Diffraction  pattern  for  a  fully  illuminated  lead  molybdate 
Bragg  cell  driven  at  400  MHz.  Undiffracted  light  is  attenuated  by 
neutral  density  filter. 
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Fig.  6.  Diffraction  patterns  for  a  tellurium  dioxide  Bragg  cell  driven 
at  500  and  300  MHz  (double  exposure). 


_ Table  l.  Beetle  Consents  (X  IQ11  W/m») _ 

PbMoO17 _ TeO;18 

cu  1.092  053T 

c  12  0.683  0.486 

ci3  0.528  0.212 

C33  0.917  1.085 

c«4  0.267  0.244 

egg  0.337  0.562 


Figure  6  is  a  photograph  obtained  by  a  double  expo¬ 
sure  of  diffraction  patterns  for  TeC>2  excited  at  500  and 
300  MHz.  The  features  noted  above  are  evident,  as  well 
as  the  fact  that  the  location  of  /(a,/3)  is  a  function  of  A 
and  that  the  singularity  function  becomes  less  highly 
curved  as  the  acoustic  drive  frequency  increases.  In  the 
limit  of  very  high  frequencies,  or  a  large  ratio  of  A/ A,  the 
diffraction  pattern  more  nearly  approximates  an  or¬ 
thogonal  set  of  sine2  functions.  From  this  figure  we  can 
also  deduce  that  the  TeO?  cell  produces  a  singularity 
function  that  opens  more  slowly  than  that  for  PbMoQ4. 
At  400  MHz,  the  acoustic  wavelength  is  9.1  fim  for 
PbMoO.*  and  10.5  *im  for  TeOj,  which  accounts  for  part 
of  the  decreased  curvature.  The  remainder  is  due  to 
different  values  for  the  relative  changes  in  velocity.  By 
measuring  the  locus  for  each  singularity  and  by  using 
Eq.  ( 19),  we  find  that  s  =  —0.168  for  PbMoO*  and  s  = 
0.24  for  TeOo.  These  values  will  now  be  compared  with 
those  obtained  from  the  elastic  constants.  From  Eq. 
U2)  we  note  thar 


— ) C  1 ,1  +  C33)(2c»4  +  Cl3  -  Cm)  ( 

2c33<cn  -  c«) 

In  Table  I  we  give  the  values  of  the  elastic  coefficients 
for  PbMoO,*  and  Te02  from  which  we  calculate  that  s 
=  —0.176  for  PbMo04  and  s  =  0.274  for  TeOo.  Our 
measured  data  for  s  are  in  reasonable  agreement  with 
those  calculated  from  Eq.  (22),  given  the  degree  of  ac¬ 
curacy  in  measuring  the  c,>. 19 

We  also  tested  a  Bragg  cell  fabricated  from  fused 
quartz,  which  is  isotropic,  and  found  the  singularity 
function  to  be  a  circle  as  predicted  by  Eq.  (10)  or  (20). 
Since  the  rate  at  which  the  singularity  function  deviates 
from  a  straight  line  is  proportional  to  (V2  —  s),  we  find 
that  the  rate  is  lowest  for  Te02  and  highest  for  PbMo04; 
the  rate  for  fused  quartz  is  nearly  halfway  between  these 
two  values. 

The  diffraction  patterns  produced  by  acoustooptic 
cells  operated  in  the  Bragg  mode  are  similar  to  those 
obtained  by  Shaeffer  and  Bergmann20  for  cells  operated 
in  the  Raman-Nath  mode.  In  the  examples  they  give, 
as  well  as  those  illustrated  in  Refs.  21  and  22,  the  dif¬ 
fraction  patterns  are  symmetric  about  the  focal  point 
of  the  undiffracted  light  because  drive  frequencies  are 
low.  These  patterns  are  also  generally  more  compli¬ 
cated  because  both  shear  and  longitudinal  waves 
propagate  within  the  cell. 

The  significance  of  the  fact  that  the  diffraction  pat¬ 
tern  for  Bragg  cells  lies  along  a  curved  singularity 
function  instead  of  a  straight  line  depends  on  the  ap¬ 
plication.  The  Fourier  transform  of  signals  introduced 
by  Bragg  cells  is  used  extensively  in  spectrum  analysis 
and  in  spatial  filtering  applications.  In  these  applica¬ 
tions  the  Bragg  cell  is  often  illuminated  by  a  line  source 
generated  by  inserting  a  cylindrical  lens  between  L  *  and 
the  Bragg  cell  (see  Fig.  2)  as  well  as  between  the  Bragg 
cell  and  lens  L2.  In  this  fashion,  a  1-D  Fourier  trans¬ 
form  is  generated  with  the  line  source  being  imaged  in 
the  orthogonal  direction.  In  some  configurations, 
however,  it  may  not  be  possible  to  use  a  line  illumination 
because  two  Bragg  cells  are  used  in  series,23  or  they  are 
used  in  a  Rayleigh  interferometer  wherein  a  single  lens 
is  used  to  ensure  that  the  Fourier  transforms  overlap. 
For  such  configurations  we  want  to  know  the  effect  of 
integrating  (detecting)  light  in  the  vertical  direction. 

In  Fig.  7  we  see  the  results  of  integrating  l(a,0)  in  the 
0  direction  to  the  tenth  zeros  on  either  side  of  the  hor¬ 
izontal  axis.  Since  Aa  as  given  by  Eq.  (21)  gives  the 
deviation  of  the  singularity  function  from  a  straight  line, 
we  define  a  parameter  D  =  (V 2  —  s)A/A;  the  integrated 
intensity  is  shown  for  D  equal  to  0, 1/16, 1/32,  and  1/64. 
The  result  for  D  =  0  represents  the  integrated  intensity 
of  the  diffraction  pattern  of  a  hypothetical  self-colli¬ 
mating  Bragg  cell  (s  =  V2)  and  forms  a  basis  for  com¬ 
paring  the  results  for  other  values  of  D.  We  compare 
the  results  for  D  =  1/64  with  the  ideal  case,  since  it 
shows  the  greatest  departure  from  the  ideal.  First,  the 
sidelobes  are  asymmetric,  being  higher  in  the  direction 
of  the  optical  axis,  and  the  nulls  are  less  well  formed. 
Second,  the  main  lobe  width  is  3%  wider  and  is  shifted 
toward  the  optical  axis  by  4%  of  the  m.->in  iooe  width  as 
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Fig.  7.  Normalized  intensity  of  diffraction  pattern  integrated  in  the 
vertical  direction.  The  result  for  D  *  0  is  a  sinc-’(a)  function  centered 
at  -A/A. 


Fig.  8.  Details  of  the  intensity  of  the  central  lobe  of  diffraction 
pattern. 


can  be  9<!en  more  clearly  in  Fig.  8.  Third,  the  peak  in¬ 
tensity  u  reduced  by  13%  because  the  diffraction  pat¬ 
tern  is  h.  ghly  curved  toward  the  source.  Note  that  the 
results  for  an  isotropic  material  (D  *  1/16)  show  a 
similar  trend;  the  best  results  are  obtained  for  an  an¬ 
isotropic  material  having  a  value  of  s  as  close  to  as 
possible. 

In  Fig.  9  we  show  the  results  when  l(a,0)  is  integrated 
to  the  filth  zeros  on  either  side  of  the  horizontal  axis. 
The  same  general  features  as  described  above  are  evi¬ 
dent  except  that  the  degradation  from  the  ideal  is  not 
as  great.  Experimental  data  obtained  from  a  TeC>2  cell 
operating  at  290  MHz  (equivalent  to  D  =  1/32)  is  also 
plotted  in  Fig.  9.  Although  the  overall  level  of  the 
sidelobes  is  higher  for  the  experimental  data  due  to 
scattered  light,  the  general  form  of  the  data  is  in  good 
agreement  with  the  theory. 

The  best  performance  is  obtained  when  we  integrate 
only  to  the  first  zeros  in  the  vertical  direction.  The 
results,  shown  in  Fig.  10  for  the  sidelobes  extending 


toward  the  optical  axis,  suggest  that  the  Fourier  trans¬ 
form  should  be  truncated  at  the  first  zeros  by  the  de¬ 
tector  to  avoid  an  increase  in  the  sidelobe  levels.  The 
penalty  for  truncation  is  a  loss  of  ~10%  of  the  available 
optical  power. 

IV.  Bragg  Call  Illumination 

In  some  applications,  the  Bragg  cell  is  illuminated  by 
a  line  source  generated  by  the  use  of  a  cylindrical  lens 
arrangement  as  described  in  Sec.  III.  We  now  compare 
the  effects  of  1-D  illumination  with  that  of  2-D  illumi¬ 
nation  on  the  diffraction  pattern.  In  this  case,  the 
diffraction  pattern  is  understood  to  exist  in  the  hori¬ 
zontal  direction,  whereas  the  light  distribution  in  the 
vertical  direction  is  the  image  of  the  line  source. 

We  use  the  same  model  for  acoustic  wave  propagation 
as  shown  in  Fig.  2.  We  now  cannot  account  for  the  fi¬ 
nite  height  of  the  transducer  through  a  multiplication 
of  the  singularity  function  by  sinc(///3/A)  in  the  Fourier 
plane;  instead  we  must  calculate  the  effect  of  the  finite 


Fig.  9.  Theoretical  and  experimental  intensity  of  diffrac  tion  pattern 
for  TeOs  <290  MHz)  integrated  to  the  fifth  zetus. 


Fig.  It).  Intensity  of  diffraction  patterns  when  integrated  in  the 
vertical  direction  to  only  the  first  zero.  Normalized  so  that  the  in¬ 
tensity  is  I  t)  at  o  =  t)  for  D  =  t). 
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Fig.  11.  Amplitude  (left-hand  scale!  and  phase  (right-hand  scale 
in  radians)  along  the  v  •  0  axis  due  to  the  entire  transducer. 


transducer  at  the  plane  of  the  Bragg  cell.  The  contri¬ 
bution  from  an  elemental  area  of  the  transducer  for  the 
anisotropic  case  is  obtained  from  Eq.  (18): 

a(xo>)  •  -  expl-j  —  (x  +  y2/2x(l  -  2a)|}  exp(-«x). 
vjAx  l  A  1 

We  convolve  this  function  with  rect (y/H)  to  obtain  the 
total  amplitude  and  phase  function  over  the  exit  face 
of  the  Bragg  cell.  The  convolution  integral  is  then 

a(x,y>»  f  )Az=  «xp  f-j  — ^ rect (u/H)du,  (23) 
J—y/JKx  l  Ax(l-2*> 

where  we  have  neglected  the  attenuation  factor  f. 

As  before,  we  assume  that  the  width  of  the  transducer 
along  the  optical  axis  is  large  so  that  end  effects  are 
negligible.  This  assumption  may  seem  to  be  less  valid 
than  before  because  the  light  is  now  focused  within  the 
cell.  The  transducer  for  the  PbMo04  cell  is  150  jxm 
high  and  600  inn  wide.  If  a  Gaussian  beam  is  focused 
within  the  cell  so  that  the  beam  waist  height  is  *  40  nm, 
the  beam  waist  does  not  grow  by  more  than  4%  at  the 
entrance  and  exit  faces  of  the  cell.24  The  line  illumi¬ 
nation  is,  therefore,  more  nearly  equivalent  to  a  colli¬ 
mated  beam  than  it  is  to  a  convergent  or  divergent 
beam. 

The  rect (u/H)  function  in  Eq.  (23)  can  be  removed 
by  restricting  the  limits  of  integration  to  |  u  |  <  H/2  so 
that2S 

alx.y)  »  'tlCfd  -v)v<?/x|  +)-S[(l  —  yly  <//x|l 

+  'yC|(  1  +  v  )y  q/x  |  +  jS[{  1  +  y  ly  q/x  |l,  (24) 

where  C(>)  and  S(-)  are  the  Fresnel  sine  and  cosine  in¬ 
tegrals.  and  q  *  H2l 2A(1  -  2 s)L.  The  variables  in  Eq. 
(24)  have  been  normalized  so  that  |y  |  <  1  and  0  <  x  < 
1;  this  result  shows  that  the  amplitude  transmittance 
along  the  Bragg  cell  at  any  value  of  y  is  a  combination 
of  Fresnel  integrals.  The  amplitude  and  phase  along 
the  horizontal  axis  (y  *  0)  are  shown  in  Fig.  11  for  a 
PbMo04  Bragg  cell  driven  at  400  MHz  with  H  *  150  inn 
and  L  ■  2.6  mm.  Since  the  Fresnel  integrals  are  func¬ 
tions  of  v  1/x.  both  the  amplitude  and  phase  oscillate 
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rapidly  at  small  values  of  x.  The  frequency  of  oscilla¬ 
tion  decreases  as  the  distance  from  the  transducer  in¬ 
creases,  and,  for  large  values  of  x,  the  amplitude  decays 
slowly.  The  maximum  dip  in  the  amplitude  occurs  at 
xp;  from  Eq.  (24)  we  find  that  this  dip  occurs  when  the 
argument  of  the  Fresnel  integral  is  *  1.88.  Therefore, 
we  find  that 


For  an  isotropic  medium,  this  value  of  x0  agrees  rea¬ 
sonably  well  with  that  predicted  by  Ingenito  and  Cooka 
who  give  a  value  of  x,)  =  H-/9AL. 

We  now  compute  the  Fourier  transform  of  a(x,>  )  to 
obtain  the  diffraction  pattern  in  the  horizontal  direc¬ 
tion.  For  convenience,  we  normalize  the  angular  spatial 
frequency  so  that  i*i  =  (tv  +  \/A)L/X: 

A'(tri,yl  »  C  aix.y >  exp(-;2»«ix  )dx  i26l 

The  final  step  is  to  integrate  the  intensity  of  the  Fourier 
transform  over  y  to  obtain  the  diffraction  pattern  that 
will  be  measured  by  a  photodetector.  Figure  12  shows 
the  normalized  magnitude  squared  of  F(«  i)  for 
PbMoO-t  as  well  as  the  diffraction  pattern  for  the  ideal 
case  of  a  self-collimating  material.  We  see  that  the 
central  lobe  is  slightly  broader  due  to  a  combination  of 
the  amplitude  and  phase  weighting  of  a(x,y).  The 
entire  diffraction  pattern  is  shifted  by  a  small  amount; 
from  Fig.  1 1  we  see  that  the  linear  component  of  the 
phase  over  the  region  0.2  <  x  <  1  is  ~0.5  rad  so  that  the 
pattern  should  shift  by  nearly  tt/10  (i.e.,  ~10%  of  the 
angular  distance  between  adjacent  nulls).  In  compar¬ 
ing  fhe  results  in  Fig.  12  with  those  in  Fig.  7,  we  see  that 
the  aidelobes  toward  the  optical  axis  are  lower,  instead 
of  higher,  than  those  due  to  a  uniformly  illuminated 
aperture.  The  use  of  line  illumination,  therefore,  causes 
the  sidelobe  levels  to  be  more  nearly  equal  to  those 
produced  ty  a  uniformly  illuminated  aperture. 

From  Eq.  (24)  and  Fig.  11,  we  see  that  the  amplitude 
transmittance  fluctuates  rapidly  in  the  region  x  <  x(). 


^SPATIAL  FREQUENCY 


Fii!  12.  Intensitv  of  diffraction  pattern  when  the  Brags  cell  is  illu¬ 
minated  by  a  line  source. 


u 


We  might  expect,  therefore,  that  the  Fourier  transform 
would  be  more  nearly  ideal  if  we  were  to  truncate  the 
illumination  so  that  this  region  does  not  contribute  to 
the  Fourier  integral.  We  calculated  the  Fourier 
transform  of  a(x  ,0)  only  over  the  range  i0  5  i  <  1  and 
noted  some  differences  in  the  structure  of  the  sidelobes 
extending  toward  the  optical  axis.  The  major  change 
is  in  the  value  of  the  second  sidelobe  and  the  depth  of 
the  second  null  in  the  Fourier  transform.  As  we  in¬ 
creased  the  range  of  y  over  which  we  integrate  at  the 
Fourier  plane,  we  found  the  differences  to  diminish  so 
that  truncation  of  the  oscillating  portion  of  a(x,y)  has 
less  effect. 

IV.  Summary  and  Conclusions 

The  diffraction  pattern  of  a  uniformly  illuminated 
Bragg  cell  lies  along  a  singularity  function  that  is  de¬ 
termined  by  the  phase  velocity  of  the  interaction  ma¬ 
terial.  For  an  isotropic  material,  this  curve  is  a  circle 
whose  origin  is  at  the  focal  point  of  the  undiffracted 
beam  and  whose  radius  is  twice  the  Bragg  angle.  For 
an  anisotropic  interaction  medium  we  used  Waterman's 
approximation  to  the  phase  velocity  profile  and  found 
that  the  singularity  function  is  a  parabola  having  its 
vertex  at  twice  the  Bragg  angle.  Experimental  mea¬ 
surements  of  s  are  in  close  agreement  with  those  ob¬ 
tained  from  the  elastic  constants. 

In  signal  processing  applications  where  the  Bragg  cell 
must  be  uniformly  illuminated,  the  curved  singularity 
function  causes  the  sidelobe  levels  in  the  direction  of  the 
source  to  increase  when  the  diffraction  pattern  is  inte¬ 
grated  along  a  vertical  line  instead  of  along  the  locus  of 
the  curve.  The  severity  of  this  problem  can  be  reduced 
by  restricting  the  region  of  integration  in  the  vertical 
direction;  the  permissible  region  of  integration  is  a 
function  of  X,  A,  and  s.  Line  illumination  of  the  Bragg 
cell  can  be  used  to  further  reduce  the  deterioration  of 
the  diffraction  pattern.  Some  residual  aberrations  still 
exist  under  these  conditions,  but  their  effects  are  small. 
A  final  technique  for  improving  the  diffraction  pattern 
is  to  use  an  aperture  stop  to  remove  the  contributions 
of  rapid  acoustic  wavs  front  oscillations  near  the 
transducer. 
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programs.  Special  thanks  to  S.  J.  Adhav  for  assistance 
in  the  experimental  work.  This  activity  was  performed 
under  a  grant  from  the  U.S.  Army  Research  Office. 

Appendix 

From  Eq.  1 18)  we  have  that  the  effective  amplitude 
distribution  at  the  Bragg  cell  is 

1  2t,  ,  ,1 

a i  x.v)  «*  ■  ■  —  exp  '  1  —  It  +  v-/2r(  1  -  2s )| )  expi  —<x  i. 

\  >Ajc  I  A  I 

and  the  region  of  integration  is  a  wedge  having  its  apex 
at  x  *  0  so  that  0  <  x  <  L  and  |y  j  <  0t)x ,  where  ©»  de¬ 
termines  the  upper  and  lower  boundaries  of  the  wedge. 
We  first  perform  the  integration  or,  y: 


/s  4o* 

(jtx.J >■  I  exp|-;2iry-1/2Axt  1  -  2s  l|  exp<-j2wd\/\)dv. 
which  can  be  solved  by  using  the  relationship  that26 
J'  exp<-ct2  -  2dt)dt  “  yi  -expi d-/cl 

x  erfl  \  ct  +  d/\  <■  l;  c  0 

If  we  make  the  associations  and  definitions  that 
c  “  ;x/Ax)I  -  2s I, 

d  ■  jxdl A,  ( A1  * 

e  »  ;t/A(  1  -  2s), 

we  have  that 


Cixji)  ■  lk\  j\x  exp \dttie\ 

X  lerl|(0,iv  «>  +  d/\  e  lv  x  \ 

-  erf|(-®o\  f  +  d/\  e)vT||.  <A2) 


We  now  account  for  the  finite  transducer  height  H  by 
multiplying  G(x,ti)  by  sind/fd/A).  The  integration  on 
x  now  becomes 


Aitt.tl)  =  '/2  sinc(H|3/A)  ^  erfjtipov  e  +  d!\  e)\  : 


x  exp|(d2/e  —  <)x|  exp 


(A3) 


minus  a  similar  integral  in  which  <t> o  is  replaced  by  —  4o- 
We  can  extend  the  upper  limit  to  infinity  if  we  multiply 
the  integrand  by  rect (r/L  —  V2)  and  account  for  the  fi¬ 
nite  limit  through  a  convolution  in  a  with  sinc(Lo/A). 
We  now  use  the  relationship  that26 


P  exp(-ft)  erfv  gtdt  »  -  ■»/  — — ;  Re<f  +  g)  >  0.  (A4) 
Ju  f  *  f  +  g 

The  presence  of  the  attenuation  factor  (  ensures  that 
Re(f  +  g)  >  0  so  that  Eq.  ( A4)  is  valid.  We  must  exer¬ 
cise  some  care  in  applying  Eq.  (A4)  in  evaluating  Eq. 
(A2)  because  the  signs  of  the  arguments  of  the  error 
function  change  as  a  function  of  /3.  Consider  the  in¬ 
terval  for  which  \d/\fe\  <  4>  )\re:  from  Eq.  (Al)  we  find 
that  this  corresponds  to  1 0|  <  ®oX/ A 1 1  —  2s ).  For  this 
range  on  /3,  the  sign  of  the  argument  of  the  second  error 
function  in  Eq.  ( A2)  is  negat  ive,  and  we  can  use  the  re¬ 
lationship  that  erf'-t)  *  -erf(t).  We  now  make  the 
associations  that 


/  *  J  —  «  4 - -  d2/f  +  f. 


\  fix  *  'i*n\  ^ 


■d.\ 


( A5I 


v  H-2  m  ^>o\  ••  -  d/v  e. 

where  we  use  g\  in  Eq.  (A3)  and  g>  in  the  companion 
integral.  We  now  have  that 

(1  ti  I  «»  .  ' 

-  ■,  — i —  +  -  \  |  •  (AH) 

/  V  ,-rx,  /V  /♦*./ 

As  explained  in  the  text,  we  expect  .4(a.d)  to  be  a 
singularity  function  in  a  and  d-  Such  a  function  is 
obtained  in  the  limit  as  f  tends  toward  zero.  From  Eqs. 
(A5‘  and  (Al)  we  find  that  >  *  0  yields 
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„  Mi*  +  V  A) 

(J2-— - -  (AT) 

A(','2  -  s) 

in  the  limit  as  f  —  0.  Letting  e  -»  0  is  equivalent  to  the 
requirement  that  the  acoustic  attenuation  be  very  small; 
in  a  practical  sense,  if  f  is  small  over  the  aperture  length 
L,  use  of  the  singularity  function  is  valid.  Equation 
(A6)  now  becomes 


A(a,3)  •  **[  sinc(H0A)3 


A(a  +  A/ A) 
AC'j-s) 


(A8) 


which  must  be  convolved  with  sinc(LaA)  to  account  for 
the  finite  length  of  the  Bragg  cell.  The  resulting  in¬ 
tensity  of  the  Fourier  diffraction  pattern  is  then 


Harf) 


sine2 


sinc2(Hd/\), 


(A9> 


where  f{a,j3)  is  the  argument  of  the  5  function  in  Eq. 
(A8). 

Equation  (7)  is  an  important  result  since  it  embodies 
many  features  of  the  diffraction  pattern.  First,  it  is  the 
standard  form  of  a  parabola  whose  vertex  is  located  at 
a  *  —A/ A.  The  position  of  the  vertex  is  inversely  re¬ 
lated  to  the  acoustic  wavelength  and  directly  related 
acoustic  frequency.  Second,  the  ratio  of  A/A  is  twice 
the  Bragg  angle;  that  is,  at  any  acoustic  frequency  u>,  the 
angular  displacement  of  I(a,j3)  at  0  *  0  is  at  2 dB  with 
respect  to  the  undiffracted  light.  If  we  had  illuminated 
the  cell  at  the  Bragg  angle  with  respect  to  the  optical 
axis,  which  requires  a  factor  exp(-/2ir0ex/A)  in  the 
integrands  of  Eqs.  (2)  and  ( 18)  in  the  text,  the  undif¬ 
fracted  and  diffracted  light  would  propagate  at  angles 
of  ±0r  with  respect  to  the  normal  to  the  Bragg  cell. 
Third,  the  latus  rectum  of  the  parabola  is  equal  to 
A/A(V2  -  s);  the  value  of  s  then  determines  the  degree 
of  curvature  of  the  diffraction  pattern  in  the  vicinity  of 
a  ■  —A/A.  Fourth,  if  s  ■  0,  the  crystal  i8  isotropic  so 
that  the  singularity  equation  becomes 


J2»  —  (o  +  A/A),  (A10) 

A 

which  shows  that  the  focus  of  the  parabola  is  midway 
between  the  vertex  and  the  axis  defined  by  the  undif¬ 
fracted  light. 
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ABSTRACT 

In  some  applications,  the  performance  of  multichannel  Bragg  cells  is 
compromised  by  the  spreading  of  the  aconstic  naves  as  they  propagate;  the 
spreading  causes  the  signals  in  the  channels  to  overlap.  The  overlapping 
oan  be  signif ioantly  reduced  by  a  spatial  filter  in  a  Fourier/ image  plane 
the  spatial  filter  is  shown  to  be  a  cylindrical  lens  whose  power  is  a 
function  of  the  distanoe  from  the  transducer.  The  effects  of  changes  in  the 
drive  frequency  as  well  as  those  of  displacements  of  the  filter  are 
calculated.  The  reduction  in  the  modulation  transfer  funotion  as  a  function 
of  propagation  distance  is  calculated  sad  some  bounds  on  the  time-bandwidth 
product  and  the  number  of  channels  are  derived.  In  general,  the  overall 
performance  can  be  Improved  by  increasing  the  center  frequency  of  the  Bragg 
cell  while  keeping  the  bandwidth  fixed. 
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la  sea*  signal  processing  applications,  multichannel  Bragg  calls  are 

required  to  handle  wideband  parallel  signals.  One  exaaple  is  that  of 

processing  signals  from  a  phased  array  antenna  wherein  we  associate  each 

channel  of  the  Bragg  cell  with  an  antenna  eleaent.  Multichannel  Bragg  cells 
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were  deweloped  in  the  early  1960’s  by  Laabert  and  his  associates  *  ,  and 
cells  with  as  aany  as  128  channels  haws  been  described3.  As  the  channels 
beeoae  aore  densely  packed,  the  transducer  heights  are  reduced  causing  the 
acoustic  energy  to  spread  ower  larger  angles  as  it  propagates  through  the 
interaction  aaterial.  The  acoustic  waves  froa  adjacent  channels  therefore 
overlap  after  a  short  propagation  distance  so  that  the  utility  of  the  device 
aay  be  less  than  that  desired. 

Bzpcriacnts  with  a  128-channel  Bragg  cell  constructed  froa  SF-8  dense 
flint  glass  and  with  a  32-channel  device  constructed  froa  TeO,  clearly  show 
the  acoustic  spreading  phenaaena.  In  the  high  speed  recording  applications 
for  which  these  devices  were  developed3,  only  the  region  near  the  transducer 
was  used  to  aodulate  the  light  besa.  Ve  also  noted,  however,  that  even 
after  the  acoustic  beaas  had  overlapped  coaqpletely,  the  diffraoted  light 
pattern  could  be  separated  into  disorete  channels  to  reseable  the  pattern 
produced  near  tha  transducer.  This  was  aceoaplished  by  foraing  a  Schlieren 
image  at  a  plane  displaced  axially  froa  that  of  the  Bragg  oell  itself.  This 
result,  in  turn,  suggested  that  the  acoustic  spreading,  in  soma  sense, 
introduces  the  equivalent  of  a  variable  foousing  power  whioh  is  a  linear 
function  of  the  propagation  distance.  If  this  focusing  setion  could  be 
coapcnsated  so  that,  at  soae  plane  in  an  optical  systea,  the  diffracted 
light  froa  each  channel  is  confined  to  that  channel,  a  aore  widely  useful 
multichannel  Bragg  cell  would  result. 

It  is  not  sufficient,  however,  to  confine  the  diffrseted  light  to  the 
channel  without  giving  consideration  to  the  is^act  that  the  coapensation  has 
on  the  aodulation  transfer  function.  As  the  scoustic  energy  propagstes,  the 
inforaation  diverges  so  that  ooapensation  causes  it  to  be  spread  in  the 
direction  of  propagation.  The  result  is  s  loss  in  the  aodulation  transfer 
funotion  that  is  a  function  of  the  propagation  distance.  The  desire  to 


confine  til*  acoustic  spreading  iut  than  bt  btltscid  by  u  leoipttbl*  lost 
la  tha  aodulstion  transfer  function. 
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Tha  aodcl  that  vc  uaa  for  analysing  tha  aoonatic  apraading  ia  tha  a  am  a 

aa  that  aaad  before*  and  ia  shown  in  Fignra  1.  An  alaotrioal  aignal  drive* 

a  piaxoalaotric  tranadncar  having  height  H  and  aa  intaraotion  width  V.  Tha 

traaadaoar  laaaehaa  an  aooaatie  wav*  within  tha  Bragg  call  whioh  changes  tha 

index  of  refraction)  this,  in  torn,  causes  tha  phase  of  light  froa  a 

coherent  source  to  be  nodulated  in  spaoa  and  tine.  If  the  drive  signal  is 

an  KF  signal  at  frequency  f  ,  the  acoustic  wavelength  is  A  -  v/f  ,  where  V 

6  0  6 

is  the  velocity  of  sound  in  the  aediua. 

The  aooustie  wavefronts  propagate  ia  a  fashion  siailar  to  optioal 
wave front a  derived  froa  a  line  souroe  (H  ->  0) .  If  1  »  H.  we  can  nodal  tha 
wavefronts  as  cylindrical  sheets  which,  for  sa  isotropic  aedina.  spread  at 
an  angle  +♦  with  respect  to  the  x-axis .  The  degree  to  which  aooustie 
spreading  occurs  is  a  function  of  the  anisotropy  of  the  aedina,  the  acoustic 
wavelength  and  the  transducer  height. 
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The  first  step  is  to  calculate  the  Fourier  transfora  of  a  single 
channel  Bragg  cell  driven  as  shown  in  Figure  1.  1*  wish  to  derive  the 

traasfora  in  the  ydirectly  only,  while  we  inage  the  Bragg  cell  in  the 
x-direotion.  Ve  begin  by  considering  the  transducer  to  be  a  infinites iaal 
line  source  and  account  for  its  finite  height  by  aultiplying  the  Fourier 
trsnsfora  by  a  sinc-fnaotion.  The  fourier  integral  to  be  solved  is* 


6(x,0)  - 


i 

35  J 

*V 


exp  [  -  j  2  xy  *  /  2  i\x  ( 1-2  s )  ]  exp  ( -  j  2  n  By  /  X )  dy , 


(1) 


4,5 


where  fi  is  the  angular  spatial  frequency,  s  is  s  paraaeter  that 
characterises  the  degree  of  anisotropy  as  given  by  the  elastic  constants, 

A  is  the  wavelength  of  the  acoustic  wave  within  the  aedina,  and  X  is  the 
wavelength  of  light.  The  region  of  integration  is  over  a  wedge  whose  apex  is 


it  i  •  0  and  0  is  the  angle  at  which  the  aconstic  beam  has  its  first  nails  ; 
o 

this  region  contains  most  of  the  optical  energy.  The  far  field  value  of  $  i 

♦  «  A  (l-2a)/H.  (2) 

o 

In  reference  4,  we  solved  (1)  in  terns  of  error  functions  having  complex 
valued  arguments.  An  equivalent  result  that  is  somewhat  simpler  to  derive 
and  provides  more  physical  insight  can  be  obtained  through  the  use  of 
Fresnel  integrals.  Ve  let 

c  ■  n/Ax(l-2s) 

(3) 

d  -  np/X 


so  that 


•i)  -  -p=  f 


«xp[-j(cy*+2dy)]  <*y 


By  cosqileting  the  square  of  the  exponential  and  by  changing  variables,  we 
obtain 


G(x.p) 


/)Ax 


exp( jd ‘/c ) sine (pH/X) /exp(-ju* )du. 


where  a  ■  x^c”  +  d/  b  *  x^e"  +  d/  JT,  and  sinc(pH/X)  is  the 
o  o 

multiplicative  factor  needed  to  account  for  the  finite  height  of  the 
transducer.  The  integral  is  a  Fresnel  integral  that  further  modifies  the 
amplitude  of  G(x,S).  Depending  on  the  values  of  the  limits,  this  function 
may  be  nearly  rectangular  (similar  to  the  near  field  diffraction  pattern  of 
a  slit)  or  nearly  a  sinc-funct ion  (similar  to  the  far  field  diffraction 
pattern  of  a  slit).  In  all  cases,  the  most  rapid  change  in  the  value  of  the 
Fresnel  integral  occurs  at  those  values  of  P  where  either  the  upper  or  the 
lower  limit  is  equal  to  zero.  By  using  (3),  we  find  that  the  Fresnel 
''window''  has  a  nearly  uniform  value  for 


IP  I  <  $  X/A( 1-2 s)  . 
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We  substitute  (2)  into  (6)  and  find  that  the  Fresnel  integral  has  a  nearly 
uniform  value  for 


Ifll  i  X/H.  (7) 

vbicb  is  the  same  as  tbe  region  occupied  by  the  central  lobe  of  the 
sinc(0H/X)  weighting  due  to  tbe  finite  transducer  height.  The  Fresnel 
integral,  then,  has  the  effect  of  suppressing  the  sidelobes  of  the 
sine(0H/X)  function  as  well  as  producing  some  low  modulation  fringes  in  the 
central  lobe.  At  small  raluea  of  z,  the  Fresnel  integral  transitions  into  a 
broad  Fraunhofer  diffraction  pattern  of  a  narrow  slit  and  has  wary  little 
impaot  on  the  amplitude  of  G(x,0) . 

If  we  denote  the  total  amplitude  weighting  of  G(x,0)  by  f(x,0),  we  hawe 

that 

!  G(x,0)  *»  f (x,0)exp[jnxA(l-2s)0*/X*] .  (8) 

Ve  now  concentrate  on  t^  phase  part  of  G(x,0) .  The  key  point  is  that  the 
phase  factor  is  quadratic  in  0  and  linear  in  x.  At  the  transducer,  where 
!  *  ■  0,  there  is  no  optical  power,  whereas  the  optical  power  is  greatest  at 

x  "  L„  If  the  optical  power  could  be  canceled,  the  inverse  Fourier 
transform  would  confine  the  light  to  reetangular  channel  with  no  spreading 
into  adjacent  channels.  In  principle,  the  optical  po%  -  can  be  compensated 
I  by  a  section  of  a  oonical  leas  whose  power  varies  linearly  from  zero,  at 

x  ■  0,  to  a  value  of  LA(l-2s) /XFl  at  x  *  L.  An  alternative  method  is  to 
construct  a  holographic  correcting  element  from  one  channel  of  the  Bragg 
cell;  such  an  element  will  then  correct  the  acoustic  spreading  for  all 
channels  simultaneously  in  the  same  fashion  as  would  a  matched  filter. 

The  use  of  a  holographic  element  for  compensating  the  acoustic 
spreading  has  been  reported  by  Vodovatov,  e£  at* .  The  derivation  given 
there  for  the  value  of  G(x,0)  does  not  agree  with  (8),  particularly  with 
respect  to  the  phase  factor  having  quadratic  dependence  on  0  and  linear 
dependence  on  x.  The  result  from  (8),  however,  is  consistent  with  the 
observation  that  the  channels  can  be  separated  by  focusing  at  a  different 


plane  using  auxiliary  optica.  Aa  an  aaide,  it  it  easy  to  show  that  the 

axial  distance  from  the  Bragg  cell  at  which  the  channels  becone  separated  is 

equal  to  x(l-2s)A/X~  Thus,  the  "focal  line"  has  a  slope  of  (l-2s)A/X  with 

respect  to  the  horizontal  axis  ;  for  typical  Bragg  cell  paraaeters,  this 

o  o 

slope  corresponds  to  an  angle  of  88  to  89  which  is  nearly  parallel  to  the 
optical  axis. 

Te  now  consider  sone  details  relating  to  the  construction  of  a 
holographic  correcting  eleaent  based  on  the  result  given  by  (8) .  This 
analysis  provides  soae  useful  insights  into  the  perforaanee  relationships  of 
nul t i-channe 1  Bragg  cells.  Ve  construct  the  eleaent  by  using  a  line 
reference  source,  as  shown  in  Figure  2,  that  lies  horizontally  in  the  plane 
of  the  Bragg  cell  and  is  displaced  vertically  by  a  distance  D.  The  line 
source  is  produced  by  a  cylindrical  lens  that  collects  part  of  the  saae 
eolliaated  beaa  that  illuainates  the  Bragg  cell.  A  spherical-cylindrical 
lens  ooabination  then  creates  a  Fourier  transfora  in  the  vertical  direction 
and  an  iaage  in  the  horizontal  direction.  The  figure  also  shows  a  second 
Fourier/ imaging  lens  ooabination  that  will  be  used  later  to  create  a 
corrected  iaage  of  the  Bragg  cell  at  the  output  plane.  The  reference  beaa 
is  aodulated  by  a  device  (not  shown)  so  that  it  has  the  saae  temporal 
frequency  as  the  light  diffracted  by  the  Bragg  cell.  The  total  light 
distribution  in  the  Fourier  plane  is  then 

A(x.p)  -  laJ2*PD/X  +  G(x.p),  (9) 

where  K  is  the  aa^litude  of  the  reference  beaa.  This  light  distribution  is 
square— law  recorded  on  a  photographic  plate  and  the  developed  plate  is 
replaced  in  the  systea.  If  we  select  the  appropriate  diffracted  order  at 
the  output  plane,  the  effective  holographic  trsnsaitttance  function  is 

H(x,0)  ■  KG*(x,0)exp[j2nf)D/Xl  .  (10) 

The  action  of  the  holographic  eleaent  is  to  correct,  or  conjugate,  the 
cylindrical  phase  factor  contained  in  G(x,0) .  Te  see  froa  (8)  that  if  the 
holographic  eleaent  is  constructed  at  an  acoustic  wavelength  A  ,  there  will 
be  soae  residual  aberrations  as  the  drive  frequency  changes  to  produce  a 


different  acoustic  wave length  A .  Ve  now  examine  the  nature  of  this  error  as 
well  as  the  effects  of  positioning  errors  of  the  hologram. 

To  simplify  the  notation,  let  a  *  r(1-2s)A/1*  and  we  ignore  any 

amplitude  weighting  terms  as  well  as  the  linear  exponential  phase  factor 

that  merely  causes  an  overall  vertical  displacement  of  the  corrected  image. 

Suppose  that  the  holographic  filter  was  made  at  an  acoustic  wavelength  Aq. 

If  the  operational  wavelength  is  at  some  other  value,  aberrations  will  arise 

whose  phase  is  proportional  to  exp[ja0*x  -  j(a+ao)0*x],  where  sq  is  due  to 

the  value  of  A-  A  .  We  can  also  determine  the  effects  of  displacement  of 
o 

the  filter  in  the  0  and  x  directions  so  that  the  residual  phase  error  6 (x,0) 
becomes 

9(x,0)  -  a0*x  -  (a+ao)(0+0o)*(x+xo) ,  (11) 


where  0q  and  xq  are  the  positional  errors.  Note  that  the  principal  terms  in 
ap*x  cancel  so  that,  if  we  arrange  the  terms  in  descending  powers  of  0,  we 
have  that 

9 (*'P>  “  *0P*  +  *0‘xo  +  *0***0 

+  2(a+ao)0o0x  +  'l<a+ao>po0xo  (U> 

+  (a+ao)0Jx  +  <a+ao)0Jxo. 

The  first  three  terms  of  (12)  are  residual  errors  in  correcting  the 

cylindrical  phase  factor  and  lead  to  a  do focusing  of  the  output;  the  next 

two  terms  are  linear  in  0,  leading  to  a  slight  vertical  displacement  of  the 

output  ;  the  last  two  terms  are  not  functions  of  0  but  one  of  them  is  linear 

in  x  which  can  produce  a  horizontal  shift  in  the  output.  These  terms 

provide  a  convenient  diagnostic  tool  for  properly  aligning  the  filter. 

Suppose  that  after  the  filter  is  placed  in  the  system,  we  illuminate  it  with 

the  channel  driven  at  the  same  constructional  acoustic  wavelength.  In  this 

case,  a  >0  and  the  residual  terms  are 
o 


(13) 


+  2*fiofix  +  2»0o0xo 
+  *  «3»o*o. 

Since  the  effects  of  the  first  ten  (defocus ing)  may  be  difficult  to  detect, 

ve  first  observe  the  point  in  the  output  for  which  z  ■  0.  The  two  key  tens 

reaaining  are  2sB  Bx  end  ap*  x  .  The  first  of  these  tvo  tens  will  cause 
0  o  0  0 

the  output  to  shift  up  or  down  depending  on  the  signs  of  p^  and  x^.  Set 
■  0  by  noting  when  the  output  ia  syaetricslly  positioned  about  the  image 
of  the  reference  beaa;  the  reaaining  tens  are  2sPqPx  and  aP*ox.  If  ve  nov 
focus  at  x  -  L  to  get  aaxiaua  sensitivity,  ve  eaa  set  a  PQ  ■  0  by  again 
causing  ayaaetry  in  the  vertical  direotion  about  the  reference  beaa. 

Since  the  filter  is  nov  properly  positioned,  all  the  error  tens  froa 
(13)  are  equal  to  zero  so  that  ve  eaa  nov  exaaine  the  effects  of  a 
wavelength  change  (ao  #  0) .  The  only  reaaining  ten  is  then 
B(x,p)  “  exp[jaQ9*x],  a  defoeusing  ten  vhieh  shovs  the  degree  to  vhich  the 
aeoustio  spreading  has  not  been  ooapensated.  At  x  *  0,  ve  see  that  this 
ten  has  no  effeot  on  the  output,  whereas  at  x  *  L  the  phase  is  bounded  by 
±aoP*L.  This  suggests  that  ve  night  be  able  to  balance  the  aberration  so 
that  it  has  equal  but  opposite  values  at  x  ■  0  and  x  *  L.  This  can  be  done 
by  focusing  the  reference  beaa  at  that  plane,  axially  displaced  froa  the 
Bragg  cell,  vhere  the  self-focusing  property  of  the  diffracted  light  causes 
the  channels  to  be  fully  separated.  If  ve  let  xx  be  the  corresponding 
horizontal  distance  at  vhich  the  tvo  beaas  are  jointly  focused,  the  residual 
phase,  for  PQ  “  x0  "  ®»  1*  given  by 

9(x,P)  ■  -aP*xx  +  »0P*X*  (14) 

This  procedure,  in  effect,  provides  s  focusing  bias  ten  vhich  offsets  the 
bias  in  x  (i.e.,  x  is  bounded  by  the  range  0  to  L) .  To  balance  the 
aberrations  st  the  ends  of  the  cell,  ve  find  that  x.  *  (a  /a)L/2,  so  that 


(13) 


«(x,p>  -  +*0P»(x-L/2) 

and  tli«  Aberration*  ar*  proparly  balanced. 

T*  now  oonaider  the  optima  rain*  of  the  aoonatio  wavelength  to  be  naed 

in  constructing  the  holographic  eleaent.  A  Bragg  cell  is  generally 

specified  by  its  center  frequence  f  and  a  bandwidth  Af.  If  f,  ■  f  -  Af/2 

o  c 

and  fs  «  f^  +  Af/2,  the  Bragg  cell  is  operated  so  that  tt  1  2fx  to  prevent 
interaodnlation  products.  Since  A.  is  inversely  proportional  to  f, 
constructing  the  holographic  filter  at  the  center  frequency  does  not  equally 
divide  the  wavelength  rang*.  Instead,  w*  construct  the  holograa  so  that 

Ao  -  V( ft+f x)/2fafx  (1«) 

This  wavelength  choice  shows  that  th*  hologrsa  tends  to  favor  correction  of 
the  aberrations  at  the  lower  drive  frequencies  where  th*  acoustic  spreading 
is  greatest. 

V*  now  coatpar*  the  corrected  wavefront  aberration  with  the  uncorreoted 

wavefront.  Without  th*  holographic  eleaent,  we  have  a  maxima  phase  error 

equal  to  (a  'fa  )p*L,  whereas  th*  aaxiaua  error  with  correction  is  +i  p*L. 

©  o 

The  ratio  of  |a  |/(a+a  )  is  equal  to  (l-f,/f.)/2 j  for  the  worst  case  of  an 
o  o 

octave  bandwidth,  th*  ratio  is  1/4.  Therefore,  the  aaxiaita  phase  error  with 
correction  is  at  least  four  tiaes  less  than  th*  worst  case  (at  fi)  without 
correction,  and  at  least  two  tiaes  better  than  th*  best  ease  (at  fa)  without 
oorrotion.  If  th*  holographic  element  is  constructed  to  talar***  the 
aberrations  along  the  cell,  the  corrected  wavefront  is  yet  another  factor  of 
two  better  than  the  uncorreoted  case. 

3.  EXPERIMENTS 

V*  constructed  a  holographic  eleaent  froa  one  channel  of  a  32-channel 
Bragg  cell.  This  cell  has  a  nominal  bandwidth  Af  ■  80  MHz  at  a  center 
frequency  f  “  133  MHz.  The  cell  is  fabricated  froa  teluriua  dioxide 
notarial  having  a  longitudinal  velocity  of  4.2  Ea/sec.  The  active  length  of 
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tli*  crystal  is  L  *  6  an  so  that  the  tine  delay  is  "l.Spsee  and  the  tine- 
bandwidth  product  is  *120.  The  transducers  have  heights  equal  to  100  p  and 
they  are  placed  on  250  ft  center  spacings.  The  acoustic  wavelength  at  the 
oenter  frequency  is  31  |i.  This  Bragg  cell  was  originally  designed  and 
successfully  used  as  a  high-speed,  aulti-ohannel  light  nodulator  wherein 
only  the  region  near  the  transducer  was  used.*  Since  this  cell  produces  a 
considerable  aaount  of  aeoustie  spreading,  it  representa  a  severe  test  of 
the  ability  of  a  holographic  clsaent  to  channelise  the  light  over  a 
significant  aaount  of  the  horizontal  aperture. 

The  light  source  is  a  10  aw  He-Ne  laser  for  which  X  ”  <32.8  na.  The 

holographic  eleaenta  were  constructed  on  SO- 120  glass  plates  with  a 

reference-to-signal  beaa  ratio  of  unity  near  0  -  0.  After  the  plate  was 

developed  and  replaced  in  the  systea,  we  positioned  it  by  using  the 

procedure  outlined  above  to  set  0  and  z  equal  to  sero.  This  was  achieved 

o  o 

by  driving  only  the  channel  used  to  construct  the  holograa  (see  Figure  2), 
while  observing  the  corrected  inage  of  the  Bragg  cell  at  the  output  plane. 

Figure  3  shows  the  uneorrected  and  the  corrected  output  when  channels 

are  driven  at  the  frequency  corresponding  to  A0.  Figure  3a  shows  the  degree 

of  acoustic  spreading  for  the  uneorrected  oase.  Of  the  twelve  channels 

available,  we  drove  the  top  channel  to  illustrate  how  rapidly  the  acoustic 

energy  spreads  and  the  bottoa  four  channels  to  illustrate  how  the  besas 

overlap  and  add  coherently.  Te  see  that  the  beans  are  well  separsted  only 

in  the  region  near  the  transducer  as  uoted  earlier.  If  we  wish  to  process  a 

significant  aaount  of  tine  history,  we  find  that  the  beaas  begin  to  overlap 

at  a  position  that  is  only  10%  of  the  available  aperture.  Figure  3b  shows 

the  corrected  output;  we  see  that  there  is  very  little  evidence  of  beaa 

spreading  and  that  the  light  has  been  confined  to  within  ohannel  heights  ss 

deterained  by  the  transducers.  The  aaount  of  beaa  spreading  is  s  function 

of  the  drive  frequency;  it  is  greatest  at  the  low  frequency  band  edge  (due 

to  a  large  A)  and  least  at  the  high  frequency  band  edge  (where  A  is  snail) . 

Furtheraore,  as  shown  by  (12),  the  output  is  fully  corrected  only  when 

a  ■  0  «  z  ■  0.  If  0^  ■  z^  ■  0,  the  residual  aberration  is 
o  o  o  o  o 

6(z,0)  "  +sq0*z,  so  that  if  we  drive  the  Bragg  cell  at  its  lowest  and  highest 

frequencies,  we  can  see  the  effects  of  the  residual  aberrations. 


I 


Figures  4a  aad  4b  show  the  corrected  oatpot  at  the  low  aad  high  drive 

frequencies.  Ve  aee  ao  appreciable  ehaage  ia  the  tpreadiag  due  to  the 

reeidaal  aberratioa  aa  a  function  of  drive  frequency  except  at  the  lowest 

frequency.  This  was  coafiraed  visually  by  sweeping  the  drive  oscillator 

over  the  frequency  band  aad  aotiag  very  little  change  in  the  channel  height 

which  is  most  evideat  at  x  ■  L.  The  reasoa  oaa  be  seea  from  the  following 

argwaents.  The  value  of  4  is  given  by  (2)  ao  that  at  x  •  L,  the  acoustic 

o 

bean  extends  a  distance  h  above  aad  below  the  center  line  of  the  channel: 

h  -  A (l-2s)L/2H  (17) 
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For  the  paraaeters  of  the  TeOx  Bragg  cell  (a  -0.2).  we  have  that  h  **  558|i 
so  that  the  aain  beaa  spreads  through  the  adjacent  channel  and  nearly  to  the 
center  of  the  third  'channel  away.  This  behavior  is  consistent  with  the 
results  shown  in  Figure  3a.  kt  the  lowest  frequency,  the  spreading  at  the 
end  of  the  cell  increases  to  792|i,  whereas  at  the  highest  frequency  the 
spreading  decreases  to  432|t.  The  holographie  eleaent  corrects  the  spreading 
coavletely  at  the  aidband  frequency;  at  the  high  sad  low  frequencies  the 
spreading  is  that  produced  by  aa  equivalent  acoustic  wavelength 
A#  -  V[(fx-fx)/2fxfx] V  When  this  value  of  wavelength  is  used  in  (17)  the 
worst  case  spreading  is  reduced  to  182p.  This  degree  of  spreading  extends 
approx iaately  to  the  center  line  between  transducers  and  is  appsrent  only  st 
the  end  of  the  Bragg  cell  reacts  from  the  transducer.  If  the  correction 
were  aade  so  that  the  spreading  is  balanced  at  both  x  ■  0  and  x  -  L.  the 
spreading  would  be  further  reduced  by  a  factor  of  two. 

4.  THE  MODULATION  TRANSFER  FUNCTION 

Te  see.  then,  that  the  holographic  eleaent  confines  the  diffrscted 
light  waves  to  channels  as  defined  by  the  transducer  height.  The  next  issue 
is  to  deteraine  the  effect  that  this  correction  has  on  the  nodulstion 
transfer  function  of  the  Bragg  cell.  Ve  can  estiaate  the  degree  of 
aodulation  transfer  loss  ss  a  function  of  the  propagation  distance  by 
considering  the  Poynting  vector  surface  as  s  function  of  x.  The  acoustic 
energy  propagating  at  an  angle  $  with  respect  to  the  x-axis  (see  Figure  5) 
travels  at  a  velocity 


V 


(18) 


Suppose  that  we  were  to  drive  the  Bragg  cell  with  as  iapslae  function  and 
obaerve  the  aurface  defining  the  acouatic  energy  at  some  time  t.  The 
acoustic  energy  will  occupy  a  finite  pulse  width,  centered  on  the  curved 
line  shown  in  Figure  5,  consistent  with  the  reciprocal  bandwidth  of  the 
Bragg  cell.  The  energy  along  the  axis  will  have  traveled  a  distance  x,-Vt 
whereas  the  energy  at  an  angle  $  will  have  traveled  a  distance  .  The 
horizontal  projection  of  this  distance  is  x,  *  V^t  c os 4> : 

x3  =  Vt[l  -  *»/2<l-2s),  (19) 

and  the  difference  between  x4  and  x(  is 

dx  -  Vt<)»*/2(l-2a) .  (20) 

If  we  use  (2)  in  (20) ,  we  have 


dx  *  xA* (l-2s)/2H* . 


(21) 


The  effect  of  imaging  the  Bragg  cell  in  the  horizontal  direction  while 
correcting  for  spreading  in  the  vertical  direction  is  to  integrate  the 
energy  function  in  the  vertical  direction.  Ve  now  need  to  form  a  basis  for 
a  reasonable  estimate  for  the  resulting  pulse  width. 


The  3-dB  bandwidth  of  the  cell  is  Af  so  that,  if  there  were  no  acoustic 

spreading,  the  equivalent  rectangular  pulse  width8 would  be  Tq  •  1/Af  =  dQ/V. 

This  result  is  obtained  by  representing  the  MTF  by  sinc[d  (f-f  )/V]  where  f 

o  c  c 

is  the  center  frequency.  The  effect  of  integration  in  the  vertical 
direction  due  to  compensating  the  acoustic  spreading  is  approximately  dx/30, 
obtained  by  numerical  integration.  The  equivalent  rectangular  pulse  width 
d,  after  correction  is  then  dQ  +  dx/30 


d 


2 


V 

AF 


xA2  (l-2s) 
60H7 


+ 


(22) 


The  argonaut  of  tha  sinc-function  representing  the  corrected  MTF  is  then 

Z  *  d.(f-f  )/  V  or 
*  c 


(23) 


where  N  *  AfT  m  T/T  is  the  tine  bandwidth  product  of  the  Bragg  cell  and 

o 

x'  *  x/L  is  a  normalized  distance.  For  x'  *  0,  there  is  no  loss  in  MIT;  at 
x*  ■  1,  there  nay  be  some  loss  in  the  MTF.  For  the  parameters  of  the  Bragg 
cell  used  in  these  experiments,  the  MTF  is  reduced  by  0.8  dB  at  x*ml  and  by 
0.4  dB  at  x  '*1/2. 


Ve  now  calculate  the  number  of  channels  that  can  be  used,  subject  to 

the  restriction  that  the  main  lobes  of  the  far-field  acoustic  waves  do  not 

overlap.  If  the  Bragg  cell  crystal  is  square,  we  have  that  L  ■  (N  -l)kfl  +  H 

c 

N  kfl,  where  N  is  the  number  of  channels  and  kH  is  the  separation  between 
c  c 

transducers.  The  condition  for  avoiding  overlap  is  to  require  that  h,  as 
given  by  (17),  be  less  than  or  equal  to  kH/2  as  shown  in  Figure  S.  By 
combining  these  relationships,  we  have  that 

N  iH/\(l-2s),  (24) 

c 


where  A1  »  V/ft  is  the  worst  case  wavelength.  For  the  corrected  case,  we 
substitute  th  equivalent  value  of  the  acoustic  wavelength  which  is 

A  -  VAf/2flf1,  (25) 

© 

so  that 

Nc  i  Vflf  (l-2s)  ‘  <“> 

The  increase  in  the  number  of  usable  channels  is  then  A^/ A  =  2fa/Af, 

provided  that  Nc  is  not  greater  thsn  L/H.  Since  the  minimum  value  of  the 

transducer  height  may  be  set  by  other  considerations,  the  value  of  Nq  as 

given  by  (26) may  not  be  achievable.  For  the  parameters  of  the  Bragg  cell  we 

used,  the  maximum  value  of  N  without  correction  is  about  4;  with 

c 


correction,  N  is  sbont  17.  To  achieve  the  full  utility  of  the  devioe,  it 
0 

is  necessary  to  correct  the  spreading  so  that  (15)  is  satisfied  ;  in  that 
oaae  Nq  is  equal  to  34.  Equation  (26)  then  indicates  how  closely  packed  the 
ohannels  can  be  without  the  residual  aberrations  causing  the  light  to 
interact,  while  (23)  indicates  the  loss  in  aodulation  transfer 
funetion.Experinental  results  oonfira  the  analysis  that  the  acoustic 
spreading  can  be  successfully  corrected  without  introducing  an  excessive 
loss  in  the  MTF.  Ve  nodulated  the  carrier  with  a  12.5  KHz  square  wave  and 
aeasured  the  increase  in  the  equivalent  pulse  width  at  the  aidpoint  and  at 
the  two  ends  of  the  cell.  The  results  were  consistent  with  those  predicted 
by  the  MTF  as  given  by  (23) . 

5.  Cn«m1tn<Bit« 

Ve  have  shown  that  the  effects  of  acoustic  spreading  within  the  Bragg 
cell  can  be  coapensated  in  the  sense  that  the  diffracted  light  can  be  aore 
nearly  contained  within  channel a  aa  defined  by  the  transducers.  The  results 
show  that  the  beat  perforaanee,  both  in  terns  of  correcting  for  the  acoustic 
spreading  and  liniting  the  loas  in  MXF,  ia  obtained  when  we  use  a  Bragg  cell 
interaction  aediua  that  haa  a  large  value  of  s.  A  value  of  s  approaching 
1/2  iaplies  that  the  acoustic  spreading  is  intrinsically  snail;  such 
aaterials  aay  not.  however,  have  the  oharacteristica  required  to  achieve 
large  bandwidtha,  high  diffraction  effieieneiea  or  low  attenuation  per  unit 
length.  For  a  given  interaction  nediua,  the  results  show  that  the 
perforaance  inproves  as  the  fractional  bandwidth  decreases.  The  correction 
for  acoustic  spreading  can  also  be  iaproved  by  a  factor  off  two  if  the 
correction  ia  optiaized  at  the  aidpoint  off  the  cell  this  optiaization  does 
not,  in  general,  affect  the  MTF. 

The  degree  of  correction  required  is  application  dependent.  If  each 
channel  is  driven  by  an  independent  wideband  signal  plua  noise  and  we  wish 
to  correlate  these  signals  with  a  reference  signal  in  a  tine-integrating 
architecture,  a  high  degree  of  correction  is  required.  The  reason  is  that 
we  do  not  want  the  inforaation  in  any  one  channel  to  spread  into  adjacent 
channels;  there  are  other  applications  where  the  Bragg  cell  is  used  ss  a 
two-diaensional  light  aodulator  that  also  require  a  high  degree  of 
correction.  In  applications  such  as  processing  bistatic  or  nonostatic  radar 


return*  using  s  phased  array  antenna,  the  data  in  adjacent  channels  doe  to 
adjacent  antenna  eleaents  generally  does  not  change  rapidly.  Nevertheless, 
in  high  performance  radars,  where  the  return  pulse  must  be  processed  through 
correlation,  the  performance . may  be  improved  by  optically  correcting  for  the 
acoustic  spreading. 

In  both  of  the  general  applications  cited  above,  a  correlation 
operation  is  performed  on  the  received  time  signal  in  a  space-plane 
implementation.  An  alternative  way  to  perform  correlation  is  through 
frequency  plane  multiplication  techniques  wherein  we  create  the  two- 
dimensional  Fourier  transform  of  the  multichannel  Bragg  cell  and  multiply  it 
by  the  Fourier  transform  of  the  reference  signal.  If,  for  example,  the 
reference  function  is  time-invariant,  the  constructed  matched  filter  will 
automatically  incorporate  the  holographic  correcting  element  that  we  have 
implemented.  In  many  cases,  however,  the  reference  funetion  must  be 
programmable  (i.e.,  it  ia  time-variant)  so  that  the  multiplication  technique 
is  less  attractive  when  applied  to  multichannel  Bragg  cells. 

The  notion  of  frequency  plane  processing  also  suggests  applications 
such  as  emitter  sorting.  In  this  case,  the  two  dimensional  Fourier 
transform  of  the  multichannel  Bragg  cell  produces  an  optical  mapping  of  the 
microwave  field.  The  frequency  of  a  particular  emitter  is  displayed  along  a 
frequency  axis  while  the  azimuth  (or  direction  of  arrival)  is  displayed 
along  an  azimuth  axis.  Results  obtained  without  the  use  of  the  holographic 
correcting  element  show  that  the  diffraction  pattern  is  not  an  exactly 
orthogonal  sine- function;  in  the  vertical  direction  the  sine-functio;i 
follows  the  locus  of  a  parabola  having  its  foci  at  the  optical  axis.4  If 
the  holographic  correction  element  is  used,  the  parabola  becomes  a  straight 
line.  Sine*  the  azimuth  information  does  not  map  directly  into  the 
rectangular  coordinates  of  the  Fourier  plane  (i.e...  the  vertical 
displacement  is  proportional  to  sinij*),  some  post  processing  of  a 
photodetector  array  output  is  needed  to  properly  remap  the  azimuth 
information.  This  remapping  could  also  include  any  geometric  scaling  needed 
to  compensate  for  the  effects  of  acoustic  spreading. 


The  use  of  a  holographic  corrector  plate  for  the  acoustic  spreading,  as 
veil  aa  a  broader  understanding  of  the  basic  difffraction  patterns,  opens 
the  way  for  new  application  of  nultiehannel  Bragg  cells.  As  the  bandwidths 
of  conunication  and  collection  systesis  increase,  the  need  for  such  devices 
will  expand  because  they  help  utilize  the  full  parallel  processing 
capabilities  of  optical  systesis. 
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FIGURE  CAPTION  LIST 


FIGURE  1.  Model  for  Acoaatic  Spreadiag  la  a  Bragg  Cell. 

FIGURE  2.  Optical  Syatea  for  Coaatraetiag  the  Holographic  Eleaent  aad 

Correotiag  the  Acoaatic  Spreadiag. 

FIGURE  3.  Maltiehaaael  Bragg  Cell  Diffractioa  Beaas:  (a)  Uncorrected 

Case  Showing  Beaa  Overlap  ;  (b)  Corrected  Beaaa  at  Optiana 
Acoaatic  Wavelength  (123  MHz) . 

FIGURE  4.  Corrected  Beaaa  Showing  Eaergy  Confineaeat:  (a)  Reaalt  for 

175  MHz  ;  (b)  Reaalt  for  95  MHz. 

FIGURE  5.  Diagraa  Showing:,  the  Sorface  of  Acoaatic  Energy  aad  the 

Traaadacer  Geoaetry. 
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Adaptive  optical  processor 


A.  Vander  Lugt 


The  transversal  filter  is  widely  used  in  digital  data  processing.  In  this  paper,  we  describe  a  method  for  using 
optically  tapped,  Bragg  cell  delay  lines  to  implement  a  general  adaptive  linear  prediction  algorithm.  Flexi¬ 
bility  is  achieved  through  changing  the  electronic  signals.  The  implementation  is  in  the  Fourier  domain  so 
that  the  wide  bandwidth  of  the  Bragg  cells  can  be  fully  used.  A  high  dynamic  range  can  be  achieved  because 
the  system  is  interferometric  and  linear  in  light  amplitude. 


I.  Introduction 

In  this  paper  we  describe  an  adaptive  optical  pro¬ 
cessing  technique  wherein  the  processing  operation  is 
performed  in  the  frequency  domain.  Adaptive  filtering 
as  implemented  by  linear  prediction  has  been  applied 
in  several  areas.  A  tutorial  review  of  linear  prediction 
has  been  given  by  Makhoul1  and  a  survey  of  its  appli¬ 
cation  to  communications,  including  adaptive  tech¬ 
niques,  was  given  by  Lucky.2  The  early  work  by  Lucky 
was  on  adaptive  redundancy  removal  in  data3  and 
adaptive  equalization  of  digital  communications  sys¬ 
tems4'5  to  increase  the  rate  of  data  transmission. 
Variations  of  these  adaptive  schemes,  in  which  the 
system  is  equalized  either  by  transmitting  a  training  set 
or  by  operating  on  the  data  itself  (decision  directed), 
have  been  reported.6-8  Adaptive  equalization  signifi¬ 
cantly  increases  the  allowable  data  rate  over  a  filtered 
channel  for  a  given  bit  error  rate  required  at  the  output 
of  a  receiver. 

Adaptive  filtering  has  also  been  applied  to  antenna 
sidelobe  weighting  by  Widrow  et  al.9  and  others.1011 
Widrow  et  al 12  have  also  described  adaptive  noise- 
cancellation  techniques  such  as  notch  filtering  for  in¬ 
terference  rejection,  adaptive  self-tuning  for  spectral 
line  enhancement,  and  spectral  analysis.  Griffiths,  in 
particular,  has  addressed  the  problem  of  estimating  the 
instantaneous  frequency  contained  in  digital  signals.13 
The  same  basic  algorithms  are  finding  widespread  ap¬ 
plication  in  speech  processing14-16  and  in  combating 
intersymbol  interference.4-817 
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The  wide  range  of  applications  of  adaptive  filtering 
is  due  partly  to  the  flexibility  of  the  transversal  filter 
that  is  a  basic  element  of  all  these  schemes.  A  trans¬ 
versal  filter  consists  of  a  tapped  delay  line  that  contains 
the  discrete  time-sampled  values  of  a  received  signal. 
The  outputs  of  each  tap  are  weighted  as  determined  by 
the  processing  operation  and  summed  to  provide  an 
estimate  of  the  signal.  If  the  weights  are  selected  to 
provide  the  best  estimate  of  the  received  signal,  the 
system  is  called  a  linear  predictor  or  estimator.  If  the 
estimated  signal  is  subtracted  from  the  received  signal 
and  if  the  residual  signal  is  used  to  control  the  tap 
weights,  the  system  is  an  adaptive  linear  predictor.  We 
would  like  to  extend  these  processing  algorithms  to 
handle  wide  bandwidth  signals. 

An  acoustooptic  cell  has  wide  bandwidth  and  behaves 
as  a  delay  line  which  can  be  tapped  optically;  we  de¬ 
scribe  here  an  optical  processing  architecture  that  can 
be  used  to  implement  operations  similar  to  those  cited. 
In  Sec.  II,  we  briefly  review  the  basic  theory  of  adaptive 
linear  prediction  and  derive  the  equivalent  frequency 
plane  operator.  We  then  describe  the  adaptive  optical 
processing  technique  and  relate  its  operation  to  those 
that  have  been  implemented  digitally.  An  important 
observation  is  that  the  optica)  system  can  operate  with 
either  analog  (continuous-time)  or  sampled  (discrete- 
time)  signals  with  no  change  in  the  basic  components 
in  the  system. 

II.  Adaptive  Linear  Prediction 

We  now  describe  the  basic  function  of  an  adaptive 
linear  predictor.  A  general  version  of  the  predictor  for 
discrete-time  signals  is  shown  in  Fig.  1.  Let  s(t )  be  the 
received  signal  that  has  been  sampled  at  the  proper  rate. 
The  delayed  samples  are  multiplied  by  the  weights  c;, 
j  =  l  to  N,  to  produce  an  estimate  s(t): 

N 

S(f  )  *  ^  CyJiU  —  }D\,  (  1 ) 

where  D  is  the  time  delay  of  each  element  of  the  tapped 
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Fig.  1.  Block  diagram  of  adaptive  linear  predictor 
using  discrete  tapped  delay  lines. 


delay  line,  and  N  is  the  number  of  delays.  The  opti¬ 
mum  weights  must  be  determined  based  on  our  criterion 
for  how  well  the  estimated  signal  s(t)  represents  s(t). 
A  commonly  used  criterioi-  is  to  minimize  the  squared 
error  of  the  residual  signal  z  [t)  =  s(t)  —  s(t). 

Let  E  be  the  total  energy  of  the  residual  signal: 

E*IWD-I  cjs(t  -  jD)|  .  (2) 

M  l  j-l  I 

where  the  energy  is  summed  over  M  samples.  We 
minimize  E  with  respect  to  a  specific  tap  weight  c, 
through  the  derivative 

SE  * 

—  *■  -2  Z  s(t)  -  Z  cjtU  -  jD)  sU  -  iD).  (3) 

<9c,  m  y-i 

If  we  set  the  derivative  equal  to  zero,  we  find  that  the 
optimum  fixed  tap  weights  are  given  by  the  solutions 
to  the  N  equations 

Z  c,R((»  -  j)D )  •  R(i£>),  i  «  1  toN,  (4) 

i- 1 

in  which  R(-)  can  be  thought  of  as  the  covariance  matrix 
if  the  sum  is  over  a  set  of  sample  values,  or  as  the  cor¬ 
relation  function  if  the  sum  is  over  time. 

In  some  applications,  the  characteristics  of  s(t )  are 
such  that  we  can  improve  the  performance  of  the  system 
by  adapting  the  weights  in  time.  If  the  rate  of  change 
is  slow  (i.e.,  the  process  is  quasi-stationary),  we  can  use 
the  method  of  steepest  descent3-9  to  determine  the  di¬ 
rection  in  which  we  need  to  adjust  the  weights.  From 
Eq.  (3)  we  see  that  the  value  of  the  derivative  can  also 
be  computed  as 

dE 

—  «-2Z  i(t)s(t-jD),  jm  ItolV.  (5) 

dCj  M 

We  can  calculate  these  derivatives  if  we  multiply  the 
residual  signal  by  the  delayed  values  of  v,he  received 
signal  and  use  this  signal  to  update  the  tap  weights. 
This  process  is  shown  in  Fig.  1,  wherein  the  predicted 
signal  i(t)  is  subtracted  from  s(t),  the  residual  z(t)  is 
multiplied  by  a  loop  gain  G,  and  the  residual  signal  is 
then  multiplied  by  the  delayed  received  signal  at  the 
output  of  each  tap  to  provide  the  derivative  for  adjust¬ 
ing  the  tap  weight  as  given  by  Eq.  (5*. 

Rhodes18  described  an  optical  implementation  of  the 
adaptive  linear  prediction  filter  shown  in  Fig.  1.  The 


implementation  was  in  the  time  domain  using  a  com¬ 
bination  of  an  electrooptic  light  modulator  and  an  ac¬ 
oustooptic  Bragg  cell  in  a  time-integrating  architecture 
to  form  the  product  given  in  Eq.  (5).  The  integration 
and  the  storage  of  the  tap  weights  was  performed  by  a 
liquid  crystal  spatial  modulator;  the  spatial  positions 
in  the  light  modulator  correspond  to  the  time  delays.  A 
second  Bragg  cell  was  then  used  to  perform  the  opera¬ 
tion  given  by  Eq.  (1)  and  a  photodetector  provided  the 
signal  s(t).  This  implementation  has  some  disadvan¬ 
tages  in  that  the  full  bandwidth  of  the  Bragg  cells  is  not 
used  and  the  spatial  modulator  has  limited  dynamic 
range.  The  implementation  can  also  be  performed  in 
the  frequency  domain  to  more  fully  utilize  the  available 
bandwidth.  The  main  purpose  of  this  paper  is  to  de¬ 
velop  such  a  technique;  its  advantages  will  become  ap¬ 
parent  later. 

To  cast  the  problem  into  frequency  domain  notation, 
we  must  first  describe  the  adaptive  linear  predictor  in 
terms  of  continuous- time  variables.  In  the  optical 
implementation,  the  delay  line  taps  are  not  discrete, 
although  the  number  of  independent  tap  weights  is 
determined  by  the  time-bandwidth  product  of  the 
Bragg  cells.  We  begin  by  calculating  the  value  of  the 
tap  weights.  Let  the  ;  th  weight  at  time  t  be  the  sum  of 
its  value  at  the  previous  update  time  t  —  tx  and  the 
value  of  the  error  derivative: 

Cj(t )  «  c/U  -  t- i)  -  G  f*  ztu)slu  -  jD)du,  (6) 

Jt-T\ 

where  T\  is  the  integration  time  of  the  accumulator 
whose  output  provides  the  updated  weight  value. 
When  theaystem  reaches  steady  state,  we  find  that  c;  (f ) 
=  Cj(t  —  T\),  so  that  the  value  of  the  integral  in  Eq.  (6) 
must  be  zero.  Since  the  integral  is  the  cross-correlation 
of  the  received  signal  and  the  residual  signal,  we  have 
the  physical  insight  that  the  system  adapts  its  weights 
until  the  residual  contains  no  components  that  are 
correlated  with  the  input;  this  is  sometimes  called  the 
correlation  cancellation  loop  method  for  implementing 
adaptive  filtering.14  This  method  for  calculating  the 
adaptive  weights  can  be  done  in  real  time  and  saves  the 
computational  complexity  of  inverting  the  matrix  for¬ 
mulation  given  by  Eq.  1 4). 

The  passage  to  a  continuous-time  representation  is 
fairly  straightforward.  We  represent  the  delayed  signal 
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as  s(t  -  t );  it  is  multiplied  by  c(r)  and  integrated  over 
the  total  time  delay  T  of  the  delay  line  to  give 

sU(  “  (If)dt  -  r|dr  (7) 

The  value  of  the  continuous-weight  function  c(r)  can 
be  found  by  repeated  application  of  Eq  (6)  or  directly 
by  noting  that 

c(t)  » G  C  ziutsiu  -  r)du.  (8) 

Jt-Tt 

If  the  accumulator  time  window  T i  is  sufficiently  long 
relative  to  the  variations  in  the  statistics  of  the  input 
signal,  the  weights  are  not  a  function  of  the  present  time 
t.  During  the  adaptation  period,  or  during  periods 
when  the  signal  statistics  change  significantly,  the 
weights  are  functions  of  both  t  and  T.  We  retain  the 
implicit  relationship  of  c(t)  to  time  t  as  in  Eq.  (8)  to  aid 
in  deriving  the  Fourier  domain  solution.  We  now 
substitute  Eq.  (8)  into  Eq.  (7)  to  get 

s(i)  »  G  C  f  z(u)s(u  -  r)s(t  —  r)dudr.  (9) 
«/0 

After  some  changes  of  variables,  we  can  express  Eq. 
(9)  as 

s(t)  -  G  x;x  zit  +  q)s(t  +  q  +  r)s(t  +  T)dqdr.  (10) 

To  assist  in  solving  Eq.  (10),  we  need  a  definition  for  the 
instantaneous  Fourier  transform  of  a  truncated  con¬ 
tinuous-time  signal.  We  define 

Friui.t)  *•  P  fiu)  exp(—juu)du,  (11) 

Jt-T 

or,  in  an  equivalent  and  sometimes  useful  form,  we 
have 


expOwf  )Fr(u,r)  «* 


rect 


fit  +  u)  exp<~juu)du. 


(12) 


We  use  the  rectangular  function  and  the  inversion  of  Eq. 
(12)  in  Eq.  (10)  to  get  that 


Ht)mG  f‘.SdqdT[S-.  Zr((<»),t)  exp[/ui(t  +  <7)|<fai 

X  sit  +  q  +  r)  |^*  SriZ.t)  exp[/{(t  +  r)|d£ 


(13) 


or  that 

i(t>  -  G  xx:  Zr,l<«i,t)Sr({,f)  exp[/(w  +  $)t] 


*  xx:  sit  +  q  +  r)  expl/fuK;  +  l-T)\dqdTd£dui.  (14) 

We  now  perform  the  integration  on  r,  using  a  change  of 
variables  in  which  f  +  q  +  r  =  u.  Further,  although  the 
range  of  integration  on  r  appears  to  be  infinite,  the 
rectangular  functions  must  be  used  again  to  limit  the 
range  on  u  to  the  interval  t  +  q  -  T  tot  +  q.  The  in¬ 
tegral  on  u  now  becomes 

*\u  >  exp (j(u  )du  exp(-;£(f  +  <7)] 

t  +  q-  T 

“  Sri-£.t  +  <7 >  exp[-;£(t  +  c?)|  (15) 


by  virtue  of  Eq.  (11).  If  the  statistics  of  s(t)  are  slowly 
varying  relative  to  the  total  delay  time  represented  by 
the  taps,  we  have  Sr(—^,t  +  q)  a  Sr(_?,f  )•  Further¬ 
more,  we  note  from  Eq.  (11)  that  Sr(-£,f)  =  SV<{,f)  if 
the  received  signal  is  real.  We  can  now  combine  Eq. 
(15)  with  Eq.  (14)  to  get 


s(t)  "  G  xx:  Zt,  (<*>,! )  |  St(  u>,t )  ( 2  expOut ) 

X  |  J*  exp(-;'<?({  -  o))|d(?j  d£du.  (16) 


The  integral  on  q  is  bounded  by  the  range  0  <  q  <  T\\ 
we  see  that  this  integral  has  the  form  sinc[(£  —  w)TJ2\. 
We  note,  however,  that  the  definition  of  )  is  such 
that  it  already  contains  the  effects  of  convolution  with 
a  sine  function  due  to  a  time  truncation.  If  T\  >  T,  the 
additional  convolution  in  Eq.  (16)  will  not  affect  the 
results  and  we  can  complete  the  integrations  to  get  the 
final  result: 


S(t)  =  G  ^  Zr,(w.D|Sr(<*>,t)|2 exp(/‘«t)dw.  (17) 

If  the  received  signal  is  bandlimited  to  a  bandwidth  W, 
the  limits  in  Eq.  (17)  are  finite  and  restricted  to  the 
frequency  range  of  the  signal. 

Equation  (17)  is  the  frequency  domain  representation 
of  the  estimated  signal  s(t).  We  carried  out  the  deri¬ 
vation  in  some  detail,  because  treating  the  limits  of  in¬ 
tegration  required  considerable  care  throughout. 
Furthermore,  the  derivation  gives  some  insight  into  how 
the  result  given  by  Eq.  (17)  can  be  implemented  by  the 
use  of  Bragg  cells. 

III.  Adaptive  Optical  Processor 

If  a  Bragg  cell  is  positioned  so  that  its  transducer  is 
located  on  the  optical  axis,  we  can  represent  a  signal 
traveling  in  the  negative  direction  as  f(t  +  x/u),  where 
u  is  the  acoustic  velocity,  and  x  is  the  spatial  coordinate. 
Equation  (10)  suggests  that  we  need  to  convolve  three 
such  functions  whose  values  are  shown  in  Fig.  2.  The 
function  z{t  +  q)  can  be  produced  by  a  Bragg  cell  whose 
transducer  is  located  at  the  q  -  0  axis,  with  the  wave 
propagating  in  the  negative  q  direction.  The  function 
s(f  +  r)  can  be  produced  by  a  Bragg  cell  whose  trans¬ 
ducer  is  located  at  the  r  =  0  axis,  with  the  wave  propa¬ 
gating  in  the  negative  r  direction.  The  function  s(t  + 
q  +  t)  must  be  multiplied  by  the  other  two  functions. 
We  see  that  if  a  third  Bragg  cell  has  its  transducer  lo¬ 
cated  on  the  q  =  — r  axis,  with  the  wave  propagating  in 
the  negative  q  +  r  direction,  we  have  formed  the  ap¬ 
propriate  products.  The  limits  of  integration  are  also 
shown  in  Fig.  2. 

We  now  show  how  to  implement  these  functions  op¬ 
tically  by  first  generating  the  product  equivalent  to  z(t 
+  q)s(t  +  t).  Consider  the  optical  configuration  shown 
in  Fig.  3.  This  subsystem  contains  two  orthogonally 
oriented  Bragg  cells  in  which  the  information  contained 
in  the  first  Bragg  cell  is  multiplied  by  the  information 
contained  in  the  second.  These  two  Bragg  cells  are  then 
imaged  onto  a  common  plane  P.\,  although  in  practice 
it  is  sufficient  that  the  optical  system  perform  the  ap- 
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Fig.  2.  Convolution  of  three  signals  necessary  to  produce  estimate 
of  received  signal. 


■use  cm 


Fig.  3.  Optical  implementation  to  produce  a  spatial  cross-correlation 
function. 

propriate  Fourier  transform  of  the  product  of  the  signals 
contained  in  the  Bragg  cells.  Each  Bragg  cell  has  a  time 
delay  T  and  a  length  L  -  vT.  The  bandwidth  of  the 
Bragg  cells  is  W,  so  that  the  number  of  independent 
time  delays  is  N  =  2 TW.  Both  Bragg  cells  are  posi¬ 
tioned  so  that  their  transducers  are  located  on  the  op¬ 
tical  axis;  this  is  done  partly  because  it  simpl  ifies  the 
mathematics  and  partly  because  it  leads  to  some  in¬ 
teresting  physical  interpretations  of  the  results. 

The  first  Bragg  cell  in  plane  Pi  is  driven  by  the  re¬ 
ceived  signal  s(t ),  which  has  been  translated  to  a  suit¬ 
able  center  frequency  <*>c.  The  undiffracted  light  is 
removed,  and  the  Bragg  cell  is  imaged  onto  plane  P3  so 
that  the  light  distribution  is  uniform  in  the  y  direction. 
In  the  x  direction  we  have 

<n(x,t)  *■  jmwfaFTL  s(t  +  x/u)  exp[/<i)c(t  +  x/u) j,  (18) 

where  mj  is  the  modulation  index  of  the  first  Bragg  cell, 
P  is  the  laser  power,  and  a  is  that  fraction  of  the  laser 
power  that  reaches  plane  P\.  This  function  is  now 
multiplied  by  the  transmittance  of  the  second  Bragg 
cell,  which  is 


Fig.  4.  Optical  system  to  produce  cross-spectrum  of  truncated  time 
signals. 


a2{y,t)  =  jmim-2  \' ctPa'i/L  r(t  +  y/v )  exp\ju>r(t  +  v /(.’)] ,  (19) 

where  m2,  the  modulation  index  of  the  second  Bragg 
cell,  is  generally  different  from  mu  and  a;  is  the  average 
signal  power  in  the  time  interval  T. 

We  now  arrange  to  take  the  2-D  Fourier  transform  of 
a1(x,t)a2(y,t)  as  shown  in  Fig.  4: 

/»0  /+0 

A(p,q,t )  =  J  J  ai(x,t)a2(.v,t)  exp[-;(px  +  qy)]dxdv.  (20) 

We  want  to  show  that  A(p,q,t)  is  equivalent  to  the 
Fourier  transform  of  the  cross-correlation  of  z(t)  and 
s(t).  As  shown  by  Eq.  (20),  the  integral  is  separable  in 
x  and  y.  We  need  to  couple  the  time-delay  arguments 
of  a i(jc ,t )  and  a2(y,t)  so  that  the  cross-correlation  is 
displayed  as  a  function  of  space  instead  of  time.  One 
way  to  do  this  is  to  rotate  the  coordinate  axes  in  plane 
P3  followed  by  a  1-D  Fourier  transform.19  An  alter¬ 
native  method  is  to  perform  the  2-D  Fourier  transform 
and  to  evaluate  it  along  the  line p  -  q.  This  is  equiva¬ 
lent  to  a  rotation  of  the  axes  in  the  Fourier  plane  and 
using  a  narrow  slit  oriented  along  the  q'  =  0  axis,  where 
p'  and  q’  are  the  new  coordinates.  We  then  have  p  = 
q  -  p'/\/2,  which  we  use  in  Eq.  (20)  along  with  the  as¬ 
sociations  that  rx  =  x/u,  Ty  =  y/u,  and  u;  =  up'/\/2- 
Equation  (20)  can  now  be  written  as 

X0  ^*0 

r  J  ^s(f  +  Tx)Z(f  +  Ty)  eXpl/uMf  +  Tx  ,  ] 

Xexp[/uic(f-  rv)|  exp[—  /a)(rt  +  ry)|orxdrv,  (21) 

where  A\  =  —  m t m 2 v  aPj;/L.  We  solve  the  integral 
on  tz  first.  We  have 

XO 

s(t  +  rx)  exoL/u)cU  +  ''J)  exp(-ju)Tx)dTx. 

Let  t  +  tx  =  u,  so  that 

A\(w,t)  »  f  s(u)  expi;u'ru)  exp[-;w(u  -  f)ldu.  (22) 
Jt-T 

Given  the  definition  of  the  instantaneous  Fourier 
transform  of  a  time  function  as  in  Eq.  ( 11 ),  we  recognize 
that  Eq.  (22)  can  be  expressed  as 

/ti(ui.t)  *  -  “Jr.t)  exppuJt).  (23) 

In  a  similar  fashion,  we  evaluate  the  integral  on  r, 
and  combine  all  factors  to  get 
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Fig.  5.  Second  leg  of  interferometric  system. 


A(tai.t)  =  v2AiSt^ui  -  uv.tiZrioi  -  wc,£)  exp{j2ut).  (24) 

This  result  shows  that  the  complex  valued  light  am¬ 
plitude  along  the  line  p  =  q  in  plane  P4  is  the  product 
of  the  Fourier  transforms  of  s(t)  and  z(£).  By  the  use 
of  the  convolution  theorem,  A(co,t)  is  also  the  Fourier 
transform  of  the  cross-correlation  of  s(t)  and  u(t)  as 
desired. 

We  now  want  to  multiply  .4(oi,£ )  by  a  third  function 
associated  with  s(t  +  q  +  r )  and  integrate  the  product 
over  the  frequency  plane.  This  can  be  achieved  inter - 
ferometrically  by  using  a  third  Bragg  cell  in  a  separate 
path  of  the  interferometer.  Consider  the  optical  system 
shown  in  Fig.  5,  which  shows  the  second  part  of  the  in¬ 
terferometer.  Since  the  important  part  of  A(p,q,t)  lies 
along  the  line  p  -  q  and  is  centered  at  p'  =  \f2w/v,  we 
orient  the  third  Bragg  cell  at  a  45°  angle  to  the  x  —  y 
coordinate  system  of  plane  P5.  To  preserve  the  fre¬ 
quency  scaling,  the  Fourier  transform  lens  L4  must  have 
a  focal  length  a  factor  of  \r2  longer  than  that  of  L3.  In 
this  path,  we  use  a  cylindrical  lens  to  image  the  line  il¬ 
lumination  of  the  Bragg  cell  at  plane  P4  to  conserve  the 
optical  power. 

The  transducer  of  the  Bragg  cell  is  located  at  the 
optical  axis  and  is  driven  by  s(t).  If  (1  —  a)P  is  the 
fraction  of  the  laser  power  available  at  plane  P5,  the 
effective  amplitude  of  the  signal  at  plane  P5  is 

bU'.t)  *  jm3  v  (1  -  a)P/L  sU  +  x'/v)  expjjwc(t  +  x'/u)\,  (25) 

where  x’  is  the  coordinate  axis  in  plane  P5  that  contains 
the  Bragg  cell.  The  1-D  Fourier  transform  is  now 

J-»o 

sit  +  tx)  expl/ujcU  +  rt)|  exp(-;<i>rx)dTt, 


where _BX  =  jm 3  Vd  -  a)P/L,  tx  =  x'h  and  u>  = 
up'/VT  By  following  the  same  procedure  as  before,  we 
find  that 

Biu.t)  =  cBtSriui  -  uv  ,t )  6xpi  jut  1.  (27* 


Suppose  that  we  now  place  a  single  photodetector  at 
plane  P4  along  the  line  p  =  q,  located  so  that  its  center 
is  at  and  its  physical  extent  corresponds  to  the  range 
|w  —  u>c|  <  W/2.  The  photodetector  current  will  be 
proportional  to  the  intensity  at  plane  P4: 

X(*»c+  W/2 

Uui,t)d  a),  (28) 

,c—  IV/2 

where  C  is  a  constant  that  includes  the  photodetector 
sensitivity  S,  expressed  in  amps/watt  as  well  as  the 
conversion  factors  necessary  to  cause  Eq.  (28)  to  be  di¬ 
mensionally  correct.  We  have 

/  (u>,£ )  =  |.4(u),()  +  B(u;.t)|2 

-  |.4(o*,t)|2  +  |B(ai,£ )|2  +  2Re)/l I u\t )B* ( u.t )].  (29) 

The  photodetector  current  can,  therefore,  be  separated 
into  three  terms  as  given  by  Eqs.  (28)  and  (29).  From 
Eqs.  (24)  and  (27),  we  note  that  the  first  two  terms  are 
not  functions  of  time;  they  therefore  contribute  only  to 
a  bias  current.  This  statement  is  not  strictly  true,  be¬ 
cause  5r(w,f )  and  Zr(uo,t)  are  functions  of  time  through 
the  definition  of  the  instantaneous  Fourier  transform. 
These  functions  art;  slowly  varying,  however,  relative 
to  the  variations  in  the  third  term  and  can  be  removed 
by  a  bandpass  filter  after  detection. 

The  main  term  of  interest  is  the  third  term  of  Eq.  (29) 
which,  when  substituted  into  Eq.  (28)  with  the  results 
of  Eqs.  (24)  and  (27),  leads  to 

J-.  U..+  W/2 

Zr.iw  “  u;c,f )|.<i27'( w  -  u>c.t)|2 

W/2 

x  expO  wOdu*.  (30) 

By  a  change  of  variables,  we  can  write  Eq.  (30)  as 

XW/2 

Zr,(a'.()|.Sr(w.t  )|2 

■W/2 

x  expOuit  *dui.  |31) 

We  see  that  the  integral  in  Eq.  131 )  is  exactly  the  one  we 
wanted  to  evaluate  [compare  Eq.  <  17)].  The  output 
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Fig.  6.  Arrangement  of  Bragg  ceil  to  generate  delayed  estimate. 


occupies  bandwidth  W  and  is  centered  at  ujc.  By 
comparison,  the  first  two  terms  of  Eq.  (28)  have  a 
bandwidth  of  <1  IT  centered  at  w  =  0.  The  third  term 
can,  therefore,  be  easily  separated  from  the  low-fre¬ 
quency  terms.  It  turns  out  that  Blw,t)  is  in  phase 
quadrature  with  A(cj,C)  as  can  be  seen  by  considering 
the  values  of  A  i  and  B\.  The  correct  phase  relationship 
can  be  obtained  by  introducing  a  tt/2  phase  shift  in  the 
carrier  used  to  drive  the  third  Bragg  cell  or  by  optically 
phase  shifting  the  two  beams  in  the  interferometer. 

The  optical  architecture  described  so  far  is  charac¬ 
teristic  of  one  in  which  we  form  an  estimate  s(t )  at  the 
present  time  t  from  some  past  history  of  s(t).  We  may 
also  wish  to  introduce  an  additional  time  delay  a  to 
arrive  at  an  estimate  i(t  —  a)  or  to  predict  the  received 
signal  at  some  time  t  +  a.20  Suppose  that  we  geomet¬ 
rically  move  the  third  Bragg  cell  by  an  amount  x  =  ua 
as  shown  in  Fig.  6.  The  signal  within  the  cell  is  now  s  ( t 
-  a  +  tx)  exp[/a)c(f  -  a  +  rx )],  and  if  we  perform  the 
Fourier  transform  as  before,  we  arrive  at  the  result 

B(o).f)  »  vB\StIw  -  uic,t)  exp[/<*)(f  -  a)],  (32) 

which,  in  turn,  leads  to 

Xwn 

wn 

X  exp[/a>(f  -  a)\d(i>.  (33) 

The  output  of  the  photodetector  then  becomes  the  es¬ 
timate  s(t  —  a).  If  the  Bragg  cell  is  moved  in  the  op¬ 
posite  direction,  the  sign  of  a  is  reversed  and  we  have 
the  prediction  §{t  +  a).  The  interesting  interpretation 
of  this  result  is  that  the  system  processes  the  received 
signal  in  such  a  way  that  the  estimate  (or  prediction) 
occurs  at  the  time  value  corresponding  to  the  physical 
position  x  =  0.  Thus,  if  the  Bragg  cell  .s  moved  so  that 
the  signal  s(t  -  «)  occurs  at  x  ■  0  (or  rr  =  0),  the  solu¬ 
tion  is  one  of  estimation.  If  the  Bragg  cell  is  moved  so 
that  the  entire  cell  is  located  in  the  negative  x  plane, 
there  is  no  received  signal  at  x  *  0,  ;ind  the  system 
predicts  a  signal  sit  +  a)  located  at  x  =  0. 

In  many  apolications  it  is  implicitly  assumed  that  the 
estimate  is  at  some  fairly  large  value  of  a  relative  to  T. 
For  example,  in  the  transversal  filtering  iiterature,  the 
center  tap  is  often  set  at  the  time  delay  77 2.  Some 
physical  phenomena  such  as  channel  distortion  or  in¬ 
tersymbol  interference  have  an  asymmetrical  response, 
so  that  better  results  can  be  achieved  by  selecting  the 
tap  position  to  weight  the  tails  of  the  response.  In 
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general,  then,  one  Bragg  cell  may  be  more  nearly  cen¬ 
tered  on  the  optical  axis;  the  mathematical  analysis  is 
thereby  basically  unaltered  except  for  a  residual  time- 
delay  factor  in  the  output. 

The  frequency  plane  notation  also  provides  a  more 
direct  insight  to  certain  filtering  operations.  Let  us 
consider,  as  an  example,  a  filtering  application  in  which 
we  want  to  remove  the  effects  of  a  strong  narrowband 
signal  that  interferes  with  the  information  signal.  F rom 
Fig.  1  we  see  that 

2(f)  »  sit)  -  sit),  (34) 

so  that,  in  the  sense  of  the  definition  of  the  instanta¬ 
neous  Fourier  transform,  we  also  have 


2t*(w)  *  Sr(u>)  —  5t(u>)  (35) 

From  Eq.  (17)  we  use  the  result  that  $7(0))  =  GZrua)- 
jS7-(w)|2,  so  that 


1  +  G|5t<uj)|2 


(36) 


It  is  tempting  to  define  the  transfer  function  for  the 
system  as  H(w)  =  ZrM/SrM.  In  general,  this  is  not 
a  valid  procedure  because  the  system  behaves  in  a 
nonlinear  fashion  during  adaptation.  Suppose,  how¬ 
ever,  that  the  system  has  reached  its  steady-state  con¬ 
dition  and  that  the  input  signal  statistics  are  not 
changing.  Then  we  can  say  that 


H(W)« - « - — — - —  •  (4.1 

At-Iui)  l+G|Sx(ui)|- 

We  see  that  the  overall  transfer  function  is  unity  where 
the  value  of  G  |  Sr  ( u; )  | 2  is  small.  The  presence  of  strong 
signal  components  around,  say,  u>0  produces  an  inverse 
filter  H(a>)  ~  1/G|St(w0)|2,  which  suppresses  the  en¬ 
ergy  near  u>„.  In  a  general  sense,  we  see  that  the 
equivalent  filter  adapts  to  whiten  the  input  signal 
spectrum. 

The  adaptation  time  is  dependent  on  the  specific 
application,  although  it  appears  that  mast  systems  cited 
in  the  literature  adapt  within  a  time  interval  corre¬ 
sponding  to  500-2000  data  samples.  The  adaptation 
time  is  determined,  in  part,  by  the  required  values  of  T 
and  T these  values  also  influence  noise  in  the  tap 
weights  and,  therefore,  the  accuracy  at  convergence. 
The  gain  factor  G  also  affects  both  the  rate  of  adapta¬ 
tion  and  the  accuracy.  Finally,  the  dynamic  range  of 
the  photodetector  determines  the  minimum  increment 


by  which  the  tap  weights  can  be  set.  Generally,  at  least 
10*  distinct  levels  provide  adequate  performance.  How 
these  various  factors  interrelate  is  highly  dependent  on 
the  specific  application. 


IV.  Summary  and  Conclusions 

The  adaptive  optical  processing  technique  described 
here  has  several  important  features.  First,  since  it  is 
the  optical  equivalent  of  an  adaptive  digital  transversal 
filter,  it  affords  a  high  degree  of  flexibility  and  is  ap¬ 
plicable  to  solving  a  wide  range  of  problems.  Second, 
implementing  the  process  in  the  Fourier  domain  allows 
the  use  of  the  full  wide  bandwidth  of  Bragg  cells.  We 
can  envision,  for  example,  the  processing  of  data  in  the 
100-1000-MHz  range.  Third,  the  processing  is  done  on 
the  amplitude  and  not  the  intensity  of  the  light,  so  that 
we  do  not  need  to  introduce  bias  terms  to  generate  both 
positive  and  negative  weights.  Because  the  system  is 
interferometric  in  the  temporal  sense,  the  dynamic 
range  is  high  and  the  system  is  insensitive  to  scattered 
light.  By  virtue  of  having  a  single  photodetector  and 
no  requirement  for  an  intermediate  storage  medium, 
both  phase  and  amplitude  errors  are  reduced. 


This  work  was  supported  by  U.S.  Army  Research 
Office. 
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In  my  paper  on  adaptive  optical  processing.1  1  was  not 
consistent  in  drawing  the  axes  for  the  figures  to  accurately 
represent  the  mathematical  functions.  If  we  replace  the 
variable  r  by  -r  and  change  the  limits  of  integration  in  Eq. 
( 10).  the  net  effect  is  to  reverse  the  orientation  of  the  hori¬ 
zontal  Bragg  cell  in  Fig.  3.  The  final  results  are  not  affected 
hv  this  change. 

Also,  in  evaluating  Eq.  1 14),  it  is  better  to  integrate  s(t  +  q 
—  r)  over  q  and  r  simultaneously  instead  of  sequentially. 
This  procedure  products  the  same  results  without  the  some¬ 
times  restrictive  assumption  that  +  q)  =■  Sr<-$,£  ). 

R«f«r*nc* 

1.  A.  VanderLugt.  Appl.  Opt.  21.  4005  1 1982). 


APPENDIX  E 


OPTICAL  TRANSVERSAL  PROCESSOR  FOR  NOTCH  FILTERING 
TO  APPEAR  IN  OPTICAL  ENGINEERING 
VOLOMME  23,  MAY -JUNE,  1984 


OPTICAL  TRANSVERSAL  PROCESSOR 


FOR  NOTCH  FILTERING 


A.  VanderLugt 

Advanced  Technology  Department 
Government  Systems  Sector 
Harris  Corporation 
Melbourne,  Florida  32901 


ABSTRACT 


A  frequency  domain  implementation  of  an  optical  transversal  processor 
has  been  described  previously.  Since  this  system  uses  Bragg  cells  both  as 
the  delay  line  and  as  the  accumulators  which  provide  the  tap  weights,  a  key 
question  is  what  effect  the  finite  integration  times  have  on  the  performance 
of  the  system.  Computer  programs  were  written  to  simulate  an  adaptive  notch 
filtering  application;  the  measure  of  performance  is  the  correlation 
coefficient  for  the  residual  signal  and  the  desired  received  signal.  The 
correlation  coefficient  was  increased  significantly  by  tapering  the 
accumulators  so  that  the  readaptation  phenomena  caused  by  large  values 
leaving  the  accumulator  are  minimized.  Several  examples  of  the  performance 
are  given  as  a  function  off  the  number  of  taps,  the  length  and  degree  of 
taper  of  the  accumulator,  the  feedback  gain  and  the  number  of  iterations. 

The  results  show  that  a  finite  accumulator  is  not  a  serious  drawback, 
particularly  for  those  applications  where  the  system  must  operate  in  a 
rapidly  changing  environment.  The  performance  of  the  system  then  approaches 
that  of  one  having  an  infinite  accumulator  with  the  gain  adjusted  to  give 
equivalent  tracking  performance. 


h 


OPTICAL  TRANSVERSAL  PROCESSOR 
FOR  NOTCH  FILTERING 


1 .  INTRODUCTION 

Adaptive  filtering,  using  transversal  tapped  delay  lines  with 
feedback,  has  been  applied  to  problems  such  as  redundancy  removal  in  data, 
reducing  intersymbol  interference  through  equalization,  noise  cancellation, 

1  _  g 

self-tuning  and  speech  processing.  The  data  rate  or  the  signal  bandwidth 

is  limited  by  the  processing  speeds  of  digital  circuits;  we  wish  to 
investigate  the  use  of  optical  processing  to  extend  the  bandwidths  of 
signals  that  can  be  processed. 

An  acousto-optic  cell  has  a  wide  bandwidth  and  behaves  as  a  delay  line 

9-10 

that  can  be  tapped  optically.  Rhodes  described  an  optical 

implementation  of  an  adaptive  linear  prediction  filter  in  the  time-domain, 

using  a  combination  of  an  electro-optic  light  modulator  and  an  acousto-optic 
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cell  in  a  time-integrating  architecture.  VanderLugt  described  an 

implementation  in  the  frequency-domain,  using  three  Bragg  cells  in  a  space- 
integrating,  interferometric  architecture.  The  main  purpose  of  this  paper 
is  to  give  the  results  of  some  computer  simulations  designed  to  study  the 
effects  of  using  an  accumulator  whose  integration  time  is  finite. 

2 .  BACKGROUND 

A  general  form  of  an  adaptive  linear  predictor  is  shown  in  Figure  1. 
Let  s(t)  be  a  sampled  signal  that  drives  a  delay  line.  The  delayed  samples 
are  multiplied  : j  the  weights  c  ^  to  form  an  estimate  s(t)  of  the  reci  /ed 
signal.  The  optimum  weights  are  determined  by  minimizing  the  mean  square 


error  in  the  residual  signal  z(t)  which  leads  to  the  relationship  that 


where  G  is  the  gain  in  the  feedback  loop,  T±  is  the  integration  time  of  the 
accumulator  and  x  represents  the  cont inuous-time  equivalent  of  the  discrete 


delay.  If  we  let  Tq  be  the  minimum  discernible  time  delay  in  a  Bragg  cell, 

then  NT  is  the  total  delay  time  T  of  the  cell.  The  estimate  of  the 
o  J 

received  signal  is 


s  ( t ) 


T 

c(T)s(t-T)dT. 

a 


(2) 


so  that,  by  substituting  (2)  into  (1),  we  have 

t  T 

§(t)  =  G  j  J  z (u) s (u-T) s (t-T)dudT. 
t-Tj.  ° 

By  a  change  of  variables  we  can  rewrite  (3)  as 


H(t) 


’ 

T 


z (t+q) s (t+q-T ) s (t-T ) dqdT 


(3) 


(4) 


which. 


in  turn,  can  be  used  to  get  the  frequency  domain  representation 
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■5  ( t ) 


oo 


|  Z  (v,t)  Is  (u),t) 
i  t  '  T 


:  j^t 
e  duj. 


(5) 


where  S^.(o),t)  and  Z^,(u.t)  are  the  instantaneous  Fourier  transforms  of  the 
most  recent  T  seconds  of  the  received  signal  and  the  residual  signal. 


There  are  several  ways  to  configure  an  optical  system,  using  three 


i 
1 1 


I 


Bragg  cells,  to  implement  (4)  or  (5);  one  of  these  is  given  in  Figure  2. 

The  functions  z(t)  and  s(t)  drive  two  Bragg  cells  in  the  directions  shown. 
These  two  Bragg  cells  are  mutually  imaged  onto  plane  P,,  with  the  diffracted 
light  from  the  first  Bragg  cell  passing  through  the  second  Bragg  cell. 

At  plane  P,  we  have  the  light  distribution  necessary  to  provide  the  tap 
weights.  This  can  be  more  clearly  seen  from  Figure  3a  in  which  the  tap 
weight  plane  has  been  rotated  by  45  degrees  for  convenience.  The  numbers  in 
each  delay  cell  represent  the  successive  time  indices  for  the  two  functions. 
If  we  integrate  the  light  in  the  vertical  direction,  as  shown  by  the  dotted 
lines,  we  obtain  the  function  c(r).  This  function  is  then  convolved  with 
the  received  signal  s(t)  which  drives  a  third  Bragg  cell  as  we  shall  show 
shortly . 

The  number  of  tap  weights  available  is  equal  to  twice  the  time- 
bandwidth  product  of  the  Bragg  cell  (Figure  3a  shows  only  a  small  number  of 
taps).  From  the  figure  we  see  that  the  central  tap  weight  is  the  sum  of  N 
products  whereas  the  end  tap  weights  contain  only  one  product.  The  end  tap 
weights  are  therefore  likely  to  be  noisy  and  reduce  the  performance  of  the 
system.  The  number  of  taps  required  depends  on  the  application.  To  reduce 
intersymbol  interference,  the  number  of  taps  required  is  determined  by  the 
extent  of  the  channel  impulse  response  which  may  be  of  the  order  o'  10-20 
symbol  periods.  For  adaptive  notch  filtering,  the  number  of  taps  required 
is  determined  by  the  number  of  frequencies  that  need  to  be  removed 
simultaneously'  we  require  at  least  two  taps  per  frequency.  For  the  notch 
filtering  simulations  reported  here,  approximately  50  taps  were  used  for 


most  of  the  tests  . 


If  the  time-bandwidth  product  of  the  Bragg  cells  is  of  the  order  of 
1000-2000,  the  use  of  only  50  taps  means  that  the  triangular  weighting  of 
c(t)  is  less  significant,  and  we  can  truncate  the  region  of  integration  to 
the  rectangular  window  as  shown  in  Figure  3b.  Each  of  the  50  taps  will  then 
contain  approximately  2000-4000  samples  which  represent  the  finite  length  of 
the  accumulators  in  each  of  the  correlation  cancellation  loops. 

Figure  4  shows  that  part  of  the  optical  system  which  simultaneously 
integrates  the  light  in  the  vertical  direction  and  creates  the  Fourier 
transform  of  the  tap  weights  in  the  horizontal  direction.  Along  the 
horizontal  axis  the  light  distribution  is  proportional  to  Z^.(a>,  t )  S^.(b>,  t )  . 

The  third  Bragg  cell  (not  shown  here),  is  arranged  so  that  S^.(u>,t)  is  added 
to  A(a>,t)  in  an  interferometer.  The  sum  is  then  square  law  detected  to 
provide  the  integrand  of  (5),  and  the  output  of  a  single  element 
photodetector  provides  the  estimate  s(t). 

The  advantage  of  implementing  the  adaptive  transversal  processor  in 
the  Fourier  domain  are  that  the  system  is  linear  in  amplitude  so  that  no 
bias  terms  are  needed.  As  a  result,  both  positive  and  negative  tap  weights 
are  generated,  and  the  dynamic  range  is  high.  The  full  bandwidth  of  the 
Bragg  cells  can  also  be  used.  A  potential  disadvantage  is  the  finite 
accumulators  that  may  reduce  the  performance  of  the  system.  Our  major 
interest  is,  therefore,  tc  simulate  the  operation  of  the  system  under  the 
constraint  of  using  finite  accumulators. 

3.  THE  SIMULATIONS 


We  wrote  a  computer  program  to  calculate  the  performance  of  the 
adaptive  system.  The  application  we  chose  to  study  is  that  of  notch 
filtering  in  which  the  received  signal  consists  of  a  wideband  signal  g(t) 
and  one  or  more  cosine  jammers: 


(6) 


s(t)  =  g  (t)  +  Z-,  A^cosok  t. 

1  =  1 

The  wideband  signal  is  either  a  very  long  psendo  noise  sequence  or  a 

Gaussian  random  signal;  each  signal  has  zero  mean  and  unity  variance  a  *  . 

M  8 

The  variance  of  the  cosine  jammer  is  a1  =  A*/2.  The  input  signal-to- 

c  jTi  J 

jammer  ratio  is  then  SNR.  =  10  log(o  */o  *)  =  -10  log(c  *). 

i  g  c  c 

Figure  5  shows  the  spectrum  |S(w)l  of  the  input  signal  when 

SNR^  *  -18.5  dB  for  a  system  having  SO  taps  and  a  fairly  short  500  sample 

accumulator.  The  spectrum  was  obtained  by  letting  the  system  iterate  for 

2000  input  data  samples  and  then  calculating  a  1024  point  discrete  Fourier 

transform  of  the  most  recent  data  samples.  The  magnitude  squared  of  the 

system  transfer  function  |H(«i>)  I  is  also  shown;  it  is  a  notch  that  has 

formed  adaptively  to  cancel  the  cosine  jammer  at  a  frequency  u>o .  If  we  let 

Tq  be  the  time  delay  between  taps,  the  frequency  of  the  jammer  shown  is 

u  *  2n/9T  .  The  spectrum  |Z(u) I  of  the  residual  signal  is  also  shown, 
o  o 

multiplied  by  a  factor  of  0.05  for  clarity.  We  see  that  the  input  and 


output  spectra  are  similar  except  in  the  region  uQ  where  we  see  some 
evidence  of  the  cosine  jammer  in  the  residual  signal. 


The  presence  of  some  jammer  energy  in  the  residual  points  out  a 
fundamental  difference  between  the  frequency  plane  implementation  of 
adaptive  processing  and  a  digital  system  where  the  accumulator  can  be  made 
arbitrarily  long.  In  a  digital  system,  the  tap  weights  converge,  after  some 
number  of  iterations,  to  stable  values  so  that  extresiely  deep  nulls  can  be 


formed.  In  the  initial  stages  of  adaptation,  the  tap  weights  change  toward 


their  final  values  at  a  rate  determined  primarily  by  the  feedback  gain  and 
the  jammer  amplitude.  In  the  latter  stages  of  adaptation,  the  new 
contributions  to  the  tap  weights  are  small  in  magnitude  and  tend  to  alternte 
in  sign.  After  convergence,  we  find  that  s(t)  is  a  good  estimate  of  the 
amplitude  and  phase  of  the  cosine  jammer  so  that  the  null  is  very  deep;  the 
null  remains  dfep  until  there  is  some  change  in  the  statistics  of  the 
received  signal. 

In  the  optical  implementation,  the  tap  weights  rapidly  assume  values 
close  to  their  final  values  and  an  increasingly  deeper  notch  is  formed  at 
the  jammer  frequency.  After  the  number  of  iterations  corresponding  to  the 
accumulator  length,  however,  the  earliest  values  in  the  accumulator  (which 
are  large  because  they  contribute  most  to  initial  adaptation)  begin  to  leave 
the  accumulator.  The  cosine  component  in  the  residual  then  increases 
momentarily  because  s(t)  is  not  as  good  an  estimate  of  the  jammer  signal. 

The  system  then  goes  into  a  "readaptation"  phase  where  the  new  contributions 
to  the  tap  weights  must  offset  the  older  contributions  that  are  leaving  the 
accumulator.  This  process  repeats  every  MTq  seconds,  where  M  is  the  number 
of  samples  in  the  accumulator,  but  with  successively  less  impact  on  the 
residual  signal.  After  enough  iterations,  the  system  converges  to  a  stsady 
state  where  the  jaauaer  signal  level  in  the  residuul  is  just  enough  to 
maintain  the  tap  weights  at  a  (nearly)  fixed  value  over  time. 

The  adaptation  and  readptation  can  be  seen  from  Figure  6,  in  which  the 
upper  trace  is  that  of  z(t)  for  the  first  2000  data  samples.  We  note  that 
adaptation  is  relatively  rapid;  after  about  50  data  samples  the  residual 
signal  is  primarily  composed  of  the  wideband  signal  g(t).  At  the  end  of  the 
500^  data  sample,  the  earliest  tap  weights  are  leaving  the  accumulator  and 
the  system  "readapts"  to  the  new  conditions.  Similar  leadapt a t  ions  occur  at 


multiples  of  500  data  samples  and  the  system  performance  is  not  optimum.  In 
an  effort  to  reduce  the  effects  of  readaptation,  we  used  a  taper  on  the 
accumulator  so  that  the  values  in  the  accumulator  are  attenuated 
exponentially  as  they  flow  through.  When  we  used  an  exponential  taper  for 
which  the  oldest  accumulator  values  are  reduced  by  a  factor  of  one-half 
relative  to  their  initial  values,  the  residual  signal  took  on  the  values 
shown  in  the  lower  trace.  There  is  now  little  evidence  of  readaptation 
except  for  a  very  slight  perturbation  after  500  data  samples.  Note  that  in 
both  cases  z(t)  has  a  zero  mean;  we  have  added  and  subtracted  a  value  of  20 
to  separate  the  two  traces . 

The  basic  cause  of  the  readaptation  is  shown  in  Figure  7,  in  which  the 
upper  trace  represents  the  values  within  the  accumulator  associated  with  the 
first  tap  during  the  first  500  iterations,  using  no  taper.  As  described 
before,  the  earliest  contribution  to  the  tap  weights  are  large  (those  values 
from  450  to  500)  which  drive  the  first  tap  weight  to  nearly  its  final  value. 
Subsequent  contributions  are  small  and  oscillate  about  zero  (a  value  of 

0.025  was  added  to  the  value  to  separate  the  two  traces).  The  lower  trace 

shows  the  results  when  an  exponential  taper  whose  final  value  is  0.5  is 
used.  The  initial  contributions  to  the  accumulator  are  large  when  they 
enter  the  accumulator,  leading  to  rapid  adaptation.  After  they  have 
propagated  to  the  end  of  the  accumulator,  their  values  have  boen  reduced  by 
0.5  so  that  when  they  leave,  their  impact  on  s(t)  is  not  so  severe.  There 

is,  therefore,  less  evidence  of  readaptation  when  the  taper  is  used. 

An  exponential  taper  can  be  implemented  optically  by  placing  a  mask  at 
the  plane  where  z(t)  and  s(t)  are  mutually  imaged.  The  transmittance  of  the 
mask  would  be  unity  for  small  signal  delays  and  decay  exponentially  to  a 
final  vaue  of  50%  transmittance  at  the  end  of  the  accumulators.  In  the 


discussion  to  follow,  D  represents  the  final  value  for  the  exponential 
taper.  Other  tapering  functions  such  as  a  Gaussian  weighting  may  be  equally 
effective,  or  it  nay  be  possible  to  implement  the  taper  through  variable 
illumination  of  the  Bragg  cells  so  that  a  separate  mask  may  be  necessary. 

One  question  relates  to  whether  D  =  0.5  is  optimum.  Other  questions 
relate  to  the  optimum  gain  as  a  function  of  D,  the  impact  that  the  number  of 
taps  or  the  length  of  the  accumulator  have  on  the  performance  of  the  system, 
the  depth  of  the  notch  and  the  output  signal-to-noise  ratio.  To  perform 
these  tradeoffs  we  need  a  measure  of  performance.  One  measure  is  the  degree 
of  distortion  between  z(t)  and  g(t).  Since  the  applications  of  adaptive 
notch  filtering  often  require  a  subsequent  correlation  between  z(t)  and  g(t) 
over  some  long  time  interval,  a  second  measure  is  the  correlation 
coefficient  for  the  signals  z(t)  and  g(t).  As  we  now  show,  these  two 
measures  are  related. 

Consider  the  function 


y(t)  ■  z(t)  -  g(t) 


(7) 


having  mean  squared  error 


ay*  =  EKi  -  |)»} 


This  function  can  be  expanded  to  give 
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where  is  the  covariance  of  z(t)  and  g(t).  We  can  normalize  the 

covariance  to  get  the  correlation  coefficient: 
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If  the  residual  signal  is  identical  to  g(t),  the  distortion  is  zero  and  the 
correlation  coefficient  is  one.  We  shall  use  the  correlation  coefficient  p 
as  our  measure  of  system  performance. 

The  next  result  shows  the  sensitivity  of  the  performance  to  the 
exponential  taper  parameter  D.  Figure  8  shows  the  correlation  coefficient 
plotted  as  a  function  of  0  for  several  values  of  the  accumulator  length  M 
(expressed  as  the  number  of  samples  stored  in  the  accumulator) .  The  number 
of  taps  in  this  example  is  50.  and  we  calculated  p  from  (10)  over  the  last 
1024  points,  having  let  the  system  operate  for  2000  iterations.  We  note 
that  the  performance  improves,  as  expected,  as  M  increases  and  that  the 
performance  is  not  a  sensitive  function  of  D  except  for  D  >  0.8  or  D  <  0.1. 
If  0  >  0.8,  we  have  the  readaptation  phenomena  described  before.  For 
D  <  0.1,  the  exponential  taper  reduces  the  effective  length  of  the 
accumulator  so  that  a  longer  accumulator  with  a  high  degree  of  taper  may 
yield  a  performance  inferior  to  that  of  a  shorter  accumulator  with  less 
taper.  In  all  cases  shown,  the  optimum  value  of  D  is  about  0.5  and  this 
value  has  been  used  in  most  of  the  subsequent  tests.  It  should  be  noted 
that,  for  these  preliminary  tests,  both  the  accumultor  length  and  the  number 
of  iterations  are  fairly  small  to  conserve  computer  run  time.  The 
performance  of  the  system  improves  as  these  two  parameters  increase  as  we 
shall  show  later:  these  truncated  tests  are  designed  to  show  the  trends  in 


performance . 


Figure  9  shows  the  value  of  the  correlation  coefficient  as  a  function 


of  gain  with  the  number  of  taps  as  a  parameter.  In  this  example,  the 
accumulator  length  was  set  at  M  =  S00,  D  =  O.S  and  the  number  of  iterations 
was  2000.  We  would  expect  that  the  gain  should  be  inversely  proportional  to 
the  number  of  taps  because  the  tap  weights,  for  an  ideal  system,  should  be 
c(r)  =  (2/T)  cos  <aot.  The  data  from  Figure  9  supports  this  result.  The 
performance  improves  as  the  number  of  taps  increases,  but  it  should  be  noted 
that  my  simulations  do  not  account  for  the  triangular  area  of  integration 
naturally  imposed  by  the  Bragg  cells  as  shown  in  Figure  3a.  In  any  event, 
the  correlation  loss  for  the  case  of  N  *  50  versus  that  of  N  *  200  is 
-0.59  dB  versus  -0.5  dB  which  is  not  a  significant  change.  If  the  number  of 
taps  is  at  least  50,  then,  the  system  performance  is  quite  good  for  notch 
filtering. 

The  depth  of  the  notch  cannot  be  accurately  read  from  the  graphs  such 

as  that  given  in  Figure  5,  partly  because  the  period  of  the  jammer  is  not  a 

divisor  of  1024  which  is  the  sample  index  for  the  discrete  Fourier  transform 

used  to  calculate  the  steady  state  transfer  function.  The  notch  depth  can 

be  obtained  by  freezing  the  tap  weights  after  adaptation  and  letting  only 

the  jammer  signal  continue.  The  notch  depth  is  the  ratio  of  the  variance  <r* 

of  the  residual,  which  is  an  attenuated  replica  of  the  jammer,  to  the 

variance  a  *  of  the  input  jammer.  For  the  example  given  in  Figure  5,  the 
c 

notch  depth  by  this  method  of  calculation  is  -31.5  dB.  In  a  similar 
fashion,  the  output  s igna 1- to-no ise  rati  can  be  estimated  by  the  ratio  of 
a*  to  a*  after  the  weights  are  frozen.  Table  1  gives  the  notch  depth 
and  output  s  igna 1- to-no i se  ratio  for  a  single  jammer,  with  N  =  48. 
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TABLE  1.  Performance  Parameters 


►  • 


SNR. 


-19.9  dB  0.02S 


-13.9  dB  0.003 


-7.9  dB  0.007 


-1.9  dB  0.008 


SNR 


o  Notch  Depth 


16.0  dB  -35.9  dB 


25.6  dB  -39.5  dB 


21.3  dB  -29.2  dB 


20.8  dB  -22.7  dB 


U 


M  **  4000,  and  D  3  0.5.  The  gain  was  adjusted  to  produce  the  maximum 
correlation  coefficient  at  the  point  where  the  tap  weights  were  frozen  for 
the  strongest  jammer.  Since  the  tap  weights  were  frozen  after  only  dOOO 
iterations,  the  system  has  not  yet  reached  its  optimum  steady  state 
performance  and  the  values  given  in  Table  1  are  understated  somewhat.  Since 
the  residual  must  always  contain  a  portion  of  the  jammer  energy,  we  see  that 

the  notch  depth  is  a  function  of  the  input  jammer  power.  As  the  input 

jammer  level  decreases,  the  notch  depth  decreases  to  allow  the  feedback  loop 
to  maintain  the  proper  set  of  tap  weights. 

The  results  of  the  simulations  can  be  compared  with  those  obtained  by 
digital  systems  having  infinite  accumulators.  The  most  obvious  difference 
is  in  the  depth  of  the  notch.  Since  the  residual  must  contain  some 
component  of  the  jammer  signal,  the  notch  depth  cannot  exceed  a  given  value. 
In  passing,  it  is  worth  noting  that  if  a  separate  noise  signal  were 
available,  as  in  some  of  the  appliations  noted  by  Widrow,  e_t  al* ,  this 
restriction  is  removed  and  deeper  notches  could  be  obtained.  Let  us  assume, 
then,  that  the  component  of  the  jammer  in  the  residual  is  of  the  form 

B  cos[<oQt  +  0]  when  the  jammer  is  of  the  form  A  cos  i»0t.  If  the  jammer  and 

wideband  signal  g(t)  are  uncorrelated,  the  steady  state  transfer  function 
for  u  near  can  be  approximated  by  H(w) ,  where 

,  H  (u))  i  *  =  1-2R  COsKio-'j)  )  NT  (1  +  4/N)/2]  sine  ((  .  )  NT  /2  ] 

o  o  J 

(H) 

+•  R^Sinc2  (  (uj-w  JNT./21  , 

where 


R  —  *  kAAoNM 


(  12) 


and  where  G  is  the  feedback  gain,  N  is  the  number  of  taps  of  delay  T  .  M  is 
the  number  of  samples  stored  in  the  accumulator  and  k  is  a  factor  which 
accunts  for  the  taper  along  the  accumulator.  If  the  taper  along  the 
accumulator  is  severe,  the  ontput  of  the  accumulator  will  be  samller  (k  <  1) 
than  that  for  an  untapered  accumulator  (k  =  1).  To  achieve  the  same  notch 
depth,  the  feedback  gain  needs  to  be  increased  somewhat  to  offset  the 
effects  of  the  taper. 

If  we  use  the  value  of  B  in  (11)  and  calculate  the  notch  depth  at 

w  ■  w  ,  we  have  that 
o 


!h(u)q) 


kGA2NM  I2 
4+kGA2NMJ 


(13) 


We  note  that  the  notch  depth  increases  as  the  jammer  amplitude  A  increases 
to  the  point  where  the  value  of  the  jammer  component  B  in  the  residnal  have 
the  value  given  by  (12).  From  our  simulations,  however,  we  find  that  the 
length  of  the  accumulator  is  the  most  important  parameter  for  obtaining  a 
very  deep  null.  As  the  factor  kGAaNM  becomes  very  large  relative  to  the 
factor  4,  the  notch  depth  increases  and  we  have  that  B  ~  4/kGANM.  For  an 
exponential  taper  in  which  D  =  0.5,  we  find  that  k  =  0.49;  this  value,  along 
with  the  values  of  G,  A,  N  and  M  from  the  simulation  gives  a  calculated 
notch  depth  of  -35.4  dB  vs  a  depth  of  -35...  dB  as  given  by  the  data. 

The  half— power  bandwidth  can  also  be  obtained  from  (11).  If  the  value 
of  Af  is  small  compared  to  fQ,  we  can  expand  each  term  of  (11)  and  arrive  at 


the  result  that 


„  _  _ 2 _ 

'  TINT 

o 

Although  the  bandwidth  is  a  function  of  the  notch  depth  through  the 
parameter  R,  we  can  simplify  (14)  for  values  of  R  >  0.9  (i.e.  for  notch 
depths  more  than  20  dB  down)  so  that 


2  -  R  -  1/2R 
4/3  -  R/3  +4/N 


Af 


1 

NT  /1+4/N 

0 


(15) 


This  result  shows  that  the  half-power  bandwidth  is  primarily  determined  by 
the  number  of  taps  used.  If  N  is  large  compared  to  4,  the  bandwidth  is 
given  by  Af  =  1/2NTq.  The  bandwidths  as  determined  by  the  simulations  agree 
closely  with  the  values  calculated  from  (IS).  Figure  10  shows  the  transfer 
function  for  both  a  finite  (4000  samples)  and  infinite  accumulator  for  a 
N  »  48  tap  situation.  The  half-power  bandwidth  is  very  nearly  the  same  for 
the  two  cases,  but  the  notch  depth  is  about  6  dB  greater  for  the  infinite 
accumulator.  The  data  was  taken  after  10,000  iterations;  the  notch  depth 
for  the  infinite  accumulator  should  continue  to  grow  as  the  nmiber  of 
iterations  increase,  if  we  assume  that  the  jammer  and  the  signal  are 
uncorrelated. 

The  adaptation  time  constant  for  low  values  of  the  gain,  an  infinite 

6  7 

accumulator  and  uncorrelated  input  signal  and  noise,  is  given  by 


or 


r  =  -1/ ln( 1  -  Gk  ) 
P  P 


(16) 


(17) 


The  misadjustment  given  by  (18)  originates  from  random  noise  in  the  tap 
weights  due  to  the  use  of  a  noisy  gradient  instead  of  the  true  gradient  in 
the  adaptation  process.  The  relationships  in  (16-18)  are  valid  only  when 
the  adaptation  process  is  slow;  this  generally  implies  that  the  gain  is 
small.  When  the  gain  is  small,  a  second  type  of  Disadjustment  may  arise  due 
to  nons t a t j onary  input  statistics;  this  has  been  called  the  m i sadj us tment 

1  4 

due  to  lag.  As  noted  by  Widrow  ,  the  best  peiloimance  is  achieved  by 
setting  the  gain  so  that  the  Disadjustments  due  to  graident  noise  and  lag 


are  equa 1 . 


In  terms  of  the  notation  used  here,  the  total  mi  sad j us tment  can  be 
represented  by  Ij.  =  (o^1  ~  s*nce  j4®®61,  power  should  be 

completely  removed  in  the  optimum  Weiner  solution.  As  we  have  noted  before 
the  finite  length  of  the  accumulator  in  the  optical  implementation  requires 
that  some  jammer  signal  be  present  in  the  residual.  The  total  misadj ustment 
is,  therefore,  the  sum  of  that  due  to  gradient  noise,  lag,  and  jammer  power 
Further,  as  we  see  from  Figure  9,  the  performance  does  not  improve 
monotonically  as  the  gain  decreases,  since  the  accumulators  have  finite 
length.  As  a  result,  the  optimum  gain  is  generally  higher  than  that  used  i 
the  accumulator  were  infinite  so  that  the  misadjustment  is  principally  due 
to  gradient  noise  and  jammer  power. 

Since  the  gain  is  a  function  of  the  accumulator  length,  the  value  of 
the  adaptation  time  constant  as  calculated  from  (16)  or  (17)  tends  to  lead 
to  erroneous  results.  Abetter  procedure  is  to  calculate  the  jammer  power 
in  the  residual  from  (12)  and  to  then  calculate  the  misadjustment  due  to 
gradient  noise  from 


A  = 


B;/2  - 


a 


(19) 


The  adaptation  time  constant  can  then  be  calculated  from  (18).  As  an 

example  of  this  procedure,  we  use  the  results  from  the  simulations  given  in 

Table  1.  We  found,  at  the  end  of  8000  iterations,  that  a  1  -  1.12  and 

z 

_  i 

=  1.  The  gain  was  G  =  1.25(10  )  so  that,  from  (12),  we  find  that 

B  =  0.234.  From  (19),  we  find  that  the  misadjustment  due  to  gradient  noise 

is  then  0.0926  and,  from  (18),  that  the  average  time  constant  is  t  =  130 

°  ave 


We  now  give  an  illustration  of  the  adaptation  time  constant  and 
determine  how  rapidly  the  system  will  readapt  to  a  sudden  change  in  jammer 
frequency.  Figure  11  shows  the  values  of  s(t)  and  s(t)  in  the  region  of  an 
abrupt  change  in  frequency  from  <ao  =  2n/8TQ  to  u»o  =  2n/12To.  The  plot  of 
s(t)  is  delayed  by  one  sample  value  for  clarity.  The  system  was  allowed  to 
adapt  to  the  first  jammer  for  8000  iterations  so  that  s(t)  is  a  good 
estimate  of  the  jammer.  At  the  8001st/  iteration,  we  switched  to  the  lower 

A 

frequency,  and  we  see  that  s(t)  initially  decays  as  the  transient  passes 
through  the  system.  After  25  iterations  it  reaches  its  smallest  amplitude  ; 
it  then  begins  to  increase  in  amplitude  as  well  as  adjusting  its  phase  to 
the  new  jammer.  Readaptation  is  nearly  complete  after  about  130  iterations 
which  is  consistent  with  the  calculation  for  the  time  constant. 

Figure  12  shows  how  the  notch  changes  its  shape  during  the 
readaptation  phase.  Just  before  the  frequency  shift,  we  have  a  well  formed 
notch  at  the  initial  frequency.  Within  five  iterations  of  the  change  in 
frequency,  we  find  that  the  notch  depth  is  about  -10  dB.  The  notch  at  the 
lower  frequency  has  not  yet  begun  to  form  because  5  iterations  is  less  than 
one-half  of  a  cycle  of  the  new  frequency.  After  50  iterations,  a  notch  is 
beginning  to  form  at  the  new  frequency;  its  depth  is  about  the  same  as  that 
at  the  initial  frequency.  After  125  iterations,  a  notch  is  reasonably  well 
formed  at  the  „ew  frequency  and,  as  the  adaptation  process  continues,  the 
notch  depth  will  increase  while  the  remnant  notch  at  the  initial  frequency 
will  continue  to  diminish. 

The  final  example  is  that  of  suppressing  two  unequal  jammers  at 
different  frequencies.  In  this  case  one  jammer  is  at  +20  dB  relative  to  the 
signal  and  the  second  is  at  +10  dB.  Figure  13  shows  the  two  notches  formed 
by  a  48  tap,  4000  sample  accumulator  setup.  The  notch  at  the  lower 


frequency  corresponds  to  the  +10  dO  jammer;  we  see  that  the  notch  is  not  as 
deep  (  »  -20  dB)  as  is  predicted  by  (13)  and  there  is  more  evidence  of  jammer 
power  in  |Z(u>)l  at  this  frequency.  Examination  of  z(t)  shows  that  the 
adaptation  constant  (t  ^  =  175)  is  primarily  due  to  the  strong  jammer,  with 
a  secondary  adaptation  constant  due  to  the  weak  jammer.  The  correlation 
coefficient  is  p  =  0.93. 

4.  SUMMARY  AND  CONCLUSIONS 

The  purpose  of  the  simulations  was  to  determine  the  effects  of  using 
finite  accumulators  as  dictated  by  the  finite  time-bandwidth  product  of  the 
Bragg  cells.  The  measure  of  performance  used  throughout  has  been  the 
correlation  coefficient  between  the  output  signal  and  the  transmitted 
signal.  The  effect  of  using  a  finite  length  accumulator  is  to  reduce  the 
correlation  coefficient  by  =0.5  to  1.0  dB  depending  on  the  system 
parameters.  Some  degree  of  taper  is  needed  on  the  accumulator  to  reduce  the 
effects  of  readaptation.  The  amount  of  taper  required  increases  as  the 
frequency  of  the  jammer  increases.  The  natural  attenuation,  the  frequency 
response  and  the  method  of  illumination  of  the  Bragg  cells  are  possible  ways 
to  implement  the  taper.  These  methods  have  the  advantage  of  increasing  the 
apparent  time-bandwidth  product  of  the  Bragg  cell  which  help  to  improve  the 
overall  performance  of  the  system. 

The  number  of  taps  available  is  probably  greater  than  that  needed  in 
many  applications.  We  found  that  increasing  the  number  of  taps  improves  the 
performance  but  also  increases  the  time  necessary  for  the  system  to  adapt  to 
changes  in  the  jammer  frequency  or  amplitude.  In  contrast  to  a  conventional 
system  where  the  feedback  gain  can  be  set  at  a  low  value  (at  the  expense  of 
slower  adaptation),  the  gain  in  the  optical  system  must  not  be  set  too  high 
because  then  the  system  will  not  be  stable,  or  set  too  low  because  then  the 


tap  weights  will  not  be  large  enough  to  remove  the  jammer  energy.  The  not 
depth  is  not  as  great  as  when  using  infinite  accumulators  because  some  of 
the  jammer  energy  is  needed  in  the  residual  signal  to  maintain  the  tap 
weights.  The  null  is  sufficient,  however,  to  improve  the  s igna 1-to-j ammer 
ratio  at  the  output  so  that  a  high  value  of  the  correlation  coefficient  is 
obtained. 
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simulation  studies,  and  E.F.  Smith  and  R.V.  Boyd  for  their  insights  and 
suggestions.  This  work  was  supported  by  the  U.S.  Army  Research  office. 
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